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Abstract 

We give a review of modern theoretical understanding of the physics of QCD at finite 
temperature. Three temperature regions are studied in details. When the temperature is 
low, the system presents a rarefied pion gas. Its thermodynamic and kinetic properties are 
adequately described by chiral perturbation theory. 

When the temperature is increased, other than pion degrees of freedom are excited, the 
interaction between the particles in the heat bath becomes strong, and the chiral theory is not 
applicable anymore. At some point T = Tc a phase transition is believed to occur. The physics 
of the transitional region is discussed in details. The dynamics and the very existence of this 
phase transition strongly depends on the nature of the gauge group, the number of light quark 
flavors, and on the value of quark masses. If the quarks are very heavy, the order parameter 
associated with the phase transition is the correlator of Polyakov loops < P*(x)P(0) >t related 
to the static potential between heavy colored sources. When the quark masses are small, the 
proper order parameter is the quark condensate < qq >t and the phase transition is associated 
with the restoration of the chiral symmetry. Its dynamics is best understood in the framework 
of the instanton liquid model. Theoretical estimates and some numerical lattice measurements 
indicate that the phase transition probably does not occur for the experimentally observed 
values of quark masses . We have instead a very sharp crossover, "almost" a second order 
phase transition. 

At high temperatures T » fJ-hadr the system is adequately described in terms of quark 
and gluon degrees of freedom and presents the quark-gluon plasma (QGP). Static and kinetic 
properties of QGP are discussed. A particular attention is payed to the problem of physical 
observability, i.e. the physical meaningfulness of various characteristics of QGP discussed in 
the literature. 



^Permanent address: ITEP, B. Cheremushkinskaya 25, Moscow 117259, Russia. 
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1 Introduction. 

The properties of QCD medium at finite temperature have been the subject of intense study 
during the last 15-20 years. It was reahzed that the properties of the medium undergo a 
drastic change as the temperature increases. At low temperatures, the system presents a gas 
of colorless hadron states: the eigenstates of the QCD hamiltonian at zero temperature. When 
the temperature is small, this gas is composed mainly of pions — other mesons and baryons 
have higher mass and their admixture in the medium is exponentially small ~ exp{ — M/T}. 
At small temperature, also the pion density is small — the gas is rarefied and pions practically 
do not interact with each other. 

However, when the temperature increases, pion density grows, the interaction becomes 
strong, and also other strongly interacting hadrons appear in the medium. For temperatures 
of order T ~ 150 Mev and higher, the interaction is so strong that the hadron states do not 
present a convenient basis to describes the properties of the medium anymore, and no analytic 
calculation is possible. 

On the other hand, when the temperature is very high, much higher than the characteristic 
hadron scale ^ihadr ~ 0.5 Gev, theoretical analysis becomes possible again. Only in this range, 
the proper basis are not hadron states but quarks and gluons — the elementary fields entering 
the QCD lagrangian. One can say that, at high temperatures, hadrons get "ionized" to 
their basic compounds. In the 0— approximation, the system presents the heat bath of freely 
propagating colored particles. For sure, quarks and gluons interact with each other, but at 
high temperatures the effective coupling constant is small ois{T) <^ 1 and the effects due to 
interaction can be taken into account perturbatively 0. This interaction has the long-distance 
Coulomb nature, and the properties of the system are in many respects very similar to the 
properties of the usual non-relativistic plasma involving charged particles with weak Coulomb 
interaction. The only difference is that quarks and gluons carry not the electric, but color 
charge. Hence the name: Quark- Gluon Plasma [QGP). 

^We state right now, not to astonish the experts, that there are limits of applicability of perturbation theory even 
at very high temperatures, and we are going to discuss them later on. 
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Thus the properties of the system at low and at high temperatures have nothing in common. 
A natural question arises: What is the nature of the transition from low-temperature hadron 
gas to high- temperature quark-gluon plasma ? Is it a phase transition ? If yes, what is its 
order ? I want to emphasize that this question is highly non-trivial. A drastic change in the 
properties of the system in a certain temperature range does not guarantee the presence of 
the phase transition point where free energy of the system or its specific heat is discontinuous. 
Recall that there is no phase transition between ordinary gas and ordinary plasma. 

Whether or not the phase transition occurs in the real QCD with particular values of 
quark masses is the question under discussion now. Our personal feeling is that the answer is 
probably negative and what really happens is not the phase transition but a sharp crossover 
— "almost" a second-order phase transition. Another discussed possibility is that a weak first 
order phase transition (with small latent heat) still occurs. A point of consensus now is that 
the real phase transition does occur in some relative theories — in pure Yang- Mills theory 
(when the quark masses are sent to infinity) and in QCD with 2 or 3 exactly massless quark 
flavors. 

There are at least 4 reasons why this question is interesting to study: 

1. It is just an amusing theoretical question. 

2. Theoretic conclusions can be checked in lattice numerical experiments. Scores of papers 
devoted to lattice study of thermal properties of QCD have been published. 

3. Perhaps, a direct experimental study would be possible on RHIC — high-energy ion 
collider which is now under construction. 

4. During the first second of its evolution, our Universe passed through the stage of high- 
T quark-gluon plasma which later cooled down to hadron gas (and eventually to dust 
and stars, of course). It is essential to understand whether the phase transition did 
occur at that time. A strong first-order phase transition would lead to observable effects. 
In particular, it could have created inhomogeneities in baryon number density in early 
Universe which would later affected nucleosynthesis and therefore leave a signature in 
the primordial nuclear abundances [0. We know (or almost know — the discussion of 
this question has not yet completely died away) that there were no such transition. But 
it is important to understand why. 

Note that there is also a related but different question — what are the properties of rela- 
tively cold but very dense matter and whether there is a phase transition when the chemical 
potential corresponding to the baryon charge rather than the temperature is increased. This 
review will be devoted exclusively to the thermal properties of QCD, and we shall assume zero 
baryon charge density. 

At present, we are able to describe the physics of thermal phase transition mainly in 
qualitative terms (we will address, however, the issue of critical exponents and the critical 
amplitudes where certain theoretical predictions can be made). As was already mentioned. 
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there are, however, two regions where many quantitative theoretical results can be obtained. 
These are the small temperature region where the system presents a weakly interacting rarefied 
pion gas and the high temperature region where the system presents a weakly interacting (and 
also rarefied in some sense) quark-gluon plasma. 

In both cases, we are facing a reasonably clean and rather interesting problem of theoretical 
physics lying on cross-roads between the relativistic field theory and condensed matter physics. 
Personally, I had a great fun studying it. Unfortunately, at present only a theoretical study 
of the problem is possible. That concerns especially the high temperature phase. Hot hadron 
medium with the temperature above phase transition can be produced for tiny fractions of a 
moment in heavy ion collisions but: 

• It is not clear at all whether a real thermal equilibrium is achieved. 

• A hot system created in the collision of heavy nuclei rapidly expands and cools down 
emitting pions and other particles. It is not possible to probe the properties of the 
system directly, but only indirectly via the characteristics of the final hadron state. 

• Anyway, the temperature achieved at existing accelerators is not high enough for the 
perturbation theory to work and there are no quantitative theoretical predictions with 
which experimental data can be compared. RHIC would be somewhat better in this 
respect, but even there the maximal temperature one could expect to reach is of order 
0.5 GeV This is above the expected phase transition temperature ~ 200 MeV (see 
a detailed discussion in Chapters 4 and 5), but is still probably not high enough for 
the effective coupling constant to be small so that perturbative calculations would be 
justified. 

Thus at present, there are no experimental tests of non-trivial theoretical predictions for 
QGP properties. The effects observed in experiment such as the famous J/'0 - suppression 
(for a recent review see [^) just indicate that a hot and dense medium is created but says little 
on whether it is QGP or something else. 

The absence of proper feedback between theory and experiment is a sad and unfortunate 
reality of our time: generally, what is interesting theoretically is not possible to measure and 
what is possible to measure is not interesting theoretically |^. 

The situation is somewhat better with numerical lattice experiments. One of the advantages 
of the lattice approach is that the question of whether thermal equilibrium is achieved does not 
arise — path integrals on an Euclidean cylinder with imaginary time size (3 = 1/T describe the 
partition function and other characteristics of the equilibrated system just by definition. A lot 
of numerical results at the temperatures below or slightly above critical temperature exist and 
in some cases they provide very useful insights which help to understand better the physics of 
the phase transition. In principle, one can also measure on the lattice the characteristics of the 

■^One of the possible exceptions of this general rule in the field of thermal QCD is a fascinating perspective to 
observe the phenomenon of disoriented chiral condensate at RHIC and we will discuss it. 
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high temperature phase. Due to a finite size of the lattice, it is more difficult, however, and 
such studies just started. A limitation here is that we can calculate numerically only Euclidean 
path integrals which means that we can study only static characteristics of the system. Kinetic 
properties (spectrum of collective excitations in QCD heat bath, the transport phenomena such 
as viscosity and electric conductivity, etc) cannot be probed in this way. (May be, it is possible, 
however, to study a nontrivial kinetic characteristics, the rate of axial charge non-conservation 
in high temperature phase, studying numerically classical Yang-Mills theory on a hamiltonian 
lattice. We will address this question at the end of the review.) 

Thus we will not attempt (or almost will not attempt) to establish relation of the results 
of theoretical calculations with realistic accelerator experiments. Comparison with the lattice 
data will be done when the latter are available, i.e. in the low and moderate temperature 
region. What we will do, however, is discussing the relation of theoretical results with gedanken 
experiments. Suppose, we have a thermos bottle with hot hadron matter or with QGP on a 
laboratory table and are studying it from any possible experimental angle. We call a quantity 
physical if it can in principle be measured in such a study and non-physical otherwise. We 
shall see later that many quantities discussed by theorists may be called physical only with 
serious reservations, and some are not physical at all. 

Many reviews on the subject have already been written (see e.g. |||], Q - lHH)- However, 
all of them were either written a considerable time ago and do not cover very important 
recent developments or have a limited scope. I will try here to fill up this gap and present 
a modern discussion which covers a reasonably broad set of questions. Of course, I could 
not cover everything. The issues which interested me more and which I am more familiar 
with are discussed at greater length and in greater details. In cases when the question is still 
controversial and no general consensus exists, I will present my own viewpoint. A bias of this 
kind is, however, inevitable. 

I had some problems in arranging the material and found no way not to break the causality: 
cross-references both back in text and forward in text will be abundant. The physics of all 
finite temperature systems has something in common, and the parallels will be drawn all the 
time. Moreover, some basic facts about quark-gluon plasma (the subject of Chapter 6) will 
be used in Chapter 3 where the physics of pure glue systems is discussed, etc. 

But, anyway, some plan should be chosen, and it is the following. The next chapter 
presents technical preliminaries. I will give a review of finite T diagram technique. A particular 
emphasis is given to the real time formalism which is most convenient when studying kinetic 
properties of the system. It is much less known than Euclidean Matsubara technique (well 
suited to study static characteristics) and has entered the standard tool kit of field thenrists 
studying thermal field theories only recently. 

Chapter 3 is devoted to physics of pure Yang-Mills theory at finite T. The system displays 
deconfinement phase transition in temperature associated with changing the behavior of the 
potential between static quark sources. In particular, we discuss the structure of the high 
temperature phase and show that there is only one such phase, not Nc different phases with 
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domain walls between them as people believed for a long time and some continue to believe 
up to now. There is no such physical phenomenon as spontaneous breaking of symmetry 
at high temperature in the pure Yang-Mills system. 

In Chapter 4 we go over to the theories involving light dynamical quarks and discuss first 
the properties of low temperature phase. We discuss at length the temperature dependence 
of chiral condensate (it decreases with temperature) and also explore the fate of pions and 
massive hadron states at small nonzero temperature. We show that the leading effect of the 
pion heat bath on massive hadrons is that the latter acquire finite width. Real part of the poles 
of corresponding Green functions is also somewhat shifted. For nucleons, this shift is rather 
tiny. For vector mesons, the shift is also tiny in the temperature region up to ~ 100 MeV. 
What happens at higher temperatures is not quite clear by now and is a subject of intense 
discussions. 

Chapter 5 is devoted to the physics of thermal phase transition in QCD. It is associated 
with restoration of chiral symmetry which is broken spontaneously at zero temperature. We 
show how the physics of the phase transition depends on the number of light quark flavors 
(For two massless flavors it is the second order while for three flavors it is the first order. We 
argue in particular that, in the theory with 4 or may be 5 massless flavors, chiral symmetry 
would probably not broken at all and, correspondingly, there would be no phase transition ). 
We discuss in details the physical mechanism of the phase transition which, in our opinion, is 
best understood in the framework of the instanton-antiinstanton liquid model. In this model, 
characteristic vacuum fields contributing in the Euclidean functional integral for the partition 
function present a collections of quasi-classical objects: instantons and antiinstantons. Each 
(anti)instanton supports fermion zero modes (one fermion and one antifermion zero mode 
for each flavor). When interaction between quasi-particles is taken into account, the modes 
shift from zero, but, as the interaction turns OTit to be strong enough, the characteristic 
eigenfunctions become delocalized and characteristic distance between neighboring eigenvalues 
is very small ~ i^yEud^ This brings about nonzero fermion condensate. When temperature 
is increased, the density of quasi- particles decreases and, which is even more important, their 
interaction decreases. As a result, instantons and antiinstantons tend to form "molecules" with 
localized fermion eigenfunctions. Characteristic eigenvalues are far from zero in this case and 
the chiral condensate is zero. Remarkably, it is basically the same mechanism which brings 
about the so called "percolation phase transition in doped semiconductors". When the density 
of impurity is high, electrons can jump between adjacent impurity atoms, acquire mobility, 
and the material becomes a conductor. In the end of the chapter we discuss the phenomenon 
of disoriented chiral condensate. 

Physics of QGP is discussed in Chapter 6. We discuss the limits of applicability of pertur- 
bation theory due to so called " magnetic screening" - the effect which is specific for nonabelian 
gauge theories. We discuss also the spectrum of plasmons and plasminos — the collective ex- 
citations with quantum numbers of quarks and gluons in QGP. We address the controversial 
issue of the plasmon damping. We show that damping depends on the gauge convention and 
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has as such Httle physical meaning. On the contrary, transport phenomena (such as viscosity, 
electric conductivity, energy losses of a heavy energetic particle passing through QGP, etc.) 
are quite physical and can be measured in a gedanken experiment. Finally, we discuss a pure 
non-perturbative effect of chirality non-conservation in quark-gluon plasma (or, if you will 
baryon number non-conservation in hot electroweak plasma). The rate of non-conservation 
is proportional to times a power of coupling constant. What is this power and whether 
a numerical algorithm exists where this quantity can be determined is still a question under 
discussion now. 



The main point of interest for us are the physical phenomena in hot QCD system. However, as 
the main theoretical tool to study them is the perturbation theory and we want in some cases 
not only to quote the results, but also to explain how they are obtained, we are in a position 
to spell out how the perturbative calculations at finite T are performed. 

There are two ways of doing this — in imaginary or in real time. These techniques are 
completely equivalent and which one to use is mainly a matter of taste. Generally, however, 
the Euclidean technique is more handy when one is interested in pure static properties of the 
system (thermodynamic properties and static correlators) where no real time dependence is 
involved. On the other hand, when one is interested in kinetic properties (spectrum of collective 
excitations, transport phenomena, etc.), it is much more convenient to calculate directly in 
real time. 



we will describe it here only briefly. Consider a theory of real scalar field described by the 
hamiltonian H[(j)(x), n(x)]. The partition function of this theory at temperature T can be 
written as 



2 Finite T diagram technique. 



2.1 Euclidean (Matsubara) technique. 

Many good reviews of Matsubara technique are available in the literature (see e.g. Q) and 




X 




where /C is the quantum evolution operator in the imaginary time (3 = 1/T : 





n 



are the eigenstates of the hamiltonian. One can express the integral in RHS of Eq.(2.1) as 



an Euclidean path integral: 




(2.3) 
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where the periodic boundary conditions are imposed 



(/<(x,r + /?) = ,^(x,r) (2.4) 
A thermal average < O >x of any operator 0[i?i(x, r)] has the form p 



<0 >T= Z-^ j \{d(l){-^,T)0[(t){^,T)]e^^\^- j^T J (ix£[(/)(x,r]| (2.5) 

One can develop now the diagram technique in a usual way. The only difference with the 
zero temperature case is that the Euclidean frequencies of the field (/>(x, r) are now quantized 
due to periodic boundary conditions (^): 

= 2mnT (2.6) 

with integer n. To calculate something, one should draw the same graphs as at zero tem- 
perature and go over into Euclidean space where the integrals over Euclidean frequencies are 
substituted by sums: 

'^'P-J{p) ^TJ2[ ^/(2vrmr,p) (2.7) 



(2^)4— Z^J (2^)3 

The same recipe holds in any theory involving bosonic fields. In theories with fermions, one 



should impose antiperiodic boundary conditions on the fermion fields ■0(x, r) (see e.g. |13] for 
detailed pedagogical explanations), and the frequencies are quantized to 

= iiT{2n + l)T (2.8) 

An important heuristic remark is that, when the temperature is very high, in many cases only 
the bosons with zero Matsubara frequencies po = contribute in < O >t- The contribution 
of higher Matsubara frequencies and also the contribution of fermions in < O >r become 
irrelevant and, effectively, we are dealing with a 3- dimensional theory. As was just mentioned, 
it is true in many, but not in all cases. For example, it makes no sense to neglect fermion 
fields when one is interested in the properties of collective excitations with fermion quantum 
numbers. For any particular problem of interest a special study is required. 

2.2 Real time (Keldysh) technique. 

Matsubara technique is well suited to find thermal averages of static operators 0{x). If we 
are interested in a time-dependent quantity, there are two options: i) To find first the thermal 
average < 0{x,t) > for Euclidean r and perform then an analytic continuation onto the real 
time axis. It is possible, but quite often rather cumbersome, ii) To work in the real time right 
from the beginning. The corresponding technique was first developed in little known papers 
1 14] and independently by Keldysh p^] who applied it to condensed matter problems with a 



9 



particular emphasis on the systems out of equihbrium. For systems at thermal equilibrium, it 
was effectively reinvented in a slightly different approach in the seventies by the name "thermo 
field dynamics" ||l6|. It was fully apprehended by experts in relativistic field theory only in 
the beginning of the nineties. 

There is no good review on real time diagram technique addressed to field theorists. The 
existing review [0] is very deep and extensive, but is written rather formally and is hard to 
read. What we will do here is in a sense complementary to the review [0]. We will not derive 
the real time technique in a quite accurate and regular way, but rather elucidate its physical 
foundations, show how multicomponent Green's functions appear, formulate and discuss the 
real time Feynman rules, and discuss also a simplistic but rather useful version of the real time 



technique due to Dolan and Jackiw |18] 



Consider a quantum mechanical system in thermal equilibrium. Suppose a thermal av- 
erage ( ^^ ) of a Heisenberg operator ^o(*) is zero. Let us perturb the system adding to the 
hamiltonian the term 

V{t) = VoBo6{t) (2.9) 

where Bq is some operator and we assume the constant Vq to be small. The Heisenberg 
operator Anit) of the full hamiltonian H = Hq + V{t) is related to the Heisenberg operator 
Ao(t) of unperturbed hamiltonian (or the operator in interaction representation) as 

Anit) = expji j'^ y(T)(ir| Ao(t) exp |-i J*^ V{T)dT^ (2.10) 

(Generally, T-ordered exponentials enter, but for an instantaneous perturbation T-ordering 
is irrelevant). In the first order in Vq 

Anit) = Ao{t)-iVoe{t)[Ao{t),Bo] (2.11) 

The thermal average is 

< Anit) >T = -iVo < e{t)[Ao{t),Bo] >t (2.12) 



The expression in the RHS of Eq. ( p. 12 ) is called the response function — it determines the 



response of the system at some time i > at the instantaneous small perturbation applied at 
t = 0. Let us make a Fourier transform and define the generalized susceptibility 



/oo 
e*^*0(t) < [Ao{t),Bo] >Tdt 

(2.13) 



-E„/T 



AnmB'mn BnmAmn 



UJ + En - Em + iO UJ + Em - En + iO 



We see that D'^^{uj) is analytic in the upper half-plane. It is just the corollary of the fact 
that D'^^{t < 0) = 0. The relation ( p. 131) is well known in statistical mechanics and is called 
the Kubo formula. 
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We have changed the sign convention compared to the original Kubo convention (see e.g. |1S[|) to make a 



generalization to quantum field theory more transparent. 
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Consider now a thermal field theory. Let it be first a theory of real scalar field with possible 
nonlinear interactions. Suppose the perturbation is coupled to (/)(x, 0) and we measure the 
response of the system at some later time t in terms of < (p{x,t) >t- The corresponding 
response function is called the retarded Green's function 

DR{t,x) = -ie{t) < [0(t,x),0(O)] >T (2.14) 

Its Fourier image Dji{uj,p) is analytic in the upper lo half-plane. One can show that the 
retarded Green's function describes a natural analytic continuation of the Matsubara Green's 
function (defined at a discrete set of points on the imaginary a; axis) on the complex co plane. 
We will be mainly interested with real to. 
On the tree level 

D%{u^^P) = , ^..,2^ 2 2 (2-15) 

(uj + ^0j^ — P rn 

Let us separate now the free hamiltonian from the interaction part and build up the pertur- 
bation theory for Dji(t,x) and other physically observable quantities in interaction represen- 
tation. It is very well known that at zero temperature the retarded propag ators like I^J^) are 



not quite convenient for this purpose. The matter is that to calculate, say, the exact retarded 
Green's function, we have to draw the loops involving virtual particles which do not necessarily 
go forward in time. In other words, exact retarded Green's function cannot be expressed into 
integrals of only retarded propagators. Advanced components also play a role and, as a result, 
we are arriving at the "old diagram technique" which is somewhat clumsy. 

Feynman showed that the calculations can be greatly simplified using T - ordered Green's 
functions. In that case, everything is expressed into the integrals of T - ordered propagators. 
We know that the Feynman technique for field theories at zero temperature is explicitly Lorentz 
- invariant. However, at finite temperature, Lorentz-invariance is lost — the reference frame 
where thermal medium is globally at rest is singled out. That gives an indication that Feynman 
trick probably would not work at finite temperature and we are bound to use a version of the 
old diagram technique which is not Lorentz-invariant. 

Another indication of the trouble is that the notion of S- -matrix which is basic at T = 
and for which the Feynman diagram technique has been constructed has absolutely no meaning 
at finite temperature (see for a formal proof). Particles interact with the heat bath all 
the time, scatter on the real particles there and just have no chance to arrive from infinity or 



to escape to infinity from interaction point (see [21| for a related physical discussion). 

What can be directly measured in the heat bath are classical response functions (see a 
detailed discussion of this issue in Chapter 6). Suppose we disturb normal plasma (or pion 
gas or QGP) with a concentrated laser beam at t = and measure the distribution of charges 
and corresponding electromagnetic fields at later time. Such a gedanken experiment is quite 
feasible. Thus retarded Green's functions (in contrast to the Feynman Green's functions) have 
a direct physical meaning and we should be able to calculate them. 
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Let us show that the Feynman program fails indeed at T 7^ Consider the exact 

Feynman propagator Df{x) = —i < T{(j){x)(j){0)} >t- Heisenberg field operators are 



(x,i) = S'(-cxD,t)(/>o(x,t)5(t,-cxD) (2.16) 



where (/)o(x, t) are the operators in the interaction representation, 

S{t2,ti) = Texp!^-i J^^\{t')dt'Y (2.17) 
and V{t) here is a nonlinear part of the hamiltonian and has nothing to do with the instanta- 



neous perturbation (|2.9D . Simple transformations give (see e.g. p3| ) 

Df{x) = -i < S-^T{M^)MO)S} >T (2.18) 

where S = S{oo, —00) is the S- matrix operator. At zero temperature we were interested 
with the vacuum average, S'-matrix gave just a phase factor when acting on the vacuum state 
(this factor was anyway cancelled out with the vacuum loops in perturbative expansion of 
T{(j){x)(j){0)S}) and, expanding S in powers of V , using Wick theorem etc., we obtained a 
standard Feynman loop expansion for the exact propagator. 

Thermal average (2.5) involves, however, also averaging over excited states on which the 
operator S acts in a non-trivial way. We are bound to take the operator 

S-^ = Texpli r V{t')dt'] (2.19) 



into account when performing Wick contraction of the field operators (T stands for anti- 
chronological ordering where the operators at later times stand on the right) 
As a result, three types of Green's functions appear: 

• The usual T - ordered ones coming from the contractions inside T{<j)Q{x)(j)Q{0)S}] 

• Anti-T-ordered Green's functions coming from the expansion of S~^; 

• Finally, there are also not ordered Green's functions < (f)Q{x)(pQ{0) >t coming from 
contracting the operators in S^^ with the operators in T{(j)Q{x)(j)Q{0)S}. 

Let us first deal with the latter and calculate 

-i<Mx)MO)>T= -iJ2e-^^- <n\Mx)MO)\n> (2.20) 

n 

where \n > present the eigenstates of unperturbed hamiltonian so that En = Z^p-^p^p 
with integer {A'^p} , and Cp = \/p^ + rr?. Introduce as usual a finite spatial volume V and 
decompose 

Mx) = E ke-p*+^P'' + a+e-p*-H (2.21) 

p V^^p*^ 
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where and Op are the creation and annihilation operators. Substitute it in ( p.20| ). At 
zero temperature only the vacuum average contributes and we can use the fact ap|0 >= 
to obtain a usual expression. At finite temperature, the excited states |n > with ap\n >^ 
contribute in the sum (2.20), and an additional contribution arises: < n\apap\n > = Np, < 
n\apap\n > = 1 + A^p. Trading as usual the sum for the integral and taking Fourier image , 
we obtain 



< 



ITT 



(1 + nBiep))5{uj - ep) + nB{ep)6{uj + ep)] 



(2.22) 



where 



nsiep) = <Np>T 



g/3ep 



(2.23) 



is the Bose distribution function. The indices 21 appear as a reminder that the operator 
standing on the left in in ( p. 20 ) and corresponding to the "final point" appeared from the factor 
S^^ in ( 2.18| ) while the operator standing on the right — from the factor T{(j)Q(x)(j)Q(0)S}. 
Similarly, a function 



D 



2H 



-UJ, 



Dl 



(2.24) 



can be defined. 

T - ordered and anti-T - ordered tree propagators will also differ from their zero temper- 
ature expressions. Proceeding in the same way, we obtain 



Z)?i(^,p) 



^22 (^^,P) 



2TTinB{^p)5{oj'^ 



ipx 



- e2 + iO 



2t: in B{ep) 5 {uj^ — 



a;2 - el 



iO 



(2.25) 



Various kinds of thermal propagators in (2.22, 2.24, 2.25| ) are conveniently "organized" in a 
2x2 matrix: 



D\^{LO,p),DUu,jy) 

Z)0l(^,P),^2°2(^,P) 



(2.26) 



To derive diagram technique, suppose first that the interaction V{t) corresponds to scattering 
on a static classical field -F(x): 



V{t) 



(l)'^{x,t)F{x)dx 



(2.27) 
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Let us consider first the Feynman propagator (2.18). Expand both S and S^^ in V{t) and 
compose various products of field operators contracting them according to Wick rules, (pix) 
(corresponding to the outgoing external leg) can be paired with an operator from the expansion 
of S or with the operator from the expansion of S^^. That gives us either the 11 - component 
of the Green's function ( |2.26| ) or the 12- component. Suppose the pairing occurred within 
S~^. The remaining field operator in ( |2.27| ) can be further paired either with an operator in 
S or with an operator in S~^. In the first case, we get 21 - component and in the second — 
22 - component. The last pairing occurs with 0(0) and we have either 21 or 11 - component 
depending on where the line came from. 
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Figure 1: Keldysh propagator in external field. 



The resultive "wandering" of the virtual particle is depicted in Fig. |^ ( Note that two hnes 
in Fig. H which, in our approach, correspond to the structures T{(p{x)4>{0)S} (the lower line) 
and (the upper line) have a direct correspondence in the path integral approach [p4| , |T7| 
which we will not discuss here.). When writing a corresponding analytical expression, we have 
to have in mind that the interaction vertex on the upper line corresponding to has an 
opposite sign compared to usual one — just because the sign of V{t) in S^^ is reversed. 

Likewise, we can calculate in any order the 22 - component of the exact Green's function in 
which case the virtual particle line starts and ends in the anti-chronological-ordered domain. 
For the mixed component component D21, the wandering starts on the lower line in Fig. |l| 
and ends up on the upper line. For the mixed component component D12 — the other way 
round. Two last statements follow from the simple fact 

< 0(x)0(O) > = < S{-oo,t)(f)o{x)S{t,-oo)S{-oo,0)(po{0)S{0,-oo) > = 
< S{-oo, t)Mx)S{t, 00)5(00, O)0o(O)5(O, -00) > = < f{Mx)S~^}T{MO)S} > (2.28) 

and similarly 

< mn^) > = < f{MO)S-'}T{Mx)S} > (2.29) 

We are ready now to formulate general Feynman rules in the Keldysh technique. Suppose 
we want to find the exact thermal Green's function in the theory with interaction Xcp"^ /6 ^ at 
finite T. Present in in the matrix form like in ( 2.26| ). Then 11 - component would stand for 



thermal average of T - ordered product of exact Heisenberg operators etc. Draw the same 
graphs as at T = 0. Each line corresponds now not to a scalar function, but to a matrix ( |2.26| ). 



^Never mind that the theory does not exist due to vacuum instability, this is only an illustrative example. 
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The vertices are now tensors Tabc where a,b,c = 1,2. The tree vertices (to be substituted in 
the graphs) are, however, very simple: 



111 



A, r 



222 



-A 



(2.30) 



and ah other components are zero. T222 has the opposite sign compared to due to the 



reversed sign of V{t) in (2.19). 
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Figure 2: One-loop Green's function in (j)^ theory. 



Take the graph in Fig.^ as an example. The corresponding analytic expression reads 

Dl,{u -z)j d^z J d^y rlaD%iz - y)Di,iz - yjF",,,, D^V " ^) (2-31) 

Summation over the indices which occur twice is assumed. 

Similarly, exact Green's functions with arbitrary number of legs can be found. Generally, 
exact vertices have non-vanishing mixed components like ri2i etc. 

Look again now at ( 2.26| ). Note that not all components in the matrix Green's function 
are independent. The relation 

+ ^22 = D^i2 + ^21 = D% (2.32) 

holds. This relation is true in a general system, not necessarily in thermal equilibrium. One 
can show that the same relation holds also for the exact Green's functions. 

One can present (and it makes a lot of sense) the thermal propagator in the form which 
involves only 3 independent components. To this end, we make an orthogonal transformation 

Dn = Q-^DQ (2.33) 

where 



We have 
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where Df> = Du — D12 and Da = Dn — D21 are nothing else as retarded and advanced 
components of the Green's function. Indeed, subtracting a straight product from the T - 
product, one gets the commutator ( 2.l4 ) for Dr and the commutator 



DA{t,x) = ie{-t) < [0(t,x),0(o)] >T 



(2.36) 



for Da- For the Fourier components the relation Da{^^, p) = D'^{u, p) holds so that Da{^^, p) 
is analytic in the lower to half-plane. 

-Dp is defined in ( p. 32 ). It is symmetric in lo due to (2.24). Like Dr and Da, it also 
has a direct physical meaning. For example, the rate of photon production in quark-gluon 
plasma is defined is given by the imaginary part of the P - component of the correlator of 
electromagnetic currents. As far as the 2-point Green's function is concerned, P-component 
generally provides the information about particle distribution on the energy levels. 

For future references, let us write down explicitly the inverse relations 



Du 
D21 



1 



(Dr + Da + Dp), D 



12 



-(Dr-Da + Dp), D: 



1 



1 
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-Dr + Da + Dp), 
-Dr-Da + Dp), 



(2.37) 



2' " " 2 

The matrix Green's functions ( 2.261) , ( 2.35| ) were invented in the first place to describe the 
processes out of equilibrium. For example, studying P - component of the 2-point Green's func- 
tion in the state which is slightly out of thermal equilibrium, one can easily derive Boltzmann 



kinetic equation [15, ^]. In the general case, Dr = D\ and Dp are completely independent 
functions. 

We are interested here in the system at thermal equilibrium. Then Dp is related to Dr 
and Da- It is not difficult to check that, at the tree level, the relation 



D%{uo,p) 



coth^(D 



D^, 



(2.38) 



holds. It can be proven that the same relation holds also for the exact Green's functions. 

As the relation ( |2.32| ) holds also for the exact Green's function, the rotated exact Green's 
function Dq retains the form ( ^.35 ) with zero in the upper left corner and involving only 
physical components. 

All the steps of our derivation are easily generalized on the fermion case. The fermion 
Green's function has a matrix form. On the tree level, 



G?i(u;,p) 
= {uj^^ — p7 + m 

G^2(^,P) 

(^7° - P7 + m) 



< T{Mx)MO)} >T e^P" d'x 
+ 2t: in p{e-p) 5 {uj^ — 



el + iO 



+ 2t: in p{e-p) 5 {uJ^ — 
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ITT 



-{uY - P7 + m)— [(1 - nF{ep))6{u; + ep) - nF{ep)S{uj - Cp) 
G°i(c^,p) = -i [ < M^)MO) >T e'P^d^x = 



-{ioj° - P7 + m)— [(1 - nF{ep))5iio - ep) - nFiep)6{u; + ep)] (2.39) 

where 

npiep) = (2.40) 

is the Fermi distribution function. 

Again, four components of the matrix Green's function satisfy the relation (|2.32| ) (which 
holds also in non-equilibrium case) and the orthogonal transformation ( p.33| ) brings the Green's 
function in the triangle form 

^ ( ) 

where Gr = and Gp are the physical components. In thermal equilibrium, the relation 

Gp(c^,p) = tanh^(GH-GA) (2.42) 

holds. 

The exact matrix Green's functions satisfy Dyson equations 

D = D^ + D^HD, G = G^ + G^llG (2.43) 

where 11 and S are boson and fermion polarization operators. Writing 11 = [D^)~^ — 
D~^, Tj = {G^)~^ — G~^ and using the relation (2.32), one can easily prove that 

^n,b = 0, Y.^ab = (2.44) 

ab ah 

Also for a general multileg one-particle irreducible vertex the relation 

E ^{a.} = (2.45) 

{ai=l,2} 

holds. 



Rotating 11 and S with the matrix (2.34), one can present polarization operators in the 
triangle form 
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where Hn = Hu + Hu, IIa = Hn + 1121, Hp = ITn + 1122 and similarly for S. The 
remarkable fact is that the Dyson equations ( ^.43| ) for the retarded and advanced Green's 
function components have a simple form 

Dr = dI + D%YIrDr, Da = D\ + D\YiADA (2.47) 



and similarly for fermions. This is best seen using the triangle form (|]3|), (^ ), 



That means, in particular, that to find the poles of the exact retarded Green's functions 
which determine dispersion laws of collective excitations in thermal medium (see Chapter 6 
for more details), one need to know only retarded components of the corresponding polarization 
operators. 

There is, however, one subtle point which we are in a position to discuss. In our derivation, 
extensively used the notion of S- matrix in spite of the fact that, as was mentioned earlier, it 
does not exist generally speaking at finite T. This principal difficulty corresponds to a certain 
technical difficulty. Consider e.g. the scattering process in Fig. [l|. The amplitude constructed 
according to Keldysh rules is well defined as long as scattering at individual vertices occurs 
at nonzero angle. If the external field momentum is zero, we cannot define individual terms 
in the elements of the matrix product of adjacent Green's functions D^{p)D^{p). Really, 
say £)]'2(p)-D2i(p) involves a product of two 5- functions with the same argument and is not 
defined. 

Still, a consistent real time technique can be defined, and the recipe is the following: i) 
Let us assume that the tree retarded and advanced Green's functions involve a small but 
nonzero damping. That means that ie with small but finite e should be substituted for iO in 



( 2.15 ) etc. ii) Calculate the graphs (better in triangle basis) with such Dr and Gr assuming 



Da = D*j^, Ga = G\ and substituting for Dp, Gp the expressions ( p. 38 ) and ( 2.42| ). Hi) Send 



e to zero in the very end. One can show that the amplitudes defined in such a way involve no 
pathologies and are well defined [0]. Anyway, this recipe is very physical: a nonzero damping 
e means that no asymptotic states really exist. 

e as defined above is a technical parameter and is eventually sent to zero. However, there 
is also a physical damping due to scattering processes which displays itself as the imaginary 
part of the pole in the exact Green's functions. A natural physical estimate for damping is 
the inverse free path time Ty^ee path (this is not always so — see a detailed discussion of the 
fascinating issue of so called anomalous damping in Chapter 6 — but let us neglect these 
subtleties here). When one calculates real time Green's functions, the scale r/^ee path shows 
up as a barrier beyond which the perturbative series explodes and higher order corrections are 
not under good control. 

It suffices often to use a simplified version of the real time technique due to Dolan and 
Jackiw It amounts to neglecting all other components of the matrix Green's function 

except Dii, Gu. Generally, it is wrong — in particular, the amplitudes involving zero angle 
scattering would be pathological. But, for many simple problems, other matrix components 
do not contribute, indeed, and using Dolan-Jackiw technique is quite justified. Then the 
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extra term ~ nB{e)6{io'^ — e^) in D^^ (and similarly for fermions) reflects the presenbe of real 
particles in the heat bath with Bose (Fermi) distribution. Sometimes it is called a "temperature 
insertion" . 

Suppose e.g. we are interested in the real part of the retarded polarization operator in 
the (j)^ theory calculated by the one-loop graph in Fig.^. We have Hr = Hu + ni2. Iii2 
involves a product of two 5 - functions. That means that both virtual lines are put on mass 
shell and hence is purely imaginary and does not contribute in the real part. As for IIn 
— it depends on the one-loop level only on T- ordered components of the Green's function. 
We have an ordinary T = graph and a graph with one temperature insertion (a graph with 
temperature insertions in both lines contributes in the imaginary part). 

Such a calculation is quite meaningful. Not in (p^ theory which does not exist, of course, 
but for hot Yang-Mills theory where we also have triple interaction vertices. Real parts of 
dispersion laws of plasmon and plasmino excitations are calculated just in this way. 

3 Pure Yang-Mills theory: deconfinement phase tran- 
sition. 

3.1 Preliminary remarks. 

We will start our discussion of the physics of hot nonabelian gauge theories with the theory not 
involving dynamical quarks — a pure Yang-Mills theory. Like QCD, this theory is believed 
to be confining (there is no, of course, direct experimental evidence as in QCD, but it follows 
from theoretical arguments and from lattice measurements). Its physics is, however, rather 
different. The spectrum does not involve here mesons and baryons we are used to, but only 
glueball states. The lowest glueball state has a mass of order 1.5 GeV. It is not seen in the 
real world due to large mixing to meson states, but in the pure glue theory it should be just 
a stable particle. 

As the temperature is increased, glueball states appear in the heat bath. When the temper- 
ature is very high, the system is better described in terms of original gluon degrees of freedom. 
The interaction of gluons at very high temperatures is weak and the system presents pure 
gluon plasma. 

We will show now that there are serious reasons to believe that, when we increase the 
temperature and go over from the low temperature glueball phase to the high temperature 
gluon plasma phase, a phase transition (the deconfinement phase transition) occurs on the 
way. That was argued long ago by Polyakov [|2^ and Susskind . On the heuristic level, the 
reasoning is the following: 

At zero temperature the theory enjoys conflnement. That means that the potential between 
the test heavy quark and antiquark grows linearly at large distances: 

T = 0: VQQ{r)r.ar, r ^ oo (3.1) 
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On the other hand, at high temperature when the system presents a weakly interacting plasma 
of gluons, the behavior of the potential is quite different: 



T FqqM ~ ^e---^ (3.2) 

Here m/j ~ gT is the Debye mass (see Chapter 6 for detailed discussion), and the potential is 
the Debye screened potential much similar to the usual Debye potential between static quarks 
in non-relativistic plasma. There is no confinement at large T. There should be some point 
Tc (the critical temperature) where the large r asymptotics of the potential changes and the 
phase transition from the confinement phase to the Debye screening phase occurs. 

These simple arguments can be formulated in a rigorous way. Consider the partition 
function of the system written as the Euclidean Matsubara path integral. Gluon fields satisfy 
periodic boundary conditions in imaginary time 

A;1(x,/3) = ^^(x,0) (3.3) 

where /5 = 1/T. Let us choose a gauge where Aq is time-independent. Introduce the quantity 
called the Polyakov loop 

P(x) = i-TV exp{iff/3AS (x)tn (3.4) 
It is just a Wilson loop on the contour which winds around the cylinder. Consider the correlator 

Ct(x) =<P(x)P*(0) >T (3.5) 



One can show p7| that the correlator (|3.5[ ) is related to the free energy of the test heavy 
quark-antiquark pair immersed in the plasma. 

Ct(x) = ^ exp{-/3F(|(r)} + ^ exp{-/?F(°^(r)} (3.6) 

where r = |x|. F^^^(r) and F^Q^r) are free energies of test quark-antiquark pairs (alias static 
potentials) in the triplet and, correspondingly, the singlet net color state. Let us take now the 
limit r — > oo. The quantity 

Ct(oo) = lim Crir) (3.7) 

r — >co 

plays the role of the order parameter of the deconfinement phase transition. At small T, 
FQq{r) grow linearly at r ^ cxd and Ct'(oo) = 0. At large T, free energies do not grow 
and Ct(oo) is some nonzero constant (if one would naively substitute in Eq.( |3.6D the Debye 
form of the potentials ( |3.2D , one would get CT>/i^^^^(cxD) = 1, but it is not quite true because 



FQQ{r) involve also a constant depending on the ultraviolet cutoff of the theory. See |28, 29 
for detailed discussion). There is a phase transition in between. 
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What are the properties of this phase transition ? There are not quite rigorous but sug- 
gestive theoretical arguments based on the notion of "universahty class" which predict 
different properties for different gauge groups. The main idea is that the pure YM theory based 
on SU{2) color group has some common features with the Ising model (with global symmetry 
Z2), the theory with SU{3) gauge group — with a generalized Ising model (the Potts model) 
with the global symmetry Z3 etc. The Ising model has the second order phase transition, and 
the same should be true for pure SU{2) gauge theory. A system with Z3 symmetry display, 
however, the first order phase transition, and the same should be true for pure SU (3) theory. 
The lattice data pl[ are in a nice agreement with this prediction. Also, the critical exponents 



of the second order phase transition were measured |32|. Their numerical values are close to 
the numerical values of critical indices in the Ising model. 

The arguments presented are rather suggestive, and the deconfinement phase transition 
is probably there, indeed. However, there is no absolute clarity here yet. We will return to 
discussion of this question later on and now let us dwell on a very confusing issue of domain 
walls and bubbles in the high - T phase. 

3.2 Bubble confusion. 

There was a long-standing confusion concerning the nature of deconfinement phase transition 
in pure YM theory. It has been clarified only recently. 

In scores of papers published since 1978, it was explicitly or implicitly assumed that one 
can use the cluster decomposition for the correlator (|3.5D at large T and attribute the meaning 
to the temperature average < P >t- Under this assumption, the phase of this average can 
acquire Nc different values: < P >t= C exp{27rik / Nc} , k = 0,1, Nc — 1 which would 
correspond to Nc distinct physical phases and to the spontaneous breaking of the discrete Zjv- 
symmetry. Assuming that , the surface energy density of the domain walls separating these 
phases has been evaluated in [^]. Then the "bubbles" of one of the high temperature phases 
inside another would be abundant in early Universe. These bubbles would be surrounded by 
the walls carrying significant surface energy. That could influence the evolution of the Universe 
and could lead to some observable effects. 

However, the standard interpretation is wrong. In particular: 



1. Only the correlator ( |3.5| ) has the physical meaning. The phase of the expectation value 
< P >T is not a physically measurable quantity. There is only one physical phase in the 
hot YM system. 



2. The "walls" found in |33] should not be interpreted as physical objects living in Minkowski 
space, but rather as Euclidean field configurations, kind of "planar instantons" appearing 
due to non-trivial tti^Q) = Zj\f where Q = SU{N)/Zn is the true gauge symmetry group 
of the pure Yang-Mills system. 

3. The whole bunch of arguments which is usually applied to nonabelian gauge theories can 
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be transferred with a little change to hot QED. The latter also involves planar instantons 
appearing due to non-trivial vri [[/(!)] = Z. These instantons should not be interpreted 
as Minkowski space walls. 

We refer the reader to our paper ||2^ where a detailed argumentation of these statements 
was given. Here we restrict ourselves by outlining the main physical points in the arguments. 
A. Continuum theory. 

A preliminary remark is that the situation when the symmetry is broken at high tempera- 
tures and restores at low temperatures is very strange and unusual. The opposite is much more 
common in physics. We are aware of only one model example where spontaneous symmetry 



breaking survives and can even be induced at high temperatures |34|. But the mechanism of 
this breaking is completely different from what could possibly occur in the pure Yang-Mills 
theory. 

Speaking of the latter, we note first that there is no much sense to speak about the spon- 
taneous breaking of Zj^j - symmetry because such a symmetry is just not there in the theory. 
As was already mentioned, the true gauge group of pure YM theory is SU {N) /Z^ rather than 
SU{N). This is so because the gluon fields belong to the adjoint color representation and are 
not transformed at all under the action of the elements of the center Z^r of the gauge group 
SU{N). 

< P >T as such is not physical because it corresponds to introducing a single fundamental 
source in the system: < P >t = exp{— /^F^} where Ft is the shift in free energy of the heat 
bath where a single static fundamental source source is immersed [^]. But one cannot put a 
single fundamental source in a finite spatial box with periodic boundary conditions (Such 
a box should be introduced to regularize the theory in infrared). This is due to the Gauss law 
constraint: the total color charge of the system "source -|- gluons in the heat bath" should 
be zero, and adjoint gluons cannot screen the fundamental source. This observation resolves 
the troubling paradox: complex < P >t would mean the complex free energy Ft which is 
meaningless. 

The "states" with different < P >t could be associated with different minima of the 



effective potential [37 



I A ■ \ ji I I 

mod.2 



12 



ttT 



(3.8) 



For simplicity, we restrict ourselves here and in the following with the SU (2) case. 

This potential is periodic in with the period 27rT/g. The minima at 
correspond to P = 1 while the minima at ^ = 27r(2n -|- ^)T/g correspond to P = —1. There 
are also planar (independent of y and z) configurations which interpolate between = 
at X = —00 and Aq = 2-KT/g at x = 00. These configurations contribute to Euclidean path 
integral and are topologically nonequivalent to the trivial configuration A^ = Q (Note that the 
configuration interpolating between = and Aq = AirT/g is topologically equivalent to the 
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trivial one. Such a configuration corresponds to the equator on 5^ = SU{2) which can be 
easily slipped off. A topologically non-trivial configuration corresponds to a meridian going 
from the north pole of the sphere to its south pole and presents a noncontractible loop on 
SU{2)/Z2 ). Actually, such configurations were known for a long time by the nickname of 't 
Hooft fluxes |3^ ]. 

Minimizing the surface action density in a non-trivial topological class, we arrive at the 
configuration which is rather narrow (its width is of order {gT)~^) and has the action density 

a'^^'^ = ^ + CgT' (3.9) 
(as was shown in [29|, the constant C cannot be determined analytically in contrast to the 



claim of 1 33 1 due to infrared singularities characteristic for thermal gauge theories [^9[). These 
topologically non-trivial Euclidean configurations are quite analogous to instantons. Only 
here they are delocalized in two transverse directions and thereby the relevant topology is 
determined by 7ri[G] rather than n^lG] as was the case for usual localized instantons. But, by 
the same token as the instantons cannot be interpreted as real objects in the Minkowski space 
even if they are static (and, at high T, the instantons with the size p ^ T^^ become static), 
these planar configurations cannot be interpreted as real Minkowski space domain walls. 

I want to elucidate here the analogy between nonabelian and abelian theories. The effective 
potential for standard QED at high temperature has essentially the same form as (|3.8| ): 

I (3.10) 



12 



V vrT 



mod. 2 



It is periodic in Aq and acquires minima at = ^imT/e. Here different minima correspond 
to the same value of the standard Polyakov loop -Pi(x) = exp{ie/3Ao(x)} = 1. One can 
introduce , however, the quantity Pi/jv(x) = exp{ze/3^o(x)/A^} which corresponds to probing 
the system with a fractionally charged heavy source : egourcc = e/N. Note that a fractional 
heavy source in a system involving only the fermions with charge e plays exactly the same 
role as a fundamental heavy source in the pure YM system involving only the adjoint color 
fields. A single fractional source would distinguish between different minima of the effective 
potential. If — > (X) , all minima would be distinguished, and we would get infinitely many 
distinct "phases". 

But this is wrong. One cannot introduce a single fractional source and measure < Pi/n >t 
as such due to the Gauss law constraint — the total electric charge should be zero and 
integer-charged electrons cannot screen a fractionally charged source. What can be done 
is to introduce a pair of fractional charges with opposite signs and measure the correlator 
< -fi/Af (^)-Pi*/Ar(0) >T- The latter is a physical quantity but is not sensitive to the phase of P. 
The same concerns the correlator < Pi/j^{xi) . . .Pi/]\f{xiy) >t which corresponds to putting 
N fractional same-sign charges at different spatial points. 

Finally, one can consider the configurations Aq (x) interpolating between different minima of 
( |3.10 ). They are topologically inequivalent to trivial configurations and also have the meaning 
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of planar instantons Q But not the meaning of the wahs separating distinct physical phases. 
The profile and the surface action density of these abelian planar instantons can be found 
in the same way as it has been done in Ref.|33] for the nonabelian case. For configurations 
interpolating between adjacent minima, one gets 



a 



u(l) 



+ CeT^ ln(e) 



(3.11) 



where C is a numerical constant which can in principle be analytically evaluated. 

There is a very fruitful and instructive analogy with the Schwinger model. Schwinger 
model is the two-dimensional QED with one massless fermion. Consider this theory at high 
temperature T ^ g where g is the coupling constant (in two dimensions it carries the dimension 
of mass). The effective potential in the constant Aq background has the form which is very 
much analogous to (|3:8|,|3l0|): 



ttT^ 



ttT 



1 2 



1 



mod. 2 



(3.12) 



It consists of the segments of parabola and is periodic in with the period 2tiT /g. Different 
minima of this potential are not distinguished by a heavy integerly charged probe but could 
be distinguished by a source with fractional charge. Like in four dimensions, there are topo- 
logically non-trivial field configurations which interpolate between different minima. These 
configurations are localized (for d = 2 there are no transverse directions over which they could 
extend) and are nothing else as high-T instantons. The minimum of the effective action in the 
one-instanton sector is achieved at the configuration |29, 40 1 

^exp |-^(x - xo)} , X < xo 
- exp |^(xo - x)|j , X > xo 



Aq{x) 



ttT 

a 



(3.13) 



The instanton (3.13) is localized at distances x — xq ~ g~~^ and has the action Sj = n^^'^T/g. 
But, in spite of that it is time-independent, it is the essentially Euclidean configuration and 
should not be interpreted as a "soliton" with the mass MgoLi = TSj living in the physical 
Minkowski space. 
B. Lattice Theory. 

The most known and the most often quoted arguments in favor of the standard conclusion 
of the spontaneous breaking of ^Tv-symmetry in hot Yang-Mills theory come from lattice 
considerations. Let us discuss anew these arguments and show that, when the question is 
posed properly, the answer is different. 

Following Susskind |26|, consider the hamiltonian lattice formulation where the theory 
is defined on the 3-dimensional spatial lattice and the time is continuous. In the standard 
formulation, the dynamic variables present the unitary matrices V{r, n) dwelling on the links 



^In the abelian case, there are infinitely many topological classes: 7ri[J7(l)] = Z. 
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of the lattice (the link is described as the vector starting from the lattice node r with the 
direction n). The hamiltonian is 

^ = E - ^ E -^{ViV.V.V,} (3.14) 

links plaq. 

where a is the lattice spacing, g is the coupling constant and E"" have the meaning of canonical 
momenta [i?"(r, n), ^(r, n)] = t"y(r, n). Not all eigenstates of the hamiltonian ( p. 14] ) are, 
however, admissible but only those which satisfy the Gauss law constraint. Its lattice version 
is 

G"(r) = ^^"(r,n) = (3.15) 

n 

It is possible to rewrite the partition function of the theory ( p. 14 , 3.15| ) in terms of the dual 
variables Or G SU{2) which are defined not at links but at the nodes of the lattice, 
are canonically conjugate to the Gauss law constraints ( 3.15| ) and have the meaning of the 
gauge transformation matrices acting on the dynamic variables ^^(r, n). fir correspond to the 
Polyakov loop operators ( |3.4[) in the continuum theory. In the strong coupling limit when the 
temperature is much greater than the ultraviolet cutoff A^itr ~ l/flj the problem can be solved 
analytically. The effective dual hamiltonian has 2 sharp minima at = 1 and ^Ij. = —1 and 
this has been interpreted as the spontaneous breaking of Z2-symmetry. 

Note, however, that the same arguments could be repeated in a much simpler and the 
very well known two-dimensional Ising model. Being formulated in terms of the physical spin 
variables a, the theory exhibits the spontaneous breaking of Z2-symmetry at low temperatures, 
and at high T the symmetry is restored. But the partition function of the Ising model can 
also be written in terms of the dual variables rj defined at the plaquette centers |41|. 

Z = Y.exp{-(3H{{ai})} = A{(3) exp{- [3* H'^^'^^ {{rj,})} (3.16) 

where (Ji are the original spin variables sitting on the nodes of the two-dimensional lattice, rji 
are disorder variables residing on the plaquettes, f3* = — l/21n(tanh f3) and the dual hamilto- 
nian H'^^°'\{r]i}) = —J2ij{VijVi+i,j + VijViJ+i) has the same functional form as the original 
one. When f3 is large, fi* is small and vice versa so that dual variables are ordered at high rather 
than at low temperatures and, accepting the dual hamiltonian /^'^""'({r/jj}) at face value, one 
could conclude that spontaneous breaking of Z2 symmetry occurs at high temperatures. This 
obvious paradox is resolved by noting that the dual variables are not measurable and have 
no direct physical meaning. The "domain wall" configurations interpolating between rj = 1 
and r] = —1 do contribute in the partition function formulated in dual terms. But one cannot 
feel these configurations in any physical experiment. 

Another example of a lattice theory which is even more similar to the lattice pure YM 
theory in question is Z2 three-dimensional Ising flux model which we will discuss in little more 
details later. Its dual is the standard 3D Ising spin model. The dual spin variables are ordered 
at high T and, treating the dual spin hamiltonian seriously, we would have two ordered phases 
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with domain walls which separate them. However, as dual spin variables are not physical and 
cannot be observed, such a physical interpretation is wrong [^ ]. 

And the same concerns the lattice pure YM theory [29, 43, 42 1. There are configurations 
interpolating between different Jlr G Z2 and contributing to the partition function, but they 
do not correspond to any real-time object and cannot be detected as such in any physical 
experiment. That means, in particular, that metastable states in electroweak theory of the 
kind discussed in |44| also do not exist. The absence of metastable states in hot QCD [in 
the theory with quarks, the effective potential V^^^ {Aq) is no longer periodic under the shift 
j4q — > Aq+2'kT / g but involves global minima with P = 1 at = AirT/g and local metastable 
minima at Aq = 27r(2n + l)T/g with P = —1] was shown earlier in [^ ]. 

To summarize, there is only one physical phase at high T. Its properties are relatively 
simple — it is the weakly interacting plasma of gluons. The description in terms of dual 
variables is useful for some purposes but one should be very careful not to read out in it 
something which is not in Nature. 



3.3 More on phase transition. 

Bubbles are not there, but is the deconfinement phase transition there ? If it is not really 
associated with spontaneous breaking of physical symmetry, what is its physics and are 
the heuristic arguments presented in the beginning of this chapter really compelling ? We will 
discuss here this question from different viewpoints giving the arguments pro and contra. 

An additional argument pro comes from the large Nc analysis. In the plasma phase, the 
energy density E{T) is proportional to (A'^^ — 1)T^ where T'^ appears by dimensional reasons 
and the factor — 1 is (half) the number of colored degrees of freedom. 

We see that the energy density becomes infinite in the limit Nc ^ 00. That means that if 
we start to heat the system from T = 0, we just cannot reach the QGP state — to this end, 
an infinite energy should be supplied ! 

This physical conclusion can also be reached if looking at the problem from the low tem- 
perature end. The mass spectrum of the theory in the limit Nc ^ 00 involves infinitely many 

^The exact calculation for the coefficient of the leading term and the perturbative corrections is presented in 
Chapter 6. 
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narrow states. The density of states grows exponentially with mass |^ 

p{M) oc e^^ (3.17) 

The contribution of each massive state in free energy density is exponentially suppressed 
cx exp{— M/T} for large M/T [cf. Eq. ( [4.21| ) in the next chapter]. But the total free energy 



j p{M)e-^'/^dM (3.18) 



diverges at T > c ^ . There is a Hagedorn temperature Th = c ^ above which a system cannot 
be heated ||. 

When Nc is large but finite, no limiting temperature exists (It is seen also from the low 
temperature viewpoint: at finite Nc the states have finite width and starting from some energy 
begin to overlap and cannot be treated as independent degrees of freedom), and one can bring 
the system to the plasma phase when supplying enough energy. But when N^. is large, the 
required energy is also large. That suggests (though does not prove, of course) that at large 
finite Nc a first order phase transition with a considerable latent heat takes place. 

This conjecture is supported by the lattice measurements |^l| . Indeed, a first order phase 
transition is observed at A'^c = 3 and a second order phase transition — at A'^c = 2. A 
first order phase transition need not be associated with an order parameter. But a second 
order phase transition usually is. The problem here that a good local order parameter is not 
available in the problem. We have seen that the Polyakov loop expectation value as such is 
not a physically observable quantity and the only order parameter available is the correlator 



of Polyakov loops (3.5) which is nonlocal. The situation reminds a little bit the Berezinsky - 
Kosterlitz - Thouless phase transition in 1-dimensional statistical systems ||5^ where also no 
local order parameter is present, but some nonlocal correlator changes its behavior at large 
distances. 

Another, a much more close analogy is the percolation phase transition which is well known 
in condensed matter physics (see e.g. |5l[]) and displays itself in the three-dimensional lattice 
Ising fiux model Let us discuss it in some more details. 

The partition function of the model is defined as 



Z= Vexp (3.19) 




^One of the way to see it is to use the string model for the hadron spectrum. A string state with large mass is 
highly degenerate. The number of states with a given mass depends on the number of ways p{N) the large integer 
N ~ M^/a [a is the string tension) can be decomposed in the sums of the form N = rii (see e.g. ^^). p{N) grows 
exponentially with N . That does not mean, of course, that in real QCD with large the spectrum would be also 
degenerate. Numerical calculations in QCD2 with adjoint matter fields show that there is no trace of degeneracy and 
the spectrum displays a stochastic behavior j4^. And that means in particular that there is little hope to describe 
quantitatively the QCD spectrum in the limit Nc — s- cx) in the string model framework. But a qualitative feature 
that the density of states grows exponentially as the mass increases is common for the large Nc QCD and for the 
string model. 
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where discrete variables 9i = 0,1 are defined on the hnks of the 3D cubic lattice. When 6i = 1, 
the link carries a flux with energy a and, when 9i = 0, there is no flux. An additional constraint 
that the sum of the fluxes 9i over the links adjacent to any node is even so that the flux lines 
cannot terminate. 

When the temperature is small, configurations with nonzero di are suppressed, only few 
links carry fluxes and they are grouped in a set of rarefied disconnected clusters. When the 
temperature is increased, the density of clusters is also increased. The phase transition occurs 
when a connected supercluster of flux lines extending on the whole lattice is formed. 

The model ( |3.19 ) (as well as the three-dimensional Z2 gauge lattice model which bears a 



considerable resemblance to the lattice SU (2) gauge theory) is dual to the 3D Ising model. Its 
partition function just coincides up to a coefficient with the partition function of the latter: 

Z = A{P) ex-p <p'^{sijSi+ij + SijSi j+i)y (3.20) 

{si} [ ij ) 

with tanh/3 = e~'^^'^ and where dual spin variables Si dwell on the nodes. Dual spins are 
ordered at large temperature T. At the critical temperature and below it, ordering disappears. 
But, as dual spins are not physically observed variables, < Si > cannot serve as a physical 
order parameter. The order parameter for the original theory is the percolating flux density 
which is easy to visualize (it is the probability that a given link belongs to the percolating 
supercluster), but no nice analytic expression in terms of original variables 9i can be written 
for this quantity. Deconfining phase transition for the model ( p.l4|) resembles in many respects 
the phase transition in this example — though the order parameter which is the string tension 
is easy to visualize, it cannot be expressed via physical dynamic variables ^(r, n). f\ 

However, an analogy is not yet the proof. What is quite definite, of course, is that the 
potential between heavy quarks grows linearly at large distances at T = and is screened in 
plasma phase. The most simple and natural assumption is, indeed, that a change of regime 
occurs at a finite critical temperature T^. On the other hand, we want to emphasize that 
it is impossible now to conclude from pure theoretical premises that Tc 7^ 0. A queer but 
admissible theoretical possibility would be that the interquark potential is screened at any 
finite temperature however small it is. 

Let us discuss in more details the lattice results which display a phase transition at finite 



Tc- According to |31| 



^ \ 0.69 ±0.02, 5C/(2) 
Vo^ \ 0.56 ±0.03, SU{3) 

where a is the string tension. The scaling (i.e. the fact that the ratio Td \fo does not depend 
on the lattice spacing when the latter is small enough) has been checked which indicates that 



^We will see in Chapter 5 that also the chiral restoration phase transition for the theory with light quarks involves 
a "percolation" in Euclidean space driven by instantons. Only in the latter case percolation occurs not at high, but 
at low temperatures. 
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the phase transition temperature ( |3.21D is for real and is not just a lattice artifact. The result 



( 3.21 ) has been confirmed later in ||5^ (more recent results exceed previous ones by 10 -15 % 
but it is not an issue here). What is rather surprising, however, is that the critical temperature 
turned out to be very small. Tc should be compared with the mass M of the lowest 0^"*" glueball 
state which has been measured to be 3.5^/a for SU{3) (see e.g. |5^). At Tc ~ 0.16M, the 
system presents a very rarefied gas of glueballs (its density being proportional to exp{— M/T}). 
They almost do not interact with each other and it is difficult to understand why the system 
undergoes a phase transition at this point. |^ 

The measurements in [31| were performed with a standard lattice hamiltonian ( p. 14 ). Note, 
however, that one can equally well consider the lattice theory with the hamiltonian having the 
same form as ( |3.14| ) but involving not the unitary but the orthogonal matrices V°''^^{r,n) G 
50(3). Both lattice theories should reproduce one and the same continuous Yang- Mills theory 
in the limit when the inverse lattice spacing is much greater than all physical parameters (As 
far as I understand, there is no unique opinion on this issue in the lattice community. If, 
however, lattice hamiltonia involving unitary and orthogonal matrices would indeed lead to 
different field theories in the continuum limit, it would mean that the Yang- Mills field theory 
is just not defined until a particular procedure of ultraviolet regularization is specified. This 
assertion seems to me too radical, and I hesitate to adopt it.). 

In the strong coupling limit T ^ ^uitr. the two lattice theories are completely different. 
The theory with orthogonal matrices has the same symmetry properties as the continuum 
theory , and there is no Z2-symmetry whatsoever. The effective dual hamiltonian depending 
on the gauge transformation matrices Or'^"' G SO (3) also has no such symmetry and there is 
nothing to be broken. 

It was observed it [^] that the phase transition in the lattice system with fundamental 
SU{2) matrices is associated with the spontaneous symmetry breaking in the dual hamiltonian 
i?^'^'^ (fijf""'^') (We repeat that such a breaking is not a physical symmetry breaking because it 
does not lead to the appearance of domain walls detectable in experiment.). In our opinion, 
however, the additional Z2-symmetry which the hamiltonian ( 3.14| ) enjoys is a nuisance rather 
than an advantage. It is a specifically lattice feature which is not there in the continuum theory. 
We strongly suggest to people who can do it to perform a numerical study of the deconfinement 
phase transition for the theory involving orthogonal rather than unitary matrices. In that case, 
no spontaneous Zjv breaking can occur. Probably, for finite lattice spacing, one would observe 
a kind of crossover rather than the phase transition. The crossover is expected to become more 
and more sharp as the lattice spacing (measured in physical units) would become smaller and 
smaller. 

It would be interesting also to try to observe the "walls" (i.e. the planar Euclidean in- 



Actually, the same lattice data which indicate the presence of the phase transition at the low temperature (3.21) 
indicate also that glueballs do interact at T = Tc. The lattice value of the energy density at T = Tc is 5 times larger 
than the energy density of free glueball gas . And it is still completely unclear why the interaction of glueballs is 
so strong while their tree level density is so low. 
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stantons) for the orthogonal lattice theory. They should "interpolate" between Or'^-' = 1 and 
Or'^-' = 1 along a topologically non-trivial path. Like any other topological effect, these in- 
stantons should become visible only for a small enough lattice spacing (much smaller than the 
characteristic instanton size), and to detect them is definitely not an easy task. But using 
the orthogonal matrices is the only way to separate from lattice artifacts. The only available 
numerical study |5^] was done for the theory with unitary matrices and too close to the strong 
coupling regime (the lattice had just two links in temporal direction |^ ) where these artifacts 
are decisive. Thereby, it is not conclusive. 

A comprehensive lattice study of the deconfinement phase transition with orthogonal ma- 



trices has been performed so far. Recently, a paper appeared where the question was 
studied in the mixed SU (2) theory involving both the standard Wilson action with unitary 
matrices and the term with orthogonal matrices: 

5mixed = -^^TV^C/(n) -^Y.^TAU{a) (3.22) 

where U{0) is the product of four unitary matrices on a given plaquette and the adjoint 
trace Tr^^ is {2Trp)'^ — 1. Surprisingly, it was observed that the phase transition which was 
of the second order at zero Pa becomes the first order at some admixture of the adjoint 
term in the action. It was argued in |58] that the authors of |57] observed actually the so 



called bulk first order transition [|59[ which is a pure lattice artifact and that a more accurate 
study indicates (not yet quite definitely) the presence of both first order bulk phase transition 
and the deconfinement second order phase transition in temperature. In principle, these two 
phenomena can be sorted out on lattices of larger size. 

In recent the question was studied in the SU{2) theory with pure adjoint action. 
Unfortunately, the numerical accuracy was not so high for the critical temperature in the 
continuum limit (as the string tension is zero in pure adjoint theory due to screening of adjoint 
sources, Tc should be measured in the units of the lowest glueball mass) to be determined. 

Thus the experimental situation is still not quite clear now. In our opinion, we would not 
be quite sure that the deconfinement phase transition occurs at a finite temperature, and that 
this temperature is indeed six times smaller than the glueball mass until it will be confirmed 
in lattice experiments with pure orthogonal action. 



4 Lukewarm pion gas. 



"'^"'^We do not discuss here the measurements of interface tension between confined and deconfincd phases in SU{3) 
at Tc which is a perfectly well defined physical quantity. These measurements were performed by many people and 
at larger lattices. 
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4.1 Chiral symmetry and its breaking. 

If the theory involves besides gluons also quarks with finite mass, the static interquark potential 
VQQ{r) does not grow at large distances anymore even at T = 0. Dynamic quarks screen the 
potential of static sources. One can visualize this screening thinking of the color gluon tube 
stretched between two static fundamental sources being torn apart in the middle with the 
formation of an extra quark-antiquark pair. Thus in QCD with quarks, the Wilson loop 
average has the perimeter rather than the area law |^. The correlator of two Polyakov loops 
( |3.5| ) tends to a constant at large distances universally at low and at high temperature, and 
this correlator cannot play the role of the order parameter of phase transition. 

Still, the phase transition can occur and does occur in some versions of the theory. It is 
associated, however, not with change in behavior of the correlator (^), but with restoration 
of chiral symmetry which is spontaneously broken at zero temperature. We shall discuss the 
dynamics of this phase transition in details in the next chapter. Here we will concentrate on 
the properties of the low temperature phase. 

Let us first remind the well known facts on the dynamics of QCD at zero temperature. 
Consider YM theory with SU{3) color group and involving Nf massless Dirac fermions in the 
fundamental representation of the group. The fermion part of the lagrangian is 

/ 

where I?^ = — igA^j^t"" is the covariant derivative. The lagrangian (|4.1| ) is invariant under 
chiral transformations of fermion fields: 

<1L,R Al,r qL,R (4.2) 

where qL,R = ^(1 i 'j^)q are flavor vectors with Nf components and ^l,_r are two different 
U{Nf) matrices. Thus the symmetry of the classical lagrangian is UL{Nf) ® Uii{Nf). Not all 
Nother currents corresponding to this symmetry are conserved in the full quantum theory. It 
is well known that the divergence of the singlet axial current j^ = Yl,f Qfliil^lf is nonzero due 
to anomaly: 

d,jl ~ g^e^'^^^G^Gl^ (4.3) 

Thus the symmetry of quantum theory is SUiiNf) SUR{Nf) ® [/^(l). It is the experimental 
fact that (for Nf = 2,3, at least) this symmetry is broken spontaneously down to Uv{Nf). 
The order parameter of this breaking is the chiral quark condensate SU{Nf) matrix 

S//' = < qLfQRf >0 (4.4) 

-"^^That does not mean that there is no confinement ~ - as earlier, only the colorless states are present in the physical 
spectrum. But the behavior of the Wilson loop is not a good signature of confinement anymore. 
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By a proper chiral transformation ( [4.2[ ) it can be brought into diagonal form |^ 

In the following, the term "quark condensate" will be applied to the scalar positive quantity 
S. 

Note that the phenomenon of spontaneous chiral symmetry breaking is specific for theories 
with several light quark flavors. In the theory with Nf = 1 , the non-anomalous part of the 
symmetry of the lagrangian is just Uv{^)- It stays intact after adding the mass term and after 
taking into account the formation of the condensate < qq >■ The condensate is still formed, 
but it does not correspond to spontaneous breaking of any symmetry. 

For Nf = 2,3, spontaneous breaking occurs and this is a non-trivial feature of QCD. 
There is no way to derive rigorously from general premises that the chiral symmetry should be 



spontaneously broken. Indeed, recent lattice measurements |62, S3| indicate that the symmetry 
is probably not broken at all in the theory with four massless quarks. (Actually, such a theory 
is easier to analyze on the lattice than the theories with 2 or 3 light flavors: four flavors arise 
quite naturally in Kogut- Susskind approach due to well-known doubling of massless fermion 
lattice species.) Quark condensate was measured to be very small and compatible with zero. 
The calculations in instanton model (see a detailed discussion of this model in the next chapter) 
also indicate that at Nf = 4 (may be, at Nf = 5) the quarks condensate disappears p^ . 

As was just mentioned, by now we are not able to derive theoretically that the symmetry 
should be broken at Nf = 3 and should not be broken at Nf = 4, 5, but the trend — the 
larger is the number of light flavors, the more difficult it is to break the symmetry — is easy 
to understand. 

Condensate and spectral density 

The argument is based on the famous Banks and Casher relation connecting quark 
condensate to the mean spectral density of Euclidean Dirac operator p{\) at A ~ 0. Let us 
explain how it is derived. Consider the Euclidean fermion Green's function < q{x)q{y) > in a 
particular gauge field background. Introduce a finite Euclidean volume V to regularize theory 
in the infrared. Then the spectrum of massless Dirac operator is discrete and enjoys the chiral 
symmetry: for any eigenfunction tpnix) satisfying the equation pipn = ^ni'n , the function 
ipn = T^V'n is also an eigenfunction with the eigenvalue = — A„. 

The idea is to use the spectral decomposition of the fermion Green's function with a small 
but nonzero quark mass 

/ >r- ^n(a;)^i(y) , . 

< q{x)q{y) >= 2^ , ^ (4.6) 



^•^Not any complex matrix Nf x Nf can be brought into the form ( |4.5| ) by a unitary transformation. But the 
condensate matrix is subject to a constraint that the vector SU{Nf) flavor symmetry is not spontaneously broken 



as it follows from the Vafa- Witten theorem All eigenvalues of an admissible condensate matrix (4^) are equal 
and, after diagonalization, it is proportional to a unit matrix, indeed. 
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Set X = y and integrate over d'^x. We have 

V < qq >= -2m ^ — ^ (4.7) 

where the chh'al symmetry of the spectrum has been used and the contribution of the zero 
modes = has been neglected (it is justified when the volume V is large enough p^]). 
Perform the averaging over gauge fields and take first the limit V ^ oo and then the limit 
m ^ 0. The sum can be traded for the integral: 

E = - <qq>= 2m r -^^^dA = 7r/)(0) (4.8) 
Jo \ + m"^ 

The rightmost-hand-side of Eg. ([4.^) is only the non-perturbative m- independent part of the 
condensate . There is also a perturbative ultraviolet-divergent piece oc mA'^^^. which is pro- 
portional to the quark mass, is related to large eigenvalues A and is of no concern for us 
here. 

Thus the non-perturbative part of the quark condensate which is the order parameter of 
the symmetry breaking is related to small eigenvalues of Euclidean Dirac operator. There 
should be a lot of them — a characteristic spacing between levels is 6X ~ 1/{T,V) which is 
much less than the characteristic spacing 5X ~ 1/L for free fermions. 

The average spectral density p{0) appears after averaging of the microscopic spectral density 

PaW = ^E'^(^-^-[^]) (4-9) 

n 

over all gauge field configurations. The weight in the averaging involves a fermion determinant 
factor 



[Det(-i p + m)ff 



A„>0 



Nf 

(4.10) 



where v is the number of the zero modes of the massless Dirac operator (by Atiyah-Singer 
theorem, it coincides with the topological charge of the gluon field). For small m and small A, 
this factor is small. Thus the configurations with small eigenvalues are effectively suppressed, 
and the larger Nf is, the more prominent is suppression. It is quite conceivable that at some 
critical Nf the suppression of small eigenvalues becomes so strong that the averaged spectral 
density at A = disappears. 

The suppression of small eigenvalues displays itself in the dip of spectral density at small 
eigenvalues which sets in starting from Nf = 3. An exact result for the spectral density in 
QCD with Nf > 2 massless quarks in the region A <^ ^thadr (/^hadr is the characteristic hadron 



scale ~ 0.5 GeV) reads [67] 



p(A) = - + — ^^^|A|+o(A) (4.11) 
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The result dOl] ) was derived assuming chiral symmetry breaking occurred. We see, however, 



that the spectral density decreases as A approaches zero, and the larger is Nf — the larger 
is the effect (for Nf = 2 the slope is zero, and the value Nf = 1 is beyond the region of 
applicability of this formula). The scenario when p{0) hits zero at Nf > 5 is quite probable. 

Quenched QCD 

It is instructive to make a digression and to discuss here the quenched case Nf = 0. In that 
case we have a pure glue theory without dynamical quarks. Thus there is no chiral symmetry 
in the lagrangian and nothing can be broken, but one is allowed to couple the quark field as an 
external probe and consider the "quark condensate" (|4.7|), (^). Again, spectral density p{0) 
is given by averaging the microscopic spectral density ([4.9| ) over gluon fields where the measure 
does not involve now the determinant factor (|4.10| ) and small eigenvalues are not suppressed 
whatsoever. We think that in the quenched case the fermion condensate is actually infinite 
and that the spectral density behaves as 

p{X) ~ 1/A (4.12) 



at small A |68|. Three arguments in favor of this conjecture can be suggested: 

• This happens in the Schwinger model. Q For standard Schwinger model with dynamical 
fermions, the condensate and p{0) have nonzero finite value. They are proportional to the 
dimensional coupling constant g with a known coefficient. But for the quenched Schwinger 
model (pure 2-dimensional photo dynamics) the condensate of massless external probes 
diverges, and the spectral density behaves as 1/A. This exact result can be understood 
looking at characteristic field configurations in the path integral for the quenched vs. 
unquenched case ||68[| . 

• /f the chiral condensate would be finite in the quenched QCD, the Gell-Mann - Oakes 



Renner relation ( [4.15 ) would imply the presence of massless pions, i.e. a power falloff 



of the correlator of quark pseudoscalar currents at large distances. In our opinion, such 
a power falloff is unnatural in the quenched case: one could expect the appearance of 
the Goldstone states when a real physical symmetry of the theory acting on the dynamic 
fields in the lagrangian is broken, but not when the symmetry acts on the fields presenting 
external probes. Indeed, in the Schwinger model, the correlator < ■ip'ip{x) TpipiO) >quenched 
grows exponentially at large distances and no trace of would be "pions" is seen. 

An additional observation is that path integral has also contributions from topologically 
non-trivial instanton-like fields. The latter involve fermion zero modes A„ = which, 
as is seen from ( [4.7] ), provide the contribution ~ l/m in the quark Green's function 
and the quark condensate. For A'^^ ^ 0, this large factor was multiplied by the fermion 



14 



The latter is similar to QCD in many respects. It involves confinement, chiral symmetry, the quantum Ua{^) 
anomaly, and presents an excellent playground. 
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determinant factor in the measure ( 4.10D which involved the smah factor my^i and there 



was no divergence. In the quenched case, the divergence ~ 1/m in the condensate which 
suggests the divergence ( 4.12| ) in the spectral density survives. 



For sure, this argument alone is not sufficient. It proves that the condensate is infinite in 
the chiral limit m ^ when the Euclidean volume V is kept fixed but cannot formally 
exclude the scenario when the condensate tends to a finite value in the thermodynamic 
limit y — > oo with fixed m (as was mentioned above, the contribution of zero modes 
becomes irrelevant in the thermodynamic limit). Such a behavior would still be very 
strange. A normal physical situation (which is realized in QCD with A'^^ > 2) is when 
the condensate, the order parameter of a spontaneously broken symmetry, disappears 
rather than grows to infinity in the chiral limit at finite volume. ( In QCD with Nf = 1 
when the condensate does not play the role of order parameter, its value is the same in 
the thermodynamic and in the chiral limit). 

All that together makes us to believe that the conjecture ( ^.12| ) is true |^. That means in 
particular that the "quenched chiral theory" |^9|, which was based on the assumption inferred 
from old lattice measurements that the quenched condensate is finite and which implies that 
the quenched pion is massless, may be is relevant for finite lattices but is not realized in the 
continuum limit. 

Recently, lattice people started to see that the condensate diverges, indeed, in the chiral 
limit. The lattice data suggest the behavior p{X) ^ l/\fX for the spectral density |7C], but 



more study is definitely needed. Also I have to mention that the instanton model calculations 
give probably even divergent, but much more mild behavior of ^quenched^^^ near zero ||7l[ [the 



prediction (4.11) is qualitatively confirmed]. May be, it is one of few places where the instanton 
model gives a qualitatively wrong behavior ? Probably, to clarify the situation, it would make 
sense to perform instanton-like simulations in the quenched Schwinger model where the exact 
result is known. 

Chirally symmetric phase. 

If the chiral symmetry is not broken spontaneously as is suggested for Nf > 5 by numerical 
calculations, the dynamics of the theory has nothing in similar with QCD. First of all, the 
absence of the condensate implies the absence of Goldstone states. Pseudoscalars are massive. 
A QCD inequality was derived which tells that the baryon mass is always larger that the 
pseudoscalar mass. Thus the baryons are also massive. Then the only massless particles 
which are required in the physical spectrum to saturate the global axial anomaly are massless 
colored quarks. The theory is not confining. Actually, the same conclusion follows also without 



""^^ Still, it does not have a rank of theorem. 

^^The inequality tob > 3/4 m^r for Nf > 6 has the rank of exact theorem |Q. The result niB > TOtt for any Nf 
follows if invoking a very natural additional assumption. A more stringent constraint ms > Nc/2 follows from 
the analysis in the framework of the constituent quark model and/or string model [[73|. 
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invoking QCD inequalities. One cannot saturate the anomaly with massless baryons because 



the latter fail to fulfill the 't Hooft matching conditions for Nj > 3 [74| and one needs massless 
quarks anyway. 

The absence of confinement in a nonabelian theory with asymptotic freedom is somewhat 
unusual, but there is neither paradox nor catastrophe here. Actually, this phenomenon has 
been known for a long time (see a recent detailed discussion of this question in [^^). The 
point is that asymptotic freedom does not necessarily implies infrared slavery. Consider the 
2-loop beta function of QCD with A'^^ massless quarks. It has the form ||7^ 

= (ll - ^^/) - 4^ - y^/) (4.13) 

Suppose Nf = 16. Then the first coefficient is negative (and that corresponds to asymptotic 
freedom) and very small. The second coefficient is positive and is not particularly small. When 
we evolve the coupling in the infrared, the evolution would be very soon affected by the second 
term in (l4.13| ). We see that the equation (3{as) = has the solution 

as = — (4.14) 
151 ^ ' 

which is rather small. At this value, the third term in the beta function is still much smaller 
than the second one [the second and the first are of the same order due to accidental smallness 



of the first coefficient in (4.13)]. This is the infrared fixed point. The coupling is frozen at 
the small value ( 4.14 ), perturbation theory works and the spectrum involves colored massless 
quark and gluon asymptotic states. No trace of confinement. Recently, a similar phenomenon 
(the so called "conformal window") was discussed in the context of supersymmetric QCD with 
sufficiently large number of light fiavors |77]. 



Unfortunately, no fixed point is seen in the perturbative beta function ( 4.13| ) for Nj = 4, 5. 



The second coefficient in the beta function is still negative here and the equation /3(qs) = has 
no solution. (One can get a solution including higher orders in the perturbative expansion of 
(3{as) in a particular renormalization scheme, but, besides it is scheme-dependent, the stable 
point would correspond to large values of where perturbative expansion makes no sense) 
Thus the absence of confinement at Nf = 5 is still surprising. More studies (lattice, instanton, 
etc.) are required. 

Effective chiral theory. 

For Nf = 3, the spontaneous breaking occurs and leads to appearance of the octet of 
pseudoscalar Goldstone states in the spectrum. Of course, the quarks are not exactly massless 
in real QCD, and the mass term is not invariant with respect to the symmetry ([4.2|). As a 
result, in real world we have the octet of light (but not massless) pseudo-goldstone pseudoscalar 
states (tTjK,!]). But the small mass of pseudogoldstones and the large splitting between the 
massive states of opposite parity {p/Ai, etc.) indicate beyond reasonable doubts that the exact 
chiral symmetry ([4.2[) would be broken spontaneously in the massless case. As the masses of 
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the strange and, especially, of u- and d- quarks are small |78] , the mass term in the lagrangian 



can be treated as perturbation. E.g. the pion mass satisfies the relation |^ 

F^M^ = (m„ + mrf)S (4.15) 

(Ftt = 93 Mev is the pion decay constant) and turns to zero in the chiral limit m^^rf 0. 
It is noteworthy that the symmetry breaking pattern 

SULiNf) SUniNf) ^ SUviNf) (4.16) 

depends crucially on the assumption that the gauge group involves at least 3 colors. For SU{2) 
color group where quarks and antiquarks belong to the same representation (the fundamental 
representation of the SU (2) group is pseudoreal: 2 = 2), the symmetry group of the lagrangian 
( |4.1| ) is much higher. It is U{2Nf) and involves also mixing between quarks and antiquarks. 
C/yl(l) - part of this symmetry is anomalous and the formation of chiral condensate breaks 
spontaneously the remaining SU{2Nf) down to a simplectic group [^ ]: 

SU{2Nf) Sp{2Nf) (4.17) 

As a result, 2N'j — Nj — 1 Goldstone bosons living on the coset space appear. For Nf = 2, we 
have not 3 as usual, but 5 "pions". Two extra "pions" are diquarks. This fact is important 
to understand for people who would wish to study numerically on lattices spontaneous chiral 
symmetry breaking with SU{2) gauge group. 

The presence of light pseudogoldstones in the spectrum is of paramount importance for the 
physics of low temperature phase. When we heat the system a little bit, light pseudoscalar 
states ( pions in the first place) are excited, while p-meson, nucleon and other massive degrees 
of freedom are still frozen. We can study analytically the properties of the system at low 
temperatures when the medium presents a rarefied weakly interacting gas of pions with low 
energies. Their properties are described by the effective nonlinear chiral lagrangian 

£ = ^F^TT{df,Ud^U^ + SReTr{7W[/t} + ... (4.18) 

where U = exp{2it"' cp"" / Ft^} is the SU{Nf) matrix {(j)"' are the pseudogoldsone fields), Ai is 
the quark mass matrix and the dots stand for higher derivative terms and the terms of higher 
order in quark masses. When the characteristic energy and the quark masses are small, the 
effects due to these terms are suppressed and a perturbation theory (the chiral perturbation 
theory |^, ^]) can be developed. 

The particular nonlinear form of the lagrangian (|4.1^ ) is dictated by its symmetry proper- 
ties: it is invariant under multiplication of U by an arbitrary unitary matrix on the right or on 
the left in the massless case . This symmetry exactly corresponds to the symmetry (|4.2|) of the 
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We assumed here that the quark mass matrix is diagonal. Otherwise, the ground state of the hamiltonian would 



correspond to a nondiagonal quark condensate matrix (4.4). See the detailed discussion of the "disoriented quark 
condensate" in the next chapter. 
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original lagrangian. One can also realize the symmetry (4.2) linearly including an additional 
a field. Actually, it is equivalent for many purposes to the nonlinear a - model approach, but 
with some particular assumptions on the coefficients of higher-derivative terms in ( [4.18 ) (they 
appear if integrating a-field out). Expanding ( [4.18 ) up to the quartic terms in 0" and using 
( [4.15[ ), we obtain for Nj = 2 

C 



lio.rf - f ( 



1 



6F2 



b\2 



+ 



24F2 



I a la \ 2 



+ 



(4.19) 



When is zero, the vertex of the four-pion interaction involves the second power of momen- 
tum. That means, in particular, that the cross section of scattering of massless pions occurs 
at P-wave and behaves as ~ E'^/F^ at smah energies. 



4.2 Thermodynamics 

After these preliminary remarks about the physics of zero- temperature QCD, we are ready 
to discuss what happens if we switch on the temperature. In this chapter, we assume the 
temperature to be small so that only the lightest degrees of freedom, the pions, are excited. 
Their density is small and their interaction is weak. The basic thermodynamic characteristic 
of a finite T system is its free energy. We have in the lowest order 



-TlnZ 



-Tin 



n 



-fSnEp 



Kn=0 



(4.20) 



where Ep = + M'^ and (Nj — l) is the number of degrees of freedom of the pseudogoldstone 
field with the mass (we derive free of charge the result for any A'^^, but, physically, one 
should put Nf = 2 or, if Nf = 3, take into account the different masses of tt, K and r/ mesons). 
The sum J2n e~^"'^p is the free energy of a single boson field oscillator with a given momentum. 
Trading the sum for the integral, we obtain for the volume density of the free energy 



1)T 



(2^ 



In 



(4.21) 



A very important corollary of this simple formula is that the order parameter of the spon- 
taneously broken chiral symmetry, the quark condensate, decreases with temperature. The 
temperature- dependent condensate S(T) is defined as a logarithmic derivative of free energy 
with respect to quark mass: 



S(T) 



1 1 dF{T) 
WfV dm 



(4.22) 
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(we endowed all quark flavors with a common small mass m to be set to zero after difl'erenti- 



ation). The expansion of ( 4.21 ) in pion mass in the region ^ T reads 



V 



+ + 0{MlT) 



90 



24 



(4.23) 



Substituting it in (4.22) and taking into account the Gell-Mann 
( 4.15 ), we obtain 



s(r) = S(0) 



1 T2 



12iVj 



0{T'/Ft] 



Oakes - Renner relation 



(4.24) 



Actually, a similar phenomenon has been known for a long time in condensed matter physics. 
At finite temperature, spontaneous magnetization of ferromagnet starts to fluctuate and its 
average value decreases with temperature. Technically, these fluctuations are described by 
massless magnon degrees of freedom which are excited in the heat bath and contribute in the 
free energy (magnons are the Goldstone particles which appear due to spontaneous breaking 
of rotational symmetry in ferromagnet by the same token as pions appear due to sponta- 
neous breaking of chiral symmetry in QCD). Only the shift in spontaneous magnetization is 
proportional to T^/^ ( see e.g. [|l9|, p. 299) rather than to as in our case. 

Also higher terms of the expansion of S(T) in T-^/F^ were found. To this end, one has 
to calculate the two loop and three loop corrections to the free energy density which take 
into account pion interactions. This was done in |^], but we choose to postpone the 
presentation and the discussion of these results to the next chapter. 



4.3 Pion collective excitations. 

At T = the particles display themselves as asymptotic states. As was mentioned before, 
asymptotic states do not exist in heat bath. Their role is taken over by collective excitations. 
This notion can be easily visualized when recalling the familiar college physics problem of 
propagation of electromagnetic waves in the medium. The frequency and the wave vector of 
such a wave do not satisfy the vacuum relation u; = c|k| anymore. A refraction index may 
appear. The amplitude of the classical wave decreases with time due to dissipative effects. 
That means that the wave frequency acquires a negative imaginary part (a positive imaginary 
part would correspond to instability). 

Technically, collective excitations display themselves as poles in the retarded thermal 
Green's functions. The position of the pole a;poie(k) is called the dispersive law of the cor- 
responding collective excitation. 

^^In the system involving free electrons (plasmas) a mass gap can develop. That means that plasma supports a 
standing electromagnetic wave with k = and nonzero frequency (it is called the plasma frequency). Quark-gluon 
plasma has a similar property, and we are going to discuss it at length in Chapter 6. 
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Let us find the dispersive law for pion collective excitations in the heat bath. To this end, 
we have to calculate the retarded polarization operator including thermal effects and solve the 
dispersive equation 



= U)' 



,2 _ k2 _ 



nH(u;,k) = 



(4.25) 



The simplest graph is shown in Fig. ^ where the 4-pion interaction vertex can be inferred from 



(4.19). The corresponding calculation was first done in |83]. As a result, it was found that the 



pion mass acquires a temperature correction 



Ml{T) = m2(0) 1 + 



rp2 



(4.26) 



24F2 



+ ... 



It is not difficult to generalize this formula on the case when we have Nf light quarks with 
the same mass (it is not so, of course, in the real world). The dependence is qualitatively the 



We see that the pion mass slightly increases with temperature, but, if the quarks masses 
and -M7r(0) would be zero, pion mass would remain zero also at finite temperature. This is easy 
to understand: massless pions are Goldstone particles appearing due to spontaneous breaking 
of chiral symmetry. But if the temperature is small, chiral symmetry is still spontaneously 
broken and the massless Goldstone particle should still be there as it is. 



same, only the coefficient of the 



term is changed from 1/24 to 1/(12A'^^). 
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Figure 3: Pion polarization operator in one loop. 



It is instructive to see how this result appears in the direct calculation of the l~loop graph 
in Fig. ^. The quickest way is to use the Keldysh technique. Actually, we do not need here the 
full- -scale Keldysh technique. The situation is the same as in the model example discussed 
at the end of Chapter 2. II/j = IIii -|- , but 12 - component is zero here and it suffices to 
consider only the temperature-dependent part in the T-ordered pion propagator in the loop. 
For massless pions, the latter is —27rinB{\p\)5{p'^). The 4-point vertex involves two momenta. 
If one or two of these momenta is the loop momentum, the corresponding integrals 

r d^p Pf,6{p'^) r d^p p'^5{p'^) 

^ J (27r)3 e/5|pl-l' ^ ^ J (27r)3 c^^IpI - 1 

are just zero. The only contribution to the polarization operator comes from the term where 
both momenta sit on external lines. This contribution has the form 

n„(.,|k|) = ^ / ^ = (c^^ (4.27) 

Substituting it in the dispersive equation ([4.25| ), we find that uj = |k| is still a solution. Only the 
residue of the pole is renormalized. (Incidentally, this residue, in contrast to the pole position, 
does not have a universal meaning. Its value depends on the particular parameterization 
U = exp{27ri(/)'^t'^/F7r} used. In Weinberg parameterization [^] it would be different.) 

Another temperature effect is the renormalization of F-,^. The invariant way to define F.,^ 
at finite temperature is via the residue of the polarization operator of axial currents. At large 
distances, the latter behaves as 

< A^(x)A,(0) >r = « F^{T) exp{-M^(r)|x|} (4.28) 
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a) 



b) 



c) 



Figure 4: Axial correlator in one loop. 



The graphs which contribute to the axial correlator on the one-loop level are depicted in 
Fig. ^. The first one was already discussed for the problem of pion mass renormalization. 
Here we need to know, however, also the renormalization of the residue of the pion propagator 
which contributes in the renormalization of the residue of the axial correlator (T) . 

The latter depends also on two other graphs in Fig.^ [which do not affect the pion propaga- 
tor]. The corresponding vertex < 0|j4^|37r > can be found by "covariantizing" the derivatives 
in dirsD in a proper way [g]: d^U ^ df,U - i{Af, + Vf,)U + iC/(V^ - A,,) where = Ap^, 

= V^t"-, and A'^, are external sources coupled to axial and vector currents. The result 
coming from the sum of all three graphs is parameterization-independent and has the form 



Nf 



+ 



(4.29) 



Note that the temperature-dependent condensate given by Eg. (PI, the pion mass given by 
Eq.( [4.26|) , and Ft^{T) given by Eq.( [4.29| ) still satisfy the Gell-Mann - Oakes - Renner relation 
( |4T5D . 
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Figure 5: Pion polarization operator in two loops. 

The next in complexity is the two-loop graph depicted in Fig. ||. The essential novel 
feature is that two-loop polarization operator has not only real, but also an imaginary part. 
Correspondingly, the solution of the dispersive equation (|4.25| ) becomes complex, the imaginary 
part of ix'poie(k) describing the damping of pion waves. |^ This damping was first found in [^]. 

Let us derive it in a slightly different way than it was done in [^^, using Keldysh technique. 
The analytic expression corresponding to the graph in Fig.|5| is 

\Tik,pi,p2,P3)\^ [Duipi)DuiP2)DuiP3) - Duipi)DuiP2)Di2{P3m.30) 

where T{k,pi,p2,P3) is the 4-pion vertex and the negative sign of the second term is due to 
the fact that the corresponding vertex comes from the anti-T-ordered exponential [cf. ( |2.30| )]. 
|Tp summed over all flavor polarizations in the loop has the form 

where s = {k — = {k — p2)'^,u = (k — ps)^ and we suppressed the trivial flavor factor 



S""""' . Substitute now the expansion ( |2.37 ) for the matrix components Dn, D12 in Eq.( 4.30 ). 



Note that the terms involving products Dii{pi)Da{p2) etc. give zero after integration. 
Really, fix p3 and consider the integral 

~ j Ai DR{pi)DA{k -P3- Pi) (4.32) 

Recalling that Da = D*^ and that Dr which has the form ( p. 15 ) is analytic in the upper u) 



half-plane, it is easy to see that the integrand in ( 4.32 ) is analytic in the upper piQ half-plane. 
Closing the contour there, we get zero. 



^^In this review, we consistently define damping (^(k) = — Im Wpoic(k) as the coefRcient determining the exponential 
attenuation of the amplitude of the classical wave with time. People often discuss also the twice as large quantity 
7(k) which is the coefficient in the exponential time decay of energy or probabifity. 
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Having that in mind, one obtains after some transformations 

ImUnik) = -- f ^ ^ i2n)U^'\k - pr - P2 - P3) \T\' 

^ 24: J (27r)4 (27r)4 (27r)4 ^ ^ ^ i^i y2 \ \ 

[Dr{pi)Dr{p2)Dr{pz) + Dr{pi)Dp{p2)Dp{p3) + Dp{pi)Dr{p2)Dp{p3) 

+ Dp{pi)Dp{p2)Dr{p3)] (4.33) 



Now we have to take the imaginary part of ( 4.33| ) and subtract its zero-temperature value (the 



latter does not affect the dispersive relation and does not concern us here). Substituting the 
expression ( 2.3^ ) for Dp and noting that Im Dr{p) = —irdijP' — m?)e{po) [e(po) is the sign 



factor], we obtain 



(27r)4 (27r)4 {2tt) 
5{pI-MI) 5{pI-MI) 5{pI-MI 



coth — i- coth — ^ — eiijji)e (uj2 ) 
2T 2T ^' 



+ ^^^^^^^^^'^ - e{u}i)e{uJ3) + coth ^ coth ^ - e{u}2)e{uJ3) 



e(a;i)e(w2)e(w3) (4.34) 



where uJi = piQ. The integral acquires contributions both from positive and negative Wj. Note 
the essential difference with usual T = Cutkosky rules where only positive energies in the 
direct channel are allowed. Physically, negative energies correspond to the particles not in 
finite but in initial state — there are a lot of them in heat bath. In our case, the contribution 
from the kinematic region where all frequencies are positive is just zero. Indeed, we want 
eventually to substitute the polarization operator ( |4.34 ) in the dispersive equation ( 4.25| ) to 



find a sma// imaginary pole frequency shift: u;pQjg(k) = ujp^^{]s.) — ^C(k)• We are allowed then 
to calculate the integral assuming uj = a;o(k) = •\/k^~+Mj. In that case kinematics does 
not allow to have all frequencies positive: a massive particle cannot decay in 3 particles of the 
same mass. A massless particle can, in principle, do it, but all final particles would go parallel 
to each other and the corresponding phase space is zero. Nonzero contribution comes from the 
kinematic region when one of the frequencies is negative and two others — positive. Choose 
for definiteness wi < and W2,3 > 0. Physically, it corresponds to scattering of the ingoing 
pion with the energy uj on the pion from heat bath with the energy —uJi producing two other 
pions with energies uj2 and ojs in the final state. Taking into account three possibilities 
(either coi or 0^2 or ^3 is negative) and using the nice relation 

11 11 11-, sinh/t 

coth Ki coth K2 + coth Ki coth K3 — coth K2 coth K3 — 1 



sinh Ki sinh K2 sinh K3 ' 



^"You may ask: how can we talk of scattering now while we stated earlier that scattering matrix is an ill-defined 
notion at finite temperature. That is true, and ^-matrix does not exist as a transition amplitude between asymptotic 
states. But now we are dealing with a "scattering on the way" . We need not follow the fate of scattered pions up 
to i = 00. The wave loses energy and never reaches infinity exactly because of the damping effects we are currently 
studying. 
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where k = k,2 + — ki, we eventually obtain 

C(k) = = S^^Mm I rf„,rf„,<i.3 ^2.f sHk+P^-P2-n) \T\\ 

= Wf AEMm ' -=4^^'' i^. = v'SjTrf (4.35) 

(all energies are positive now). 

One readily checks that, if the gas is dilute which is true in the region T <^ M^r, this 
formula reduces to the relativistic expression for the colUsion rate 

^^""^ = \j (^exp(-^i/Ty°*^-i(l-vvi) (4.36) 

where v, Vi are the velocities of the two collision partners and where v"^^^ is the velocity of one 
of them in the rest frame of the other. (7*°* is the total cross section summed over the flavors 
of the initial pion from the heat bath. 
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Figure 6: Pion polarization operator in the first order in density. Vertical dash stands for the 
temperature insertion. 

This result may be obtained by the virial expansion over the pion density which corresponds 
to taking into account the temperature correction —2TrinB{Ei)5{pf — M^) in only one of the 
virtual pion propagators in the loop in Fig. ^. The sum of the graphs in Figs. |3|, |5| and many 
other graphs with only one temperature insertion in a virtual pion line can be drawn as in 
Fig. ^ where the grey blob stands for a zero- temperature amplitude of forward tttt scattering 
T^'^(s) summed over the isotopic species of one of the pions. We have 

n/?(u;,k) = - J nsiE) T-^s) + 0(n|) (4.37) 

Imaginary part of T'^'^{s) gives the total cross section and we arrive at the result ( [4. 361) (when 
the density nsiE) is small, it coincides with the Boltzmann distribution exp{—E/T}). 

Substituting |rp in the form ( 4.31| ) in the general expression ( [4. 351 ), an angular integral 



can be done and we obtain |87] 



^^^^ ^ sinh(o;/2r) vr 



8kuj (47rF^)4 

^ Jau ^ smh{Ei/2T) sinh(^2/2T) sinh[(a; + Ei - E2)/2T] ' 

h = {k+ - k^)[2{u + E,f - 9il<] - ^{kl - kl){u: + E,f + ^(fc^ - fc^ ) , 

h = iq+ -q-)iu;- E^f - \{ql - ql){u - E,f + ^iql - d) , (4.38) 

where 

fc+ = mm{k+pi, P2+P3), k- = max(|A;-pi|, \p2-P3\), 
q+ = mm{k+p2, pi+ps), Q- = max(|A;-p2|, Ipi-Psi), (4.39) 

/c = |k| and Pi = |pj|. This complicated expression can be greatly simplified in the theoretically 
interesting region T ^ u> ^ M^^. Assuming also Ei ^ T ^ co and expanding over to/Ei, uj/T, 
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we obtain 

1,2 foo rE\ f7'2_i_9p2 



^""^^^ 192tt3TF^Jo "^^Vo '^^'sinh(Si/2r) sinh(^2/2r) smh[(£;i - £;2)/2T]^^"^°^ 
The integral here involves a logarithmic divergence at E2 ~ or at i?3 ~ Ei — E2 ~ 0. 



This divergence was absent in the original expression (4.3S) and is actually an artifact of our 
assumption £'2,3 ^ ^- It is effectively cut off at -E2,3 ~ 1^ ~ A;. Finally, we obtain with a 
logarithmic accuracy 

Achiral/^N ^ JLj^ 1 i / ^ "'^ = ^ ^ ln±- (4 41) 

^ ^' 247r3Fy kJo smhHx/2) ISnFy k ^ ' 

This is not, however, the end of the story. At very large characteristic times (much larger than 
the relaxation time), the system is in hydrodynamic regime and the propagation of a classical 
pion wave is described by an effective wave equation with dissipative term 

vf — Att — CAtt + higher derivatives terms = (4.42) 

From this, one gets the behavior ( = Ck"^ at small k and there should be no logarithmic 
singularity. The matter is that, in the region k < ({k), one is not allowed to restrict oneself 
with the graph in Fig. |^, and higher order corrections are important. Physically, one cannot 
assume in this region that the imaginary parts of loop propagators in Fig. ^ are (5-functions, 
but should take into account the "smearing" of (5-functions due to the nonzero imaginary part 
CiPi) of the poles of the propagators. For very small k it is the damping itself rather than k 
which provides an infrared cutoff in the integral (4.40). The result ( 4.41D has a limited range 
of applicability 

rp5 



24Fi 



« T < ^hadr (4.43) 



where / {2AF^) is the approximate value [^] of the damping (^{k) averaged over the thermal 
distribution. At still lower k, the damping is estimated with logarithmic accuracy as 

2k'^T^ 2 2F 

Note that a similar damping-induced cutoff plays an important role in some kinetic problems 
for quark-gluon plasma. We will discuss this issue in Chapter 6. 

In the general case when the temperature is not assumed to be much larger than cj and 
the latter is not assumed to be much larger than the pion mass (and, as far as the real QCD 
system is concerned, the region T ^ M-,^ is well above the phase transition temperature, and 
the soft pion approach is not justified there), the integral in ( |4.38 ) can be done numerically. 



The results of Ref. [^] for the absorption length A = '^^^^ ^ i.e. the distance over which 
the energy density of the classical pion wave is attenuated by the factor 1/e , are displayed in 
Fig. 0. 



"^The term "mean free path" used in is not quite exact. The latter is usually applied not to the problem of 
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Figure 7: Absorption length as a function of pion momentum. The dashed hues describe the 
behavior in chiral hmit, the sohd hues correspond to the experimental value of and the dotted 
line gives the qualitative behavior in the theory. 

Also the real part of the pion dispersion relation is modified on the two-loop level. To study 
this, one has to calculate the real part of the polarization operator ( [4.33| ) from the graph in 
Fig. |5| and also take into account the graph in Fig. ^ where the box stands for the vertices 



involving four momenta in the derivative expansion of the effective chiral lagrangian (4.18). 



absorption of tfie wave with a definite value of /c, but ratfier to the problem of relaxation of a wave packet with a 
broad momentum distribution. The latter depends not on the total cross section as the absorption length, but on 
the so called transport cross section ct*'' ~ J dcr(l — cos^). An additional factor 1 — cos 6* suppresses the contribution 
of scattering at small angles. In the problem under consideration, the differential cross section (4.31) has a broad 
angular distribution, and the mean free path is of the same order as the absorption length. In other systems, in 
particular in quark - gluon plasma where the Rutherford differential cross section is strongly peaked around 6 = 0, 
the difference of these two quantities is very essential. We will discuss it in details is Chapter 6. 
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Figure 8: Contribution of four-derivative vertices in the pion polarization operator. 

An explicit calculation along these lines has not yet been done so far. A. Schenk |^8| cal- 
culated the thermal shift of Re uj(k) in a "mixed approach" : he used the virial expansion 
formula ( 4.37| ) and substituted there the forward scattering amplitude T'^'^(s) (and a similar 3 



3 forward scattering amplitude which enters the term of the second order in pion density) 
as it follows from the chiral perturbation theory. Numerical plots for Re w(k) at different tem- 
peratures assuming nonzero pion mass were drawn. Unfortunately, the results were presented 
in a way which makes an analytic analysis difficult. 

In the old paper by Itoyama and Mueller and in the recent [^], the same problem was 
studied in the linear sigma-model framework. A beautiful qualitative effect was discovered: 
the dispersion law of slow massless pions is modified to 

Rea;(k) = c(r)|k|, 

where M„ is the mass of cr-meson which was assumed to be small compared to 27rF^. A 
"refraction index" appears. It is not clear for us now whether this result is reproduced also in 
the nonlinear sigma-model [If it is, the correction should be of order T^/F^ with a possible 



dependence on the coefficients of higher-derivative terms in ( 4.18 )]. But I do not see a reason 
why it should not. 



In pO| also the two-loop contribution to the residue F^^iT) in (4.28) was calculated and 



the validity of the Gell-Mann - Oakes - Renner relation ( [4.15| ) to the order ~ was checked. 



4.4 Nucleons. 

Not only pions but also other hadron states change their properties when temperature is 
switched on. Nucleons, vector mesons etc. are not easily excited when the temperature is 
small, but still it makes sense to study the problem of propagation of a massive state in pion 
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gas and find out how the presence of the thermal heat bath affects its dispersive properties. 



Such a study was first carried out for nucleons in |91]. 

We restrict ourselves with the case when the nucleon is at rest and will be interested with 
the temperature dependence of the pole position of the nucleon Green's function $7^(0) = 
MJj — i(Jf. The real part of the pole will be called the temperature-dependent nucleon mass 
and the imaginary part Cn{T) — the nucleon damping rate. 

What we need to know is the nucleon polarization operator at finite temperature. Virial 
expansion in the pion density is very handy here. In the first order in density, the polarization 
operator is given by the graph of the same form as in Fig. and, quite analogously to Eq. 



( 4.37 ), we obtain 



E = \/M'^ + p-^. vrA^ scattering is very well studied in experiment. We may use the available 



phenomenological information |95, Q on the vrA^ scattering amplitude and need not come to 
grips with evaluating it in the chiral perturbation theory framework. The only assumption in 
(4.46) is that the heat bath include only pions and that their density is small enough for the 



virial expansion to be justified. In reality, these conditions are fulfilled up to T 130 MeV. 
The limitations are the same as for the mesonic sector. 

Consider first the damping rate. According to the optical theorem, the imaginary part of 
the forward scattering amplitude is related to the total cross section. 



Imr^iv(^) = 2M%^E^-M^a''^{E) (4.47) 

We have 

" 4^ Jm '^^ - I K+p(^) + ""^''P^^^ + ^^-p{E)\ (4.48) 

To evaluate the integral, we insert the experimental cross section, which in the energy region 
of interest is dominated by the A-resonance. In the narrow resonance limit, this gives 

C5 = (4.49) 

where Tr = 115 MeV is the width of the resonance and E^ ~ 330MeV is the corresponding 
pion lab energy. At the peak of the resonance, the cross section saturates the unitarity limit, 
^TT+p _ ^^i^^q^y.^ where qr is the cm. momentum at resonance. Using the isospin ratios 



^^Note that at finite temperature, alternative definitions of mass can be considered which in general are not 
equivalent to the one adopted here (see e.g. jo^, ^). In particular, the correlation length at large space-like 
distances is not given by the inverse of Af^. The two quantities differ even for free fcrmions where the pole position 
is temperature independent while the correlation length is given by (M^ + tt'^T^)^^/'^. [according to (2.8), the 



Euclidean fermion frequencies are quantized to an odd multiple of ttT] . It is important to understand when theoretical 
predictions are compared with the lattice data [p3 . 
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afP : afP : < p = 1 : | : |, this implies ct^ ~ 380 mb. Eq. (|^) is plotted in Fig. | 
(dashed line) together with a numerical evaluation of Eq. ( 4.48[ ), based on the Karlsruhe- 
Helsinki data (full line). Damping rapidly grows with temperature. 



53 



MeV 



200 



100 



150 



50 




T (MeV) 



50 



100 



150 



Figure 9: Damping rate 7^ = 2^^ of nucleon collective excitations as a function of temperature. 
The dashed line indicates the contribution due to the A- resonance. 

It is instructive to compare this result with a situation where the physics is weh understood: 
atoms exposed to a heat bath of photons. The virial expansion in the photon density can be 



cohision of an atom with a photon from the heat bath can be written. This cross section grows 
as at smah frequencies. As a resuh, the atom damping grows as T^. This estimate is correct 
while the temperature is stih much less than the characteristic distance between atomic levels 
~ TTiga^. If it is not the case, excitation and ionization processes set it. In other words, when 
the temperature of the heat bath is high enough, the atom does not live in the ground state 
most of the time, but is constantly excited by thermal photons and gets eventually ionized. 

Also in our case, when the temperature is high enough, the nucleon does not stay in the 
ground state, but is constantly converted into A and higher excited baryon states, absorbing 
and emitting pions. Our calculation of the nucleon damping makes contact with reality up 
to the temperatures of order T ~ 150 MeV. At this temperature, the ratio of probabilities 
for finding the nucleon in excited state and in the ground state, respectively, is given by 
4exp[— (Ma — MAr)/T] ~ 0.6. Also, at temperatures of order 150 MeV and higher, the medium 
contains an increasing number of massive excitations and is not adequately represented by a 
dilute gas of pions. 

Consider now the shift of the real part of the pole compared to its zero temperature 



carried out and a formula similar to Eq.(4.48) but involving the total cross section a^A for the 
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value. We use again the virial expansion formula ( [4.46 ). The real part of the forward irN 



scattering amplitude involves the Adler zero: in the chiral limit, it behaves as oc E'^ when 
the energy of an ingoing pion is small. As a result, the thermal shift of the nucleon mass is 
proportional to at small temperatures. 

Note that the term oc is absent here in contrast to what we had for the quark condensate 
( 4.24 ) and F^^ ( [4.29 ). It particular, early attempts to use directly the loffe formula for the 



nucleon mass [^8[ derived from QCD sum rules 

Mn ^ 1.2(47r2s)^/^ (4.50) 



and substitute there the temperature dependent condensate ( [4.24 ) are not justified. ^ Speak 



ing of the pion mass, its relative shift (MJ — M^)/M^ is of order T'^/F^, but, as M^- is small, 
the absolute value of the shift is also small, and there is no shift at all in the chiral limit. 

Like the damping rate, the shift of the nucleon mass can be determined from the formula 
(4.46) without invoking chiral theory, but using the available experimental information on the 



ttN forward scattering amplitude [[95| . The result is shown in Fig. IC. For comparison, we also 
plot the chiral perturbation theory result for the temperature dependence of the pion mass 

(EH- 



^•^There are two reasons why, contrary to original hopes QCD sum rules method cannot be straightforwardly 
generalized to the finite temperature case. First, sum rules are contaminated by new parameters: thermal expectation 
values of Lorentz - non- invariant operators. The second and even more important reason is that the standard 
"resonance + continuum" model for the spectral density of a current correlator does not hold at finite temperature 
anymore. At T 7^ 0, the spectral density acquires contributions due to scattering of ingoing current on the pions from 
the heat bath |10C|. One can show that this additional contribution to the spectral density plus the temperature 
renormalization ( 4.53| ) of the nucleon residue exactly match the thermal shift of the quark condensate ( 4.24| ) in the 

See [103 for a related 



theoretical side of the corresponding sum rule [ [LOl| (the same is true for vector mesons [ 102 
discussion. 
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Figure 10: Shift in the effective nucleoli and pion masses caused by temperature. 



The main feature here is that the shift in the real part of the nucleon pole position is 
small compared to the shift in the imaginary part in the region of applicability of the whole 



approach. This finding of |91] was confirmed in recent |104]. 



Pion exhibits itself as a pole in the axial current correlator ( 4.28 ). Likewise, a nucleon 
exhibits itself as a pole in a correlator of currents with nucleon quantum numbers. If we 
restrict ourselves with the currents not involving derivatives of the quark fields, two choices 
are possible [p^], 



mix) 
mix) 



^abc 



{ [iO^cd 



c 



iuD^Cdi] ul} 



(4.51) 



where C is the charge conjugation matrix, 7!^ 
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Figure 11: Renormalization of the nucleoli residue. Dashed hues stand for pions. 

Finite temperature brings about the renormahzation of the residues of these currents in 
the nucleon states A by the same token as it brings about the renormahzation of the residue 
of the axial current F-,^. (while the current 7?i(x) plays a preferred role in the context of QCD 
sum rules [^] , the temperature effects on the residue of these two currents are the same) . The 
corresponding graphs, which exactly correspond to the graphs in Fig. ^ in the pion case, are 



shown in Fig. 11 . Using the Goldberger-Treiman relation and the current algebra relation 



for the soft pion vertex < 0\rj\TT7rN >, 

< Olr^lA^vr'^Tr'' > = --— ^ < 0\r]\N >, (4.52) 

and calculating the graphs in Fig. ^ in the same way as we did it earlier for pions, we arrive 
at the low temperature theorem [^l| 

A- - - + ..] (4.53) 

Like Ft^, the residue drops with temperature. 

One can also derive similar formulae for other hadrons. For example, -F^ = Fk{1 — 
rV32F2 + ...), etc. 

4.5 Vector mesons. Experiment. 

Let us ask now whether the theoretical results presented in the two last sections can be 
directly confronted with experimental data in heavy ion collisions. Unfortunately, the answer 
is negative. 

We hasten to comment that the pole positions of the collective excitations with pion and 
nucleon quantum numbers are quite physical quantities and can be measured in a gedanken 
experiment. 



57 



Suppose we study the spectrum of invariant masses of two 7 emitted from hot hadronic 
matter. The spectrum has a sharp peak associated with the vr'^ 27 decay. For free pions, the 
width of the peak is very small ~ 8 eV. But thermal pions have a finite width F^ = 2(^J 
where is the thermal average of the expression ( [4.38| ) for CttC^)- Also the maximum of the 
resonance is shifted towards larger mass values according to Eq.( ^.26| ). Likewise, for nucleons. 
Proton is believed to decay eventually, and most of the decay modes (like e~^7r^) involve hadrons 
which are stuck within the hadronic medium and do not go out undisturbed, but there is also 
a mode e'^j. For the argument sake, let us assume that e"*" and 7 do not interact with the 
heat bath, but only with detector. Measuring the distribution in invariant masses of 6^7 
and focusing on the "nucleon resonance", one could, in principle, extract information on the 
position of the nucleon pole. 

Of course, these gedanken experiments are quite unrealistic. Not even speaking of the "nu- 
cleon resonance experiment" which is a pure science fiction, also the studying of 27 spectrum 
in heavy ion collisions would provide little information about pion properties in heat bath. 
The matter is that tt^ lifetime is at least six orders of magnitude larger than the lifetime of 
the fireball produced in the collision of heavy nuclei . Pions would decay on flight in vacuum 
and their mass and width would be the same as in any other experiment. 

As far as experiment is concerned, the situation is much better for vector mesons, especially 
for /3-meson. Its vacuum lifetime is pretty small, almost all p-s produced in the collision decay 
inside the fireball, and measuring the spectra of e~^e~ or fJ^^fJ-^ pairs can provide an information 
about the properties of p in hot hadron medium. Things are not so good with uj and (j) mesons. 



Lifetime of hadronic fireball is estimated to be 25 fm/c or less [105|. It is of the same order 
as the CO lifetime and almost 2 times less than (p lifetime. A considerable fraction of cos and 
most of 0-s would decay outside the fireball. 

Unfortunately, Life and Nature usually ingeniously resist our attempts to eliminate Trouble 
and put things in Order. A manifestation of this general Law in this particular case is that the 
problem to describe p-meson in hot matter involves much more theoretical uncertainties than 
in the pion or nucleon cases. The most straightforward way to study the position of p pole 
is to perform virial expansion in the pion density and to write the formula like (|4.37| ), ([4.46| ) 
involving the forward up scattering amplitude T^^p. The Trouble here is that the experimental 
information on T'jYp is not available. We have to calculate it in some theoretical model, and 
that brings about uncertainties. 

Something can be said, however. First of all, the amplitude of pion forward scattering on 
any hadron h involves the Adler zero due to Goldstone nature of pion. In the chiral limit, this 
amplitude summed over all pion species behaves as 

T^+h + T^oh + T^-h{E) ~ -^(A + iB^) (4.54) 

/^hadr ^ /^hadr / 

where ^hadr is some characteristic hadron scale, and we assumed that the hadron h has a 
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nonzero isospin. |^ It applies also to p. And that means that, like in the nucleon case, the 
real part of the p mass shift does not have a contribution ~ in the chiral limit. 

The same conclusion can be reached considering the correlators of vector and axial currents 
in pion heat bath. Some complication there is that the correlators themselves are shifted in 
the order ~ T^, but an accurate analysis displays that this shift occurs not due to a shift of 
p and Ai poles, but due to admixture of the graph with Ai pole in the vector correlator and 
the admixture of the graph with p pole in the axial correlator [106]. 

However, p-meson mass is shifted in the order ~ T^. It cannot be excluded right now that 
its temperature dependence is more profound than for the nucleon mass (Fig. |l^ shows that 
the latter is almost constant in the region T 150 MeV where the whole approach based on the 
virial expansion is justified), and it can, in principle, be essentially modified at temperatures 

T^lOO 150 MeV. Whether it is modified, indeed, and, if it is, does it eventually grow with 

temperature |107| or decreases | 108 , 109 1 — is the question under discussion now. 

Experimental data on dilepton mass spectra in heavy ion collisions ItLlC 



111] favor more the 



possibility that the effective mass of p decreases with temperature. Experimentalists observe 
an excess of dilepton pairs in the small invariant mass region in nucleus-nucleus collisions 
compared to the spectrum in proton-nucleus collisions. It seems that this excess cannot be 
explained in a conservative framework, with the processes of direct dilepton production in the 
hadronic fireball like vrvr — > pl^l~ or up ttI^I~ [112, 113]. The assumption that Mp drops 
with temperature can explain it |l0|, llTl. 



5 Chiral symmetry restoration. 

5.1 General considerations. Order of phase transition and crit- 
ical behavior. 

In QCD with 2,3, and possibly 4 massless quark flavors (see the discussion in Sect. 4.1), 
chiral symmetry is spontaneously broken at zero temperature. In the preceding chapter, we 
studied the dynamics of the system at comparatively low temperatures when chiral symmetry 
is still broken, and the spectrum of the system involves massless Goldstone states. But a spon- 
taneously broken symmetry is usually restored under sufficient heating, like a spontaneously 
broken rotational symmetry in ferromagnet is restored above the Curie point. 

There are few exceptions of this general rule. First of all, it does not apply to supersym- 
metry. Supersymmetry is always broken spontaneously at nonzero temperature irrespectively 
of whether it is broken or not at T = 0. The most immediate way to see it is to notice that dif- 

^*If the isospin of the target is zero, the imaginary part of the amplitude involves an additional suppression. The 
situation is quite parallel to the problem of soft photon scattering. If the target is charged (Compton scattering), 
the total cross section is constant in the limit uj ^ 0, and if the target is neutral (like the atom), the cross section 
behaves as ^ oj'^. 
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ferent boundary conditions in Euclidean time in the Matsubara formalism, periodic for bosons 
and antiperiodic for fermions, are not invariant under supersymmetry transformations which 
mix bosonic and fermionic fields. This is also seen in the real time hamiltonian formalism. The 
standard density matrix exp{—PH} is transformed under supersymmetry transformation into 
exp{—(3{H — Qe — eQ)} where the parameters of SUSY transformation e, e have the meaning 
of chemical potentials corresponding to conserved supercharges Q, Q [114]. Supersymmetry is 
spontaneously broken by the same reason why the Lorentz symmetry is spontaneously broken 
at finite T. In the latter case, the density matrix is non-trivially transformed under the Lorentz 
boost and acquires the form exp{— — vP)/\/l — v^} where P are conserved momenta and 
the velocity v plays the role of corresponding chemical potentials. 

Supersymmetry is a special case, but there are also some systems where spontaneous break- 
ing of a usual (bosonic) global symmetry persists until arbitrary high temperatures. Moreover, 
heating the system can even induce the spontaneous breaking which was absent at zero tem- 
perature [^]. Such systems are, however, rather exotic, and there are some special dynamic 
reasons for this peculiar phenomenon. There are no such reasons in QCD, it behaves as 
most of the physical systems do, and that means that a critical temperature Tc exists above 
which the chiral symmetry is restored and the fermion condensate < qq >t is zero. This is 
the temperature of phase transition and < qq >t is the order parameter associated with the 
transition. 

Note that the statement of the existence of the phase transition point does not apply to 
the theory with only one massless flavor. As was already noted in Sect. 4.1, quark condensate 
does not signal there a spontaneous symmetry breaking [C/^(l) symmetry is anyway broken 
explicitly by anomaly], is not an order parameter, and need not vanish at high temperature. 
So, it does not. At high temperatures when the effective coupling is small, it can be evaluated 
semi-classically in the instanton approach |5|, 115], and one can show that it falls down as a 
power of temperature and never reaches zero. 

Let us return to the case of several massless flavors when spontaneously broken chiral 
symmetry is restored at some temperature. The first question to be asked is what is the order 
of this phase transition, i.e. whether the order parameter < qq >t itself, or only its derivative 
d < qq >T /dT are discontinuous at the phase transition point. 

This question can be studied theoretically. There are rather suggestive arguments which 
indicate that the phase transition is of the second order for 2 massless flavors. When Nf > 3, 
the phase transition is probably of the first order. (If it is there, of course. We mentioned 
earlier that, for Nf > 4, 5, chiral symmetry is probably not broken at all, and there is nothing 
to be restored.) The arguments are the following [116]: 

The starting point is the observation that, in theories involving scalar fields, phase tran- 
sition of the first order often occurs when the potential involves a cubic in fields term. One 
can recall in the first place a cubic Van-der- Vaals curve P{p, T) which describe the first order 
water vapor phase transition. The simplest field theory example is the theory of real scalar 
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field with the potential 

V{ck) = \{ct>^-v^f-^i4>^ (5.1) 

Assume for simplicity <C Xv. At T = 0, the potential has one global minimum at (/> « 
V + 3;u/8A and a local minimum at (/> « —v + 3/i/8A. At nonzero temperature, the term 
~ \T'^4P' is added to the effective potential. At high temperature T ^ v, the only minimum 
occurs at (/) = 0. One can be easily convinced that, as the temperature increases, the left 
local minimum moves on the right and reaches the point i;^ = (after which it does not move 
anymore) at some temperature T* when the global minimum at positive (j) still exists. When 
the temperature is further increased, the right minimum moves up and becomes degenerate 
with the minimum at (p = at some temperature Tc- The former disappears altogether at a 
still larger larger temperature T** . In a certain temperature range, two minima of the potential 
coexist, one being a metastable state with respect to the other. This is exactly the physical 
situation of the first order phase transition. 

Let us go back to QCD. A direct application of this reasoning is not possible because the 
QCD lagrangian does not involve scalar fields. The effective chiral lagrangian (4.18) is also of 



no immediate use because higher-derivative terms which stand for dots cannot be neglected 
in the region close to critical temperature. Suppose, however, that in the region T ~ Tc some 
other effective lagrangian in Ginzburg-Landau spirit can be written which depends on the 
composite colorless fields 

^//' = QRfQLf (5.2) 

A general form of the effective potential which is invariant under SUiiNf) SUR{Nf) is 

V[<^] ~ 5iTr{^>$t} + c,2(Tr{$^>t})2 

+gsTr{m^m''} + 54(det ^> + det ^>1') + . . . (5.3) 

(the coefficients may be smooth functions of T). Now look at the determinant term. For 
Nf = 2, it is quadratic in fields while, for Nf = 3, it is cubic in fields and the effective 
potential acquires the structure similar to Eq.(^.lD which is characteristic for the systems with 



first order phase transition. A more refined analysis |116] shows that the first order phase 
transition is allowed also for Nf > 4, but not for Nf = 2 where the phase transition is of the 
second order. 



Most of the existing lattice data [117, 118, 119| indicate that, at Nf = 3, the phase 



transition is of the first order, indeed, while, at A'^^ = 2, it is of the second order. The latter is 
not, however, a generally accepted statement yet, and the discussion on this issue is still alive 
(we will discuss lattice data in some more details a bit later). 

Critical behavior 

A second order phase transition is characterized by critical behavior, i. e. a power-like 
scaling behavior of the specific heat, the expectation value of the order parameter, and the 
correlators. Different critical exponents can be defined. 
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Let us list them briefly here (for more details, see e.g. Chapter XXIV in |120[ |). Suppose 
we have a system with an order parameter < > which turns to zero at the second order 
phase transition point. Introduce also a weak external field h which is canonically conjugate to 
(j). For standard ferromagnet, < (p > \s magnetization and h is the magnetic field. For QCD, 
< (j) > is the quark condensate and h is the small common quark mass. 

The exponent a is related to the scaling behavior of the specific heat cp{T) slightly above 
or slightly below a critical temperature: 

cp{T) = ^±|r-r,|-° (5.4) 

The coefficients are called critical amplitudes and the index + or - refers to whether we 
are dealing with the system in the disordered phase T > Tc or in the ordered phase T < Tc. 

The scaling behavior of cp at T = Tc but at nonzero external field h is determined by the 
critical exponent e: 

cp{h) = A^h-' (5.5) 

When we approach from below, the order parameter tends to zero as a power 

<cP>T= Ap{T,-Tf (5.6) 

At r = Tc and h ^ 0, the order parameter also exhibits a power-like behavior 

<^>h= Ash'/^ (5.7) 

When T ^ Tc, the shift in < > provided by the external field h is proportional to h. The 
proportionality coefficient is called the generalized susceptibility x which scales as 

X = A^\T-nr (5.8) 

At T = Tc and h = 0, the mass gap in the physical spectrum is absent and the correlation 
length is infinite. The correlator of the order parameters falls down at large distances as a 
power 

< (/>(x)0(o) >T=n,h=o = ^cl^r^''"'^^^ (5-9) 

where d is the spatial dimension. At T ^ Tc and/or h ^ 0, the mass gap appears. It scales as 

Mgap = Af.h", T = Tc, h^O (5.10) 

and 

Mgap = AtlT-Tcl", T^Tc, h = (5.11) 
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Eight critical exponents a, 7, 5, e, /x, z/, are related by six scaling relations (so that only two 
of them are independent parameters): 



a + 2/3 + 7 = 2, 135 = P + -f, 
e(/3 + 7) = a, KP + 1) = 

z.(2-C) = 7, ud = 2-a (5.12) 

A very important property of the critical exponents is their universality. That means that 
the values of a, (5 etc. do not depend on the details of microscopic hamiltonian, but only on 
gross symmetry features of the theory. It follows basically from the fact that critical behavior 
is determined by the dynamics of the theory at distances which are much larger than the 
microscopic scale (much larger than 1//Xhadr in our case). At the critical point, the effective 
theory describing large-distance dynamics is a conformal theory not involving any dimensional 
parameter. At the vicinity of critical point, it is a perturbed conformal theory corresponding 
to adding to the effective lagrangian an energy operator or an order parameter operator with 
small coefficients. The effective theory thus obtained involves now a mass scale , but as long 
as it is small compared to a characteristic energy scale in the microscopic hamiltonian (/^hadr); 
this effective theory is still universal. 

And that means that, in order to determine the critical exponents in QCD with 2 light 
quarks at the vicinity of the second order phase transition point, we need not to deal explicitly 
with QCD which is a complicated theory, but are allowed to consider some other theory 
belonging to the same universality class, i.e. with the same symmetry breaking pattern. The 
exponents will be same. The universality class of QCD with Nj = 2 corresponds to the 
symmetry breaking pattern SUl{2) SUr{2) — > SU{2) which is the same as 0(4) — > 0(3). 
A rather simple example of the theory in this universality class is the so called Heisenberg 
0(4) magnet (Imagine a 3-dimensional cubic lattice involving 4-dimensional spin variable 
at each node. The hamiltonian is the interaction S^S'^ for all nearest neighbors). At low 
temperatures, spins are ordered < S"^ >x=o^ 0. When temperature is increased, the system 
undergoes a second order phase transition and goes over into the disordered phase. The 
values of critical exponents in this model are known [ [121| (they were obtained in e expansion 
technique). The universality argument tells us that their values in QCD are exactly the same 
m, |12|. Specifically, 

a = -.19 ±.06, /3 = .38 ±.01, 7 = 1.44 ±.04, 6 = 4.82 ± .05, 

e = -.10 ±.03, C = .03 ±.01, fi = .40 ±.02, = .73 ± .02 (5.13) 

A distinct feature of the model is the negative value of a so that the specific heat does not 
diverge at the critical temperature, but rather goes to zero with a cusp. 

The predictions ( 5.13| ) can be confronted with experimental lattice data. At the moment, an 



experimental situation is uncertain. The values of the exponents determined in |124] agree with 



the theory rather well. In particular, their result 6 ^ = .24 ± .03 coincides with the theoretical 
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prediction within the error bars. This finding was confirmed in recent [125] by another group. 
However, a preferred value for 6^^ given in a later work of the Bielefeld group done on larger 
lattices |126| ] is close to zero which actually indicates that the phase transition observed is of 
the first order and which sharply contradicts theoretical expectations. Also nobody has seen 
the predicted cuspy zero for the specific heat so far. One has to mention also the recent work 
|127| ] where critical exponents were measured in 2+1 - dimensional Gross-Neveu model. (Like 
QCD, it involves light fermions and enjoys a chiral symmetry. Like in QCD, chiral symmetry 
is broken spontaneously at zero temperature and is restored at high temperature) The values 
of the exponents measured in [127| do not conform with the predictions based on universality 
and the e-expansion (which is rather surprising as these predictions work remarkably well in 
the condensed matter systems where laboratory experimental information is in abundance), 
but are instead rather close to the values predicted in the Landau mean field model. 

Being not a lattice expert, I cannot judge, of course, who is right and who is wrong. 
Four-dimensional calculations are complicated. The most modern, specially dedicated for this 
purpose computers should work for hundreds of hours to produce statistically reliable results. 

However, there are also systematic errors in any experiment. Taking them under control 
depends on the skill of an experimentalist, and a layman has no say in that. The situation is 
much better when the experiment can be gauged, i.e. performed in the region where some solid 
theoretical results exist. If something known is reproduced in experiment with good accuracy, 
one can trust the experiment also when Unknown is measured. Fortunately, such a possibility 
(which has not been yet fully exploited) does exist for lattice numerical calculations. I mean 
here certain 2-dimensional models involving fermions where exact theoretical predictions can 
be made and where, on the other hand, lattice calculations are much more easy and require 
much less computer time. 

That refers in particular to Schwinger model (i.e. 2-dimensional QED) with several ( Nf > 
1) light fermion flavors. Much like as QCD, this model enjoys the chiral symmetry SUiiNf) 
SUR{Nf). In two dimensions, a continuous symmetry cannot be broken spontaneously due to 
the Mermin - Wagner - Coleman theorem |125]. So it does not and the quark condensate is 
zero in the massless theory. It can be shown, however, that at small temperatures and small 
fermion masses, the system displays a critical behavior like a system with second order phase 
transition at temperatures slightly above critical. One can say that the second order phase 
transition occurs at zero temperature (so that there is no ordered phase, indeed) ]129[ . The 
exact values of critical exponents can be determined. They are |130| ] 



a = -l, 7 = —, S = ^^, /.= ^^, u=l, C = 2-— (5.14) 

When Nf = 2, also, critical amplitudes entering ( |5.4| ) - (5.9) can be exactly calculated |131]. 
In particular, 

< V^iV'i >T=o = < ^2?/'2 >T=o = --388 . . . rri^/^g'^/^, 

Mgap = 2.008... m^/s^i/^ (5^5) 
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There is only one comparatively old lattice paper where Schwinger model with 2 flavors was 
studied [|132(| . The measured values of the critical exponents agree well with theoretical pre- 
dictions. However, the experimental values of the critical amplitudes exceed their theoretical 
values by 35% for the condensate and by 25% for the mass gap. It would be very interest- 
ing, indeed, to repeat this calculation on modern computers and with different modern lattice 
algorithms (Wilson vs. Kogut-Susskind fermions, perfect vs. improved vs. simple-minded ac- 
tions etc.). Confirmation of the theoretical predictions ( ^.14 , ^.15 ) with good accuracy would 
provide an efficient test of lattice methods. 

5.2 Insights from soft pion physics. Large A^^- 

We have already seen an indication for the chiral symmetry restoration in the previous chapter: 
according to Eq.( 4.24| ), the order parameter S(T) = — < qq >t decreases as the temperature 



is increased. That was a one-loop calculation which had only a limited range of validity 
T < 80 — 100 MeV where the temperature correction is still small. One can, however, 
improve the accuracy and find out the temperature dependence of the condensate on the 2- 



loop and 3-loop level in chiral perturbation theory. This calculation done in |83, qJ] allows 
one to monitor melting down of the quark condensate at larger temperatures and make a 
reasonable estimate for the temperature of phase transition where the condensate disappears. 
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Figure 12: Two loop contribution in the pion free energy. 





a) 



b) 



Figure 13: Some graphs contributing in the free energy on the three-loop level. 



Let us describe the results of |8J] where the temperature dependence of < qq >t has been 
determined on the 3-loop leveL To this end, one has to calculate the free energy density on the 
three loop level at small nonzero quark masses and differentiate the result over the quark mass. 
There is only one relevant two-loop graph depicted in Fig. |l^. By a simple power counting 
[We have T^M^ ~ T,m{T'^/F^) from the pion mass dependent term in the free energy and an 
additional factor ~ T'^/F^ from the vertex ( 4.19| ) involving two derivatives, i.e. the square of 
a characteristic momentum Pchar ~ T of pions in the thermal loop ], it gives the correction 
~ T'^/F^ in the condensate. There are many 3-loop graphs. We will depict two of them. The 



graph in Fig. 13a involves only the standard 4-pion vertex ( 4.19 ) and depends only on F.^ (it 



depends also on the ultraviolet cutoff, but disregard it for a moment). The graph in Fig. 13b 
contains something new. It looks as a two-loop graph, but contributes actually in the order 
~ because the square vertex in this graph stands for the terms in the effective lagrangian 



( 4.18 ) involving four derivatives. There are several chiral invariant structures of this kind. One 
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of them has the form 

= a(Tr{d^U^d^u}y (5.16) 

The dimensionless coefficient a (another name for it is Li [^^) is a new independent parameter 
of the lagrangian. In this order, there are 4 more such new dimensionless coefficients i>2,3,4,5 
which are relevant to our problem. All of them affect the temperature renormalization of 
the condensate in the order ~ T^. These new parameters are largely fixed, however, from 
experimental data on pion (and also K, r] ) interactions at intermediate range of energies 
which are also affected by Li [^]. |^ The final result for the condensate in massless theory in 
soft pion approximation (i.e. when the effects due to K and r] are disregarded together with 
effects coming from other resonances) has a rather simple form 



< 



>T=< 



>0 



rp4: 



rp2 



8F2 384Fi 



s-ln— + 

288F6 T 



(5.17) 



Here all effects from the higher-derivative term are described by the constant A. Experimental 



data on pseudogoldstone interactions give the value A ~ 500±100 Mev. The dependence ( 5.17 ) 
together with the curves where only the 1-loop correction oc or also the 2-loop correction 
(X are taken into account (please, do not put attention to the "technical" curve marked 



02 = 0) is drawn in Fig. 14 taken from Ref. 



Actually, Li appear not only from the higher- derivative vertices in the tree lagrangian, but also from countert- 
erms which are required to renormalize ultraviolet divergent pieces in the loop graphs of the kind depicted in FigJT^a. 
The effective lagrangian involves also many terms with six derivatives etc. On each new level of the expansion in 
momenta and in pion masses, many new independent constants appear, and their total number is infinite (which is 
not surprising — after all, the lagrangian ( 4.1§| ) is not renormalizable and involves an infinite number of countert- 
erms and subtraction constants.). The point is, however, that if we restrict ourselves by not too high energies (and 
we have to do it anyway: at higher energies, neglecting other than pions degrees of freedom is not justified), the 
contribution of these still higher derivatives terms is suppressed. Chiral perturbation theory is not the expansion 
over some dimensionless lagrangian coupling as usual, but rather an expansion over the parameters ~ Af^/(27ri^7r)^, 
~ I {2iT Ft^Y where p is a characteristic momentum scale of the process under consideration. 
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T 




Figure 14: Temperature dependence of the condensate in three loops. 

The expansion in the parameter T"^ /SF"^ makes sense when this parameter is smah, i.e. 
when T ^ 100— 150 Mev. Strictly speaking, one cannot extrapolate the dependence ( ^.17 ) for 



larger temperatures, especially having in mind that, at T > 150 Mev, the heat bath includes 
a considerable fraction of other than pion hadron states. But as we know anyhow that the 
phase transition with restoration of chiral symmetry should occur, the estimate of the phase 
transition temperature (i.e. the temperature when < qq >t hits zero, and we assume here 
that Nf = 2 so that the transition is of the second order) based on such an extrapolation is 
not altogether stupid. This estimate is 



Tc ^ 190 Mev (5.18) 

or, roughly, one inverse fermi. A more accurate treatment which takes into account nonzero 
rriu^d and also the presence of other mesons in the heat bath gives practically the same estimate 
for the phase transition temperature (assuming, of course, that the phase transition is there 
which, as we will see in the next section, is probably not the case): nonzero quark masses 
smoothen the temperature dependence of the condensate and push Tc up while excitation of 
other degrees of freedom sharpen the temperature dependence and pulls it down. These two 
effects practically cancel each other. 
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Large Nc- 

It is interesting to discuss what happens in the large N^. hmit |^^. As usual, we assume 
Nc ^ oo, — > such that the product g'^Nc stays constant. For simplicity, we also assume 
zero quark masses. 

In this limit, the spectrum of QCD presents a set of infinitely narrow meson resonances. 
In the first place, we have massless Goldstone states. An additional Goldstone particle is 
77'. When Nc is finite, ij' is massive because C^a(1) axial symmetry of the QCD lagrangian is 
broken explicitly by anomaly. However, in the limit Nc ^ 00, the effects due to anomaly are 
suppressed, Ua{^) becomes a good symmetry of the full quantum theory, and its spontaneous 



breaking by quark condensate brings about a new massless Goldstone state [133|. When Nc is 



large but not infinite, the mass of ry'-meson has the order M^' oc l/^/Nl |134|. Besides that, 
there are also massive meson excitations. The mass of the first such excitation = //hadr does 
not depend on Nc- The density of higher excited states grows exponentially with mass. 

Low energy properties of A'^^ Goldstone states are described by a modified chiral lagrangian 
which has exactly the same form as ( [4.18| ), only U is now a full U{Nf) rather than SU{Nf) 
matrix (there is also a determinant term |134| ], but it is suppressed in the large Nc limit). 
The temperature corrections to the quark condensate can be calculated by the same token 
as earlier, only the coefficients are modified a little bit. In particular, the first temperature 
correction to the condensate (4.24) involves now the factor N'j rather than NJ — 1. A glance on 



the large Nc analog of the expression ( 5.17 ) displays an apparent paradox. Indeed, scales 
as ^/Wc in the large Nc limit. And that seems to mean that the temperature corrections to 
the condensate become essential at T cx \/Ac/ihadr- Thereby, the condensate should turn to 
zero and chiral symmetry be restored at Tc ~ -v/AcA'hadr- 

This estimate for the phase transition temperature contradicts, however, physical intuition. 
The relevant physical scale in large Nc theory is not \/iVc/ihadr! but just /ihadr- In particular, in 
the infinite Nc limit, the temperature Tc ~ \/Ac/ihadr is much larger than the Hagedorn max- 
imal temperature Th ~ /Uhadr and just cannot be reached (see the discussion in the beginning 
of Sect. 3.3). 

If Nf > 3 (when the phase transition is of the first order), there is a way out. We may 
suppose that the condensate stays practically unchanged until Tc and then abruptly jumps 
to zero. But in the theory with Nf = 2, the phase transition is of the second order, the 
condensate should approach zero when we approach Tc from below, and we do not see yet how 
is it possible in the soft pion framework. 

The resolution of this paradox will be given in a moment, but first note that, after all, 
the fact that at, say, T = Tc/2 the condensate is practically not changed in the large Nc 
limit compared to its zero temperature value is quite natural and reasonable. Eventually, the 
condensate melts down due to nonlinear pion interactions. The latter can be estimated as 
being due to a resonance exchange, say, tttt p ^ tttt. Three-meson vertices scale as l/^/Wc 
in the large Nc limit (this suppression brings about the suppression oc 1/Nc for resonance 
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widths). Thus four-pion vertex scales as 1/Nc and, indeed, it is seen also from Eq.( [4.19 ). On 



the other hand, when the temperature is increased up to ~ /^hadn the density of thermal pions 
n ~ becomes so large that the characteristic distance between pions is of the same order 
as their size (the latter has the order ~ ^^adr ^"^^ does not scale with Nc). Pions overlap, and 
something is bound to happen at this point. 

And it does. We have seen that the correction cx in ( 5.17] ) involves, besides F.,^, also a 



new constant A which describes the contribution of the graphs of the kind drawn in Fig. |l3p 
involving higher-derivative terms in the chiral lagrangian. The corresponding dimensionless 
constants Lj, i = 1,. . . , 5 scale as Nc (Indeed, at Pchar ~ ^hadr, all terms in the pion 

lagrangian should be of the same order, and hence all coefficients in C'^'^\ C'^^\ etc. scale in the 
same way as F.,^ oc Nc- ). The contribution in S(T)/S(0) due to the graph in Fig. ^b is 
estimated to be ~ NcT^/F^, i.e., for large Nc, it dominates over the contribution coming from 
the graph in Fig. Oa depending only on F^. Staying on the 3-loop level, we would conclude 



1/3 1/2 

that the correction is of order 1 when T ^ Nc /^hadr which is better than T Nc ^hadn but 
is not satisfactory yet. 

However, one should, of course, include in the estimate also the contributions due to still 
higher order terms C^^\ C^^\ etc. Drawing the same graph as in Fig. but with a 4-pion 
vertex coming from (the corresponding coefficients scale as A^c//^hadr)' arrive at the 
estimate 



S(0)-S4-i-P(T) _ 1 ^ 1 iVc ^6 t 

^hadr 



Here the first factor is [S(0)]^^, ^hadr^^ comes from the expansion of the free energy in pion 
mass Fm oc M^T'^ and its subsequent differentiation over quark mass m, 1/F^ comes from 
the expansion of U up to the fourth order in cj)"', Nc/^i\^^^ is the estimate of a coefficient of 
six-derivative term, and is the sixth power of a characteristic momentum. 
Summing up the leading in Nc terms, we obtain 



s(r) ~ s(o) 



1 - ]v^^(-) 



(5.20) 



where F{x) is some function of a dimensionless parameter x = T'^ / fi'^^^.. We do not now how 
the function F{x) behaves at x ^ 1, but a quite natural assumption is that the Taylor series for 
this function has a finite radius of convergence so that F{x) becomes infinite at some xq. After 
all, the condensate is the derivative of free energy density over quark mass, and free energy 
definitely becomes infinite at Hagedorn temperature [see Eq.( 3.1^ )]. Let us assume further 



that F(x) has the simple pole at x = xq. Now if we keep temperature fixed and send Nc to 

^^If the lagrangian is written in terms of only 7V| — 1 Goldstone fields, a certain coefficient in called Ly 
scales as ||8^. But i) L-^ does not contribute in T,{T); ii) If it would, that would only help us to trace back the 
origin of large temperature correction; in) Anyway, we are thinking in terms of the effective lagrangian involving N'j 
Goldstone mesons. 
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infinity, the correction to the condensate vanishes. If, on the other hand, we keep Nc fixed, 
the quark condensate ( 5.20| ) would vanish and the phase transition would occur at some point 



Xcrit somewhat below xq, xq — Xcrit ~ l/-^c • The characteristic width of the region where S(T) 
is essentially changed is also estimated then as 

AT 1 , , 

If one assumes that F{x xq) ~ l/(x — xq)^, the width of the transition region would 
be somewhat larger (x l/Nc- If the singularity is weaker than l/(x — xq), the width of the 
transition region would be smaller (but is always suppressed at large Nc). 

Thus the paradox is resolved. We see that, at Nf = 2 and large Nc, the condensate is 
practically not changed, indeed, up to very vicinity of the phase transition point where it 
drops sharply and eventually hits zero. At Nf > 3, there is no "need" for the condensate 
to approach zero smoothly when the temperature approaches Tc- It can still be essentially 
changed in the transition region below Tc, approach some half-way value of order S(0)/2 and 
then abruptly jump to zero at the phase transition point. On the other hand, it can remain 
constant until the very phase transition point and then jump to zero "the whole way" . Finally, 
it can also approach zero when T — > Tc and the phase transition can still be the second order. 
(Ginzburg-Landau arguments presented above allow the first order phase transition at Nf > 3 
but do not dictate it to be the case.) Currently, no further conclusions on this point can be 
drawn. 

5.3 The real world. 

Up to now, we discussed only QCD with massless quarks. But the quarks have nonzero masses: 
rriu « 4 Mev, ~ 7 Mev, and mg ~ 150 Mev |7^]. The question arises whether the nonzero 
masses affect the conclusion on the existence or non-existence and the properties of the phase 
transition. 

Two different experimental (i.e. lattice) works where this question was studied are available 
now, and the results of these two studies drastically disagree with each other. The results of 
Columbia collaboration [117] display high sensitivity of the phase transition dynamics to the 



values of light quark masses. In Fig. 15, a phase diagram of QCD with different values of 



quark masses and niu = md as drawn in Ref.| 117 ] is plotted. 



71 




Figure 15: Phase diagram of QCD according to Ref. ||117|| . Mass values (measured in the units of 
the inverse lattice spacing for which the transition is and is not seen are denoted respectively 
by solid circles and squares. Dashed circle corresponds to the physical values of masses. 

Let us discuss different regions on this plot. When the quark masses are large, quarks 
effectively decouple and we have pure YM theory with SU{2>) gauge group where the phase 
transition is of the first order. When all the quark masses are zero, the phase transition is also 
of the first order. When masses are shifted from zero a little bit, we still have a first order 
phase transition because a finite discontinuity in energy and other thermodynamic quantities 
cannot disappear at once when external parameters (the quark masses) are smoothly changed. 

But when all the masses are nonzero and neither are too small nor too large, phase transition 
is absent. Notice the bold vertical line on the left. When = = and is not too 
small, we have effectively the theory with two massless quarks and the phase transition is of 



the second order. The experimental values of quark masses (the dashed circle in Fig.|l5|) lie 
close to this line of second order phase transitions but in the region where no phase transition 
occurs. It is the experimental fact as measured in Ref. |117|]. 



This statement conforms nicely with a semi-phenomenological theoretical argument of Ref. 



| 135 ] which displays that even z/the first order phase transition occurs in QCD, it is rather 



weak. The argument is based on a generalized Clausius-Clapeyron relation. In college physics, 
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it is the relation connecting the discontinuity in free energy at the first-order phase transition 
point with the sensitivity of the critical temperature to pressure. The Clausius-Clapeyron 
relation in QCD reads 

1 dT 

disc < qq >t„ = — 7r~~disc e (5.22) 

Tc dlTLq 

where disc e is the latent heat. The derivative can the estimated from theoretical and 
experimental information of how other essential properties of QCD depend on and from 
the calculation of T - dependence of condensate at low temperature in the framework of chiral 



perturbation theory (see Fig. 14 and the discussion thereof). The dependence on quark masses 
is not too weak: dTc/dnig ~ 0.9 — 1.0. From that, assuming that the discontinuity in quark 
condensate is as large as < qq >o, we arrive at an estimate 

disc e < 0.4 GeV/fm^ (5.23) 

This upper limit for disc e should be compared with the characteristic energy density of the 
hadron medium at T ~ Tc ~ 190 MeV. The latter cannot be calculated analytically and we 
have to rely on numerical estimates. The lattice study in [|136| (for the theory with 2 massless 
flavors) gives a rather large value e{T ~ 200 MeV) ~ 3 GeV/fm^. May be, this value is even 
too large, it is of the same order as the Stefan-Boltzmann energy (6.1, 6.3, |6.5|) of free quarks 



and gluons in QGP phase. This is surprising as, at such low temperatures, quarks and gluons 
are certainly not free and the interaction effects are important — see a detailed discussion 
in the next chapter. Anyway, one can safely conclude that the limit ( 5.23| ) for the latent 



heat is several times smaller than e(Tc). In reality, the discontinuity m < qq > at the phase 
transition point (if it is there) is, of course, much smaller than the value of the condensate 



at zero temperature. It is seen from the graph in Fig. 14. The condensate drops significantly 
still in the region where K and rj degrees of freedom are not yet effectively excited. These 
degrees of freedom (and also p and to degrees of freedom) become relevant only at T ^ 150 
MeV when the tendency of the condensate to drop down is already well established. Moreover, 
as was mentioned earlier in the paragraph after Eq.( 5.18l ), taking into account these degrees 



of freedom sharpens the condensate dependence rather than smoothens it. This is also seen 



from Eq.(4.24): The more is the number of massless quark flavors Nj, the sharper is the 
temperature dependence of the condensate. 

Thus latent heat of the first order phase transition in the theory with 3 massless quarks 
must be rather small [significantly smaller than the estimate ( |5.23| )] which means that the phase 
transition is likely to disappear under a relatively small perturbation due to nonzero rns- The 
numerical findings of Ref.[ [L17| are thereby rather appealing from theoretical viewpoint. 

However, as far as experiment is concerned, the question is far from being completely 
resolved by now. A recent lattice study |119[| done with Wilson rather than Kogut-Susskind 



fermions displayed quite a different picture shown in Fig. 16 
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Figure 16: Phase diagram of QCD according to Ref.[ll£]. First order phase transition is seen at 



sohd circles and is not seen at blank circles. The star marks out the physical values of quark masses. 
The bend of the second order phase transition hne separating the region with first order 



phase transition and the region with no phase transition reflects a recent finding [137] that 
the shape of this hne near the tricritical point nis = m*,mq = rriu = rrid = 0, where m* is the 
smahest value of the strange quark mass where we stih have a second order phase transition 
in a theory with 2 massless quarks is 

ml-nis ~ m^/^ (5.24) 



What is very surprising is the large value of m* ~ 450 MeV as measured in Ref.|119|. If they 



are right, the point corresponding to our physical world lies well inside the first order phase 
transition region in which case the first order phase transition should be of strong variety 



which contradicts Ref.[117] and also the theoretical arguments above. 



Again, we do not know what of these two measurements is more correct. Possibly, intrinsic 



lattice artifacts are more dangerous when an algorithm with Wilson fermions is used 1 137 , 135 ] , 
but I do not have my own opinion on this point. To resolve the controversy, it would make 
sense, first, as was already mentioned, to test lattice methods on a two-dimensional playground 
and, second, to perform a quantitative comparison of the temperature dependence of the quark 
condensate as measured on lattices with the theoretical curve in Figjl^. Such a comparison 
has not been done so far. 



^^This point is usually called tricritical because we can continue the phase diagram in the region of negative m^. 
Then the point (to*, 0) signifies the meeting of three regions: the region with first order phase transition, the region 
with no phase transition for positive m^, and the region with no phase transition for negative to^. Generally, the 
ratio vrLq/ms can be a complex number, and it is better to think of the point (m*, 0) as of the point where the surface 
of second order phase transitions in (|tos|, mq/mg) 3-space degenerates into the line. 
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To summarize, at the current level of understanding, the picture where the hadron gas 
goes over to quark-gluon plasma and other way round without any phase transition looks 
more probable. We have instead a sharp crossover in a narrow temperature range which is 
similar in properties to second-order phase transition (the "phase crossover" if you will). 



5.4 Instantons and percolation. 

In the discussion of the properties of the system at the vicinity of phase transition in this 
chapter, we relied so far on the fact that chiral symmetry is broken at zero temperature. It 
is an experimental fact in real QCD, but it is important to understand from pure theoretical 
premises why it is broken and what is the mechanism of its restoration at higher temperatures. 

A completely satisfactory answer to this question has not yet been obtained. The problem 
is that QCD at zero temperature is a theory with strong coupling and it is very difficult 
(may be impossible) to study the structure of QCD vacuum state analytically from the first 
principles. However, a rather appealing qualitative physical picture exists ||8|, |13S|] which is 
based on the model of instanton-antiinstanton liquid and on the analogy with the so called 
percolation phase transition in doped superconductors |[5l[| . 

The starting point is the Banks and Casher relation ([4.S|). As we have seen, nonzero quark 
condensate means nonzero spectral density of the Euclidean Dirac operator at A = and that 
implies the presence of rather small characteristic eigenvalues 

A ~ l/(i;y) < Afrec ~ 1/i (5.25) 

in the spectrum. The question is what is the physical reason for these small eigenvalues to 
appear. 

As far as we know, the first pioneer paper where a mechanism for generating small eigenval- 
ues was proposed is Ref. [14C] where small eigenvalues appeared as zero modes of monopole-like 



gauge field configurations. The disadvantage of this model is that the monopole configurations 
are static whereas it is natural to expect that characteristic gauge fields contributing to the 
Euclidean path integral at T = are more or less symmetric in all four directions with no 
particular axis being singled out. The model of instanton-antiinstanton liquid formulated in 



1 141 , 142| and developed later in | 143 | is much better in this respect. We refer the reader to 
a recent comprehensive review |13£] for details and just elucidate here main physical point of 
the reasoning. 

The basic assumption of the model is that a characteristic gauge field contributing in QCD 



path integral is a medium of instantons and instantons as shown in Fig. 17. It is not a "gas" of 



Callan, Dashen, and Gross |144| because the interaction between instantons and antiinstantons 



bringing about a short-range correlations between instanton positions and orientations cannot 



be neglected. A "liquid" is a more proper term |141] 
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Figure 17: Instanton-antiinstanton liquid. 

Each individual instanton and antiinstanton involves a fermion zero mode [ |145U . Assuming 
the constant density of quasi-particles cx ^^^dr ' total number of zero modes in the Euclidean 
volume V is N = Nj + Na ~ ^t^tadr- However, these are not exact zero modes. They are 
shifted from zero due to interaction between instantons and antiinstantons. In other words, 
individual instanton and antiinstanton zero modes ipj^\{x) are not eigenfunctions of the full 
Dirac operator = ^^{d^ — iA^f^ liquid^ ^ where 

Ni Na 
^iAiiquid(^) = ^^/(^_^^) + Y^A^ix-x^) (5.26) 
i=i j=i 

If choosing ip^j''\x — Xi) and ip^^ ~ ^j) ^ basis in the Hilbert space, the Dirac operator 
presents a matrix with off-diagonal elements 

~ J d^x {x - Xj) -pi^f^ {x - Xi) 

~ / - X,) Pi^fix - Xj) (5.27) 

Diagonalizing this matrix gives eigenvalues which are not zero anymore, but are spread 
over the characteristic range AA ~ ^^hadr- Assuming a uniform spread, the volume density of 
quasi-zero modes is p{&) ~ N/{V/S.\) ~ l^hadr- -^^^ t° Eq.( |4.8| ), a nonzero quark condensate 
appears. ^ 

^^The assumption of quasi- uniform spreading of eigenvalues is not so innocent. It probably holds only in the theory 
with several light dynamical quarks Nf < 4, 5, but not in the quenched theory {Nf = 0) where it is natural to expect 
a singular behavior of the spectral density near zero: p{X) ^ 1/A so that the "fermion condensate" (i.e. the vacuum 
expectation value < qq >o where quark fields are treated as external sources) is infinite In the theory with 

Nf > 4, we expect on the contrary a suppression of the spectral density at A ~ (see Sect. 4.1 for detailed discussion 

)■ 
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This picture is rather similar to what happens in a doped semiconductor with high enough 
doping. When a characteristic distance between individual atoms of the admixture is not 
large, the wave functions of outer electrons of these atoms overlap, and the electrons can jump 
from site to site. If the set of atoms of admixture with a noticeable overlap of wave functions 
forms a connected network in the space, the electrons can travel through this network at large 
distances and the specimen is a conductor. Note that it is not a standard metal mechanism of 
conductivity when the medium is a crystal, has the long- range order, and the electron wave 
functions are periodic Bloch waves. Here the distribution of the dope whose electrons are 
responsible for conductivity is stochastic and wave functions are complicated. The essential is 
that they are delocalized. 

Thus one can say that the vacuum of QCD is the "conductor" in a certain sense. For 
sure, there is no conductivity of anything in usual Minkowski space-time. Only the Euclidean 
vacuum functional has "conducting" properties. In principle, one can introduce formally the 
fifth time and write an analog of Kubo formula for conductivity in QCD, but the physical 
meaning of this "conductivity" is not clear. It is sufficient to say that, in a characteristic 
Euclidean gauge field background, the eigenfunctions of Dirac operator corresponding to small 
eigenvalues are delocalized. It would be interesting to check the latter statement explicitly for 
the lattice vacuum configurations or in the instanton model framework. Instanton study of 



this question is now in progress |14G]. 

What happens if we heat the system ? The effective coupling constant g'^{T) decreases, the 
action of individual instantons S = Stt'^ / g^{T) increases, and the density of quasi-particles oc 
exp{— S} decreases. Correspondingly, the characteristic value of off-diagonal matrix elements 
( |5.27D in Dirac operator decreases. Note also that the overlap integrals ( 5.27| ) decrease as 
the temperature increases even if the distance between instanton and antiinstanton is not 
changed. The matter is that the individual fermion zero modes of thermal instantons |147] ] 
fall down exponentially ~ exp{— vrTr} at large spatial distances, not just as a power ~ 
as is the case at T = 0. Both effects lead to suppression of overlap matrix elements at finite 
T. (Concrete numerical calculations |64, 139 1 show that the second effect is more important 
in the temperature region T ~ 200 MeV. Speaking of the instanton density, it falls down only 
by a factor 2 in this region. That means, in particular, that the value of the gluon condensate 
< (as/7r)G^ >T at the critical temperature is only twice as small as at zero temperature, and 
the gluon condensate is far from being melted down. This is another manifestation of the fact 
that the gluon condensate is not an order parameter of the chiral restoration phase transition. 
Neither it is an order parameter of the deconfinement phase transition in pure Yang-Mills 
theory.) 

Let us look first at our doped semiconductor when we decrease the density of admixture. 
Below some critical density, the set of atoms with essential overlapping of wave functions 
does not form a connected network in 3-dimensional space anymore. The electrons can no 
longer travel far through this network, wave functions become localized, and the conductivity 
drastically falls down. This is called the percolation phase transition (see for detailed 
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discussion). 

Likewise, there is a critical temperature in QCD above which instantons and antiinstantons 
do not form anymore a connected cluster with an essential overlap of individual fermion zero 
modes (what overlap is "essential" and what is not is a numerical question. For condensed 
matter systems, computer estimates for the critical admixture density were performed long 
time ago. Corresponding calculations in thermal QCD have been done only recently |6^]). 
At high temperatures, few remaining quasi-particles tend to form "instanton- antiinstanton 
molecules". These molecules tend to be "polarized" in imaginary time direction (the overlap 
integrals and, correspondingly, the "binding energy" for a molecule polarized along imaginary 
time axis is larger as the zero modes involve an exponential suppression factor only in spatial, 
but not in Euclidean time direction). 

A characteristic Euclidean path integral configuration at high temperature (here "high" 
means just several hundred MeV) is drawn schematically in Fig.^. In this picture, the indi- 
vidual zero modes are not spread out uniformly in the range AA ~ fJ-hadr after diagonalization, 
as was the case at zero temperature where instantons and antiinstantons formed an infinite 
cluster, but are just shifted by the value ~ fJ-hadr due to interaction in individual molecules. 
Small eigenvalues in the spectrum of Dirac operator are absent and the fermion condensate is 
zero. 
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Figure 18: Gas of instanton-antiinstanton molecules (high T). 



It is worthwhile to emphasize once more that this scenario of percolation phase transition 
leading to the molecular high-temperature phase is expected to hold only at Nf > 2. For 
Nf = 1 with arbitrary small but nonzero fermion mass, molecules are effectively ionized even 
at very high temperatures and the "medium" presents a very dilute instanton-antiinstanton 
gas — the instanton density involves a product of two small factors: exp{— 87r^/<^^(T)} and 
the fermion mass m. Differentiating logZ over m and sending m to zero, one obtains a small 
but nonzero quark condensate H, 115]. Cf. the analogous situation in the Schwinger model 

lie 



^^Whcthcr molecules are ionized or not depends on the particular quantity we are interested in. An accurate 
statement is that, for Nf — 1, the gas component gives a non- vanishing contribution to the quark condensate in 
the chiral limit. Thermodynamic quantities like the energy density would not depend on the gas component at all 
in the limit m — > and would be determined by "molecular" topologically trivial configurations. For Nf = 2, the 
gas component does not contribute to the condensate in the cliiral limit, but contributes to the so called scalar 
susceptibility which is related to the double derivative of the free energy over mass. The susceptibility is small at 
high temperatures but never vanishes. In the theory with 3 massless flavors, also the susceptibility vanishes in the 
high temperature phase, but the expectation value < uu dd ss >t is still nonzero. Gas component is responsible for 
the explicit Ua{^) breaking via anomaly. If we are specifically interested in L'^(l) breaking effects, we always have 
to take it into account. Disregarding it can lead to wrong results |148|. 
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5.5 Disoriented chiral condensate. 



When we talked in previous sections about "experimental" tests of theoretical predictions, we 
meant numerical lattice experiment. However, speaking of the particular problem of the phase 
transition in QCD associated with chiral symmetry restoration, an intriguing possibility exists 
that a direct experimental evidence for such a transition can be obtained at the high-energy 
heavy ion collider RHIC which is now under construction. 

After a head-on collision of two energetic heavy nuclei, a high temperature hadron "soup" 
is created. We do not call this soup the quark-gluon plasma because, even at RHIC energies, 
the temperature Trhic ~ 0.5 GeV Q would not be high enough to provide a sufficient small- 
ness of the effective coupling g^{T) and to make the perturbation theory over this parameter 
meaningful. (As will be discussed in details in Chapter 6, a proper normalization point for the 
coupling constant is probably 2tiT rather than T. Still perturbative corrections to the physical 
quantities like the free energy density appear to be rather large at T ~ .5 GeV, and the per- 
turbative series does not display a sign of convergence.) What is important, however, is that, 
at RHIC energies, the temperature of the soup would be well above the estimate ( ^.18 ) for the 



phase transition temperature. The high-temperature state created in heavy nuclei collision 
would exist for a very short time, after which it expands, cools down and decays eventually 
into mesons. 

Let us look in more details at the cooling stage. At high temperature, the fermion conden- 
sate is zero. Below phase transition, it is formed and breaks spontaneously chiral symmetry. 
This breaking means that the vacuum state is not invariant under the chiral transformations 
( |4.2| ) and a direction in isotopic space is distinguished. What particular direction — is a matter 
of chance. This direction is specified by the condensate matrix 

^ff'=<qLfqRf'> (5.28) 

For simplicity, we have assumed up to now that the condensate matrix is diagonal '^jf = 
— ^Sff- But any unitary matrix can be substituted for Sff (of course, it can be brought 
back in the form 6ff' by a chiral transformation ). In different regions of space, cooling occurs 
independently and flavor directions of condensate are not correlated. As a result, domains 



with different directions of condensate shown in Fig. 19 are formed (cf. cooling down of a 
ferromagnetic below the Curie point). 
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Figure 19: Domains of disoriented chiral condensate in cooling hadron soup. 



In our world, we do not observe any domains, however. The direction of the condensate 
in all spatial points is identical. This is a consequence of the fact that u- and d- quarks have 
nonzero masses which break chiral symmetry explicitly, the vacuum energy involves a term 

Evac ~ Tr{A4"^S} + c.c. (5.29) 



and the only true vacuum state is ( |4.5| ) (in the basis where the quark mass matrix Ai is 
diagonal) . 

However, the masses of u- and d- quarks are rather small and one can expect that the 
domains with "wrong" direction of the condensate are sufficiently developed during the cooling 
stage before they eventually decay into true vacuum (^]^) with emission of pions. |^ 

This is a crucial assumption. A theoretic estimate of the characteristic size of domains 
they reach before decaying is very difficult and there is no unique opinion on this question 
in the literature. But if this assumption is true, we can expect to observe a very beautiful 



effect 1 149, |123[ | (for a recent review, see |137| ]). From the true vacuum viewpoint, a domain 
with disoriented Sjj/ is a classical object — kind of a "soliton" (quotation marks are put here 
because it is not stable) presenting a coherent superposition of many pions. The mass of this 
quasi-soliton is much larger than the pion mass. The existence of such multi- pion coherent 



states was discussed long ago in pioneer papers [ 150 | but not in relation with thermal phase 
transition. 

Eventually, these objects decay into pions. Some of the latter are neutral and some are 
charged. As all isotopic orientations of the condensate in the domains are equally probable, 
the average fractions of vr*^, vr"'', and vr" are equal: < /^o >=< >= | as is also the case 



■^°Fig.|T^ implies the existence of several domains and describes better the physical situation immediately after the 
"phase crossover" in early Universe. Probably, the size of the hot fireball produced in collision of two nuclei is too 
small and the cooling occurs too fast for several domains to be developed. The popular "Baked Alaska" scenario 



1 149 1 implies the formation of only one domain with (generally) wrong fiavor orientation. 
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for incoherent production of pions in, say, pp collisions where no thermalized high-T hadron 
soup is created. 

But the distribution P{f) over the fraction of, say, neutral pions is quite different in the case 
of incoherent and coherent production. In incoherent case, P{f) is a very narrow Poissonic 
distribution with the central value < f^^o >= 1/3. The events with f^o = or with f^o = 1 
are highly unprobable: P{0) ~ ^(1) ~ exp{— CA^} where ^ 1 is the total number of pions 
produced. 

For coherent production, the picture is quite different, '^ff is proportional to a SU{2) 
matrix. Factorizing over ?7(1), one can define a unit vector in isotopic space € S"^. The fraction 
of vr^ produced would be just / = cos^ 6 where ^ is a polar angle on 5^. The probability to 
have a particular polar angle normalized in the interval < ^ < 7r/2 [ the angles Q > tx jl 
do not bring about anything new as /(vr — 0) = f{6) ] is P{9) = sm.0. After an elementary 
transformation, we get a normalized probability in terms of /: 



p{m = 




(5.30) 



As earlier, < / >= 1/3, but the distribution in / is now wide and the values / = and / = 1 
are quite probable. 

Thus a hope exists that, in experiments with heavy ion collisions at RHIC, wild fluctuations 
in the fractions of neutral and charged pions would be observed. That would be a direct 
experimental indication that a quasi-phase-transition occurs where domains of disoriented 
chiral condensate of noticeable size are developed in a cooling stage. One can recall in this 
respect mysterious Centauro events with anomalously large fraction of neutral or of charged 
particles observed in cosmic ray experiments [ [151| . Cosmic ray experiment is not a heavy ion 
collision experiment and there are no strong theoretical reasons to expect the formation of 
hot hadron matter there (see, however, [|149| ), but who knows, may be that still was the first 
experimental observation of the QCD phase transition ? 



6 Quark— gluon plasma. 

The last chapter of the review is devoted to the properties of high-temperature phase, the 
quark-gluon plasma. It is a very rich and interesting physical system which attracted lately 
a considerable attention of theorists. The number of the papers in SLAC archive involving 
the words "quark- gluon plasma" in the title exceeds seven hundred (not yet exceeds right 
now when I write it, but will certainly pass this mark by the time you read it). Thereby, this 
chapter is longer than others and still many issues were left undiscussed. In particular, we 
will touch only superficially a fascinating, but technical subject of effective lagrangian of soft 
modes and of nonabelian Vlasov equations. Instead, we will discuss at some length different 
physical phenomena which are specific for QGP with an emphasize on the problem of physical 
observability of various characteristics of QGP discussed in the literature. 
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6.1 Static Properties of QGP: a bird eye's view. 



We start with discussing static (thermodynamic) characteristics of QGP and find it useful to 
do it in two "rounds". In this section, we describe gross physical features of the system and 
go down to details in the next one. 

Thermodynamics 

The basic thermodynamic characteristic of a finite T system is its free energy. The lowest 
order results are derived exactly in the same way as for free pion gas (see Sect. 4.2). We have 
for a pure gluon system 



fP9 

— = 2{Nl - l)T 
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(6.1) 



where 2{N'^ — 1) is the number of degrees of freedom of the gluon field (the factor 2 comes 
due to two polarizations). This is nothing else as the Stefan-Boltzmann formula multiplied by 
the color factor N"^ — 1. The quark contribution is obtained quite similarly. Taking the Pauli 
principle into account, we have instead of Eq.( 4.2C| ) 



pi 



-TlnZ 



-Tin 



-Mp\ 



ANcNf 



(6.2) 



where ANcNf is the number of degrees of freedom. Trading the sum for the integral, we obtain 
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Notice that in our world with Nc 



N 



f 



3, the lowest order quark contribution to the free 



energy is roughly 2 times larger than the gluon one. 

All other thermodynamic quantities of interest can be derived from the free energy by 
standard thermodynamic relations. For example, the pressure just coincides with the free 
energy with the sign reversed. The energy density is 



E 



dT 



(6.4) 



For massless particles in the lowest order the relation 

E = 3P = -3F (6.5) 

holds. 

Debye screening. 

Consider the gluon polarization operator in QGP with account of thermal loop corrections. 
It is transverse, /c^n^,y(A;) = 0. At zero temperature, transversality and Lorentz-invariance 
dictate the form n^iy(/c) = n(/c^)(5f^,y — k^ky/k"^). At finite T, Lorentz-invariance is lost and 
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the polarization operator presents a combination of two different (transverse and longitudinal) 
tensor structures. Generally, one can write 



Uoo = Ui{io,\k\) 

n,o = ^nK..,|k|) 

n,, = -Uticu, \kmj - kikj/k^) + "^^mcu, |k|) (6.6) 

Consider first the longitudinal part of the polarization operator in the kinematic region where 
CO is set to zero in the first place after which k = \k\ is also sent to zero. By the reasons which 
will be shortly seen, we denote this quantity 

ml = limUi{0,k) (6.7) 

To understand the physical meaning of this quantity, consider the correlator 

< ^g(x)^^(0) >~ 5"'' / ^ e'*^^ oc e-^^l'^l (6.8) 

Thus niD coincides with the inverse screening length of chromoelectric potential Aq. 

There is a clear analog with the usual plasma. A static electric charge immersed in the 
plasma is screened by the cloud of ions and electrons so that the potential falls down expo- 
nentially oc exp{— r/r£)} where the Debye radius in an ordinary nonrelativistic plasma is 
given by the expression (see e.g. p^) 

^.,^W^4™(Z^ (0.9) 

where n is the electron and ion density, Z is the ion charge (so that the second term describes 
the ion contribution), and, to avoid unnecessary complications, we assumed that the electron 
and ion components of the plasma have the same temperature T Calculating the one-loop 
thermal contribution to the gluon polarization operator (see Fig.^), one can easily obtain an 
analogous formula for QGP: 

ml = ^d' = ^ + ^) (6.10) 

where the first term describes the screening due to thermal gluons and the second term — 
the screening due to thermal quarks. The result ( |6.10|) was first obtained by Shuryak | 152 |. It 



has the same structure as (^) (Note that the density of particles in ultrarelativistic plasma 
is expressed via the temperature, n cx T^). It is worth mentioning that the Debye screening 
is essentially a classical effect and not only quarks but also gluons induce screening rather 
than antiscreening. That should be confronted with the famous antiscreening of the charge in 
Yang-Mills theory at zero temperature due to quantum effects. 
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p 



a) b) c) d) 

Figure 20: Gluon polarization operator in one loop. 



In QED, the correlator < Aq{x)Aq{0) >t is a gauge invariant object and also the notion of 
charge screening is unambiguous and well defined — there are classical electric charges and one 
can measure the potential of such a charge immersed in plasma by standard classical devices. 
Not so in QCD. We do not have at our disposal classical color charges due to confinement and 
the experiment of measuring the chromoelectric field of a test color charge cannot be carried 
out even in principle. Also the gluon polarization operator which enters the definition (|6.8| ) of 
a Debye screening mass is generally speaking gauge-dependent. Thus a question of whether 
the Debye screening mass is a physical notion in a nonabelian theory is fully legitimate. 

The answer to this question is positive. But with reservations. 

Note first of all that though the gluon polarization operator is generally gauge-dependent, 
the Debye screening mass defined in ( |6.7D does not depend on the gauge in the leading order 
(we shall discuss also nonleading corrections in the next section) . That suggests that the result 
( |6.10 ) has an invariant physical meaning. 



Indeed, one can consider a gauge-invariant correlator 

C(x) =< P(x)P*(0) >T (6.11) 

where 

P(x) = ^Trexp|i^''ylo(x,T)(ir| (6.12) 



is the Polyakov line [both the Polyakov line and the correlator ( |6.11 ) were discussed at length 



in Chapter 3]. If r = |x| is not too large (the exact meaning of this will be specified later), the 
connected part of the correlator Cc(r) = C(r) — C(oo) (as was explained in Chapter 3, C(c«) 
is a nonzero constant) has an exponential behavior 

Ce(r) oc [Goo(r)]2 oc ^ex^{-2mDr] (6.13) 

The correlator ( |6.1l| ) can be attributed a physical meaning. — TlnC(r) coincides with the 
change of free energy of the system when putting there a pair of heavy quark and antiquark 
at distance r, the averaging over color spin orientations being assumed p^]. 
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Also we shall see in the next section that the perturbative corrections ^ g^T^ to a perfectly 
well defined and physical quantity — the free energy — depend directly on the value of m^i. 

Thus the notion of Debye screening mass in QGP is physical, though, unfortunately, not to 
the same extent as it is in the usual plasma. The correlator ( 6.11| ) with the correlation length 



of fractions of a fermi cannot be directly measured even in a gedanken laboratory experiment 
— at least, we cannot contemplate such an experiment. But it can be measured in lattice 
numerical experiments which is almost as good. We shall see later that a number of other 
characteristics of QGP have a similar semi- physical status. 

Magnetic screening. 

Debye mass describes the screening of static electric fields. In abelian plasma, magnetic 
fields are not screened whatsoever. Such a screening could be provided by magnetic monopoles, 
but they are not abundant in Nature. The absence of monopoles is technically related to one 
of the Maxwell equations diBi = 0. In the nonabelian case, the corresponding equation reads 
"DiBi = where 1? is a covariant derivative. Thus gluon field configurations with local color 
magnetic charge density pm ~ d^Bi (with usual derivative) are admissible. The presence of 
such configurations in the gluon heat bath results in screening of chromomagnetic fields. 

Let us see how it comes out in a perturbative calculation. Consider the one-loop graph in 
Fig.pO^ which contributes to the polarization operator of the spatial components of the gluon 
field k) oc 11^(0, k). At high temperatures, it suffices to take into account only the lowest 

Matsubara frequency (with w = 0) in gluon propagators. We have 

Ut{0, k) ~ g^k^T J p,^^'^^^, ~ g^Tk (6.14) 

The numerical coefficient can be explicitly calculated, but it depends on the gauge and makes 
as such a little sense. The loop integral is determined by the low momentum region |p|char ~ 
k <ti T (as we are interested in the large distance behavior of the gluon Green's function, we 
keep k small). Note that the Feynman integral for 11^(0, k) has a completely different behavior 
being saturated (in the leading order) by the loop momenta |p| ~ T — it is a so called hard 
thermal loop. 

The transverse part of the gluon Green's function is 

We see that at ~ g^T the one-loop contribution to the polarization operator is of the same 
order as the tree term k'^. One can estimate a two- loop contribution to llt{0,k) which is of 
order (g'^T)'^. The factor here comes from two loops (we use the rule ( |2.7| ) and take into 
account only the contribution of the lowest Matsubara frequency in each loop). Similarly, a 
three-loop contribution to n((0,A;) is of order {g'^T)^/k, a four-loop contribution is of order 
{g'^T)^/k'^ etc. (The growing powers of k in the denominator are provided by infrared 3- 
dimensional loop integrals. Infrared integrals may also provide for a logarithmic singularity in 
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external momentum k in the two- loop and higher loop contributions, but weak logarithmic 
factors are of no concern for us here). At A; ~ g'^T, all contributions are of the same order and 
the perturbation theory breaks down (3^]. 

There is an alternative way to see it. Let us write down the expression for the partition 
function of the pure glue system at finite temperature: 

Z = J ndA-(x,r)exp|-^ J^T J d^C)'} (6-16) 

When T is large and P small, the Euclidean time dependence of the fields may be disregarded. 
Also the effects due to Aq{x,t) can be disregarded — time components of the gluon field 
acquire the large mass ruD ^ gT ^ g'^T and decouple. We are left with the expression 

Z = j J]d^^(x)exp|-^ j I (6.17) 



A theory with quarks is also reduced to ( 6.17 ) in this limit — the fermions have high Matsubara 



frequencies ~ T and decouple. The partition function ( |6.17 ) describes a nonlinear 3D theory 



with the dimensional coupling constant g^ ~ g'^T No perturbative calculations in this 
theory are possible. 

In QED, there is no magnetic screening (and the effective 3-dimensional theory is trivial). 
Presumably, perturbative calculations in hot QED can be carried out at any order, though 
this question is not absolutely clear. 

What is the physical meaning of this nonabelian screening ? Can it actually be measured 

? 

The way we derived it, the magnetic screening shows up in the large distance behavior of 
the spatial gluon propagator Gt{r) oc exp{ik*r} where k* is the solution to the equation 

k'^ + Ut{0,k) =0 (6.18) 

There is no reason to expect that A: = is a solution — the behavior of nj(0, k) in the limit 
A; ^ where the perturbation theory does not work is not known, but one can tentatively 
guess that it tends to a constant ~ (g'^T)^ (As we have seen, g^T is the only relevant scale 
in this limit). If so, k* ~ ig^T and Gt{r) falls down exponentially oc exp{— Cg^Tr} at large 
distances. Hence the term "magnetic screening". 

Unfortunately, the gluon polarization operator is a gauge- dependent quantity and, would 
the God provide us with the exact expression for nj(0, /c) in any gauge, even He could not 
guarantee that the solution of the dispersive equation (|6.18| ) would be gauge-independent 0. 
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One should be careful here. It would be wrong to use the expression (6.17) for calculating, say, the free energy 



of QGP at large T. The latter takes the contributions from hard thermal loops (involving also fermions !) with 



momenta of order T which are not taken into account in Eq. (6.17). Eq. (6.17) should be understood as an effective 
theory describing soft gluon modes with momenta of order g^T . 
^^cf. Eq.( 6.74| ) and the related discussion in sect. 6.4. 
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What one can do, however, is to consider the correlator of chromomagnetic fields Cafc(x) = 
< B"(x)B^(0) >T- This correlator also depends on the gauge in a nonabelian theory, but the 
gauge dependence amounts to rotation in color space: Ca6(x) Cla'ai^)C"' ^ (x)J7{,/b(0) and 
cannot affect the exponential behavior of the propagator. This is the way the magnetic mass 
is usually measured on the lattices. 

One can do even better considering gauge-invariant correlators, the simplest one is < 
G^,^(x)G^y(0) >T- Then the quantity 

/i = -^h^ Jin < G2,(x)G2,(0) >t (6.19) 

provides an invariant definition of the magnetic screening mass. The indefinite article is crucial 
here. Choosing other gauge-invariant correlators, one would get other invariant definitions. 
For example (and this will be important in the following discussion), the true correlation length 
of the correlator of Polyakov loops ( 6.11| ) is also of order {g'^T)~^. The asymptotics ( |6.13D is 



an intermediate one and holds only in the region 

{gT)-^ « r « (g^T)-^ (6.20) 

What was wrong with our previous derivation ? The matter is we took into account earlier 
only the thermal corrections to the gluon propagator and tacitly neglected corrections to the 
vertices. Ab accurate analysis |27| shows that this is justified in the range (la2C| ) but not 



beyond. An example of the graph providing the leading contribution in < P(x)P*(0) >t at 
r > {g'^T)-^ is given in Fig.^ 



The asymptotics of the correlator of Polyakov loops at large enough r is determined by the magnetic photon 
exchange also in abelian plasma. Two magnetic photons can be coupled to P{x) via a fermion loop. In abelian case, 
magnetic photons are massless and, as a result, the correlator has a power asymptotics cx [153|. Physically, 

it corresponds to Van-der-Vaals repulsion between the clouds of virtual electrons and positrons formed near two 
heavy probe charges. Seemingly, the asymptotics cx for the Polyakov lines correlator in QGP found in |154| has 
this origin. But in nonabelian case, the asymptotics becomes exponential when taking into account the magnetic 
screening effects. 
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Figure 21: Correlator of Polyakov loops at large distances. Wavy lines are electric gluons and 
dashed lines are magnetic gluons. 



Thus the "experimental status" of the magnetic mass ( |6.1£ ) and of other similar quantities 



is roughly the same and even better than for the Debye mass. The exponential behavior 
oc exp{— /ir} of the correlators is expected to hold at any large r and the value fj, can in 
principle be determined in lattice experiments with any desired accuracy. On the other hand, 
the Debye screening mass cannot be determined with an arbitrary accuracy due to the finite 
range ( |6.20D where the asymptotics ( |6.13 ) holds (being sophisticated enough, an invariant 



definition of Debye mass still can be suggested for A'^c > 3 — see the discussion in the following 
section). And, as was also the case for the Debye mass, we cannot invent any laboratory 
experiment where the magnetic mass could be directly measured. 

Speaking of the lattice experiment, it can be done and actually has been done very recently 
by the Bielefeld group. Electric and magnetic screening lengths were measured in the high 
temperature region up to T ~ lO^GeV. Preliminary data [155| show a slow falloff of the 



quantities rriu/T and fimagn/T with temperature in a qualitative agreement with the expected 
behavior mo/T ~ g{T) and Umagn/T ~ 5^(T). 

6.2 Static Properties of QGP: perturbative Corrections. 

Debye mass. 

To provide a smooth continuation of the discussion started at the end of the previous section, 
we consider first higher- loop effects in the Debye mass. Note first of all that the definition ( |6.7] ) 
is not suitable anymore. Higher- loop corrections ArriD ~ g'^T to the Debye mass as defined in 
Eg. (|6?^) depend on the gauge ||156|] . A better way is to define the Debye mass as the solution 
of the dispersive equation k"^ + 11/(0, k) = 0. In other words, we define 



mj)= lim Ui{0,k) (6.21) 
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The longitudinal part of the gluon Green's function 

Gi{0,k) 



has then the pole at k = imij. Pole position has the "tendency" to remain gauge- invariant 
even though the Green's function itself is not. We will see in the following section that there 
are cases when formal arguments displaying gauge-invariance of the pole position fail in higher 
orders of perturbation theory due to severe infrared singularities. Speaking of Debye mass as 
defined in Eq.( |6.21 ), its gauge invariance can be shown in the order ~ g^T\n{C/g), but the 
constant under the logarithm is not computable in perturbation theory in principle, and it is 
difficult to pose a question whether it is gauge-invariant or not. In the next-to-leading order, 
the (gauge- invariant) result is [157] 

NcQ^T 
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(6.22) 



with logarithmic accuracy. As was mentioned, the constant C is not computable by per- 
turbation theory. We see that the correction is non-analytic in coupling constant. The 
non-analyticity appears due to bad infrared behavior of the loop integrals — they involve 
a logarithmic infrared divergence and depend on the low momenta cutoff which is of order of 
magnetic mass scale g^T. Thus we have 
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The magnetic infrared cutoff is actually provided by higher- order graphs — the orders of 
perturbation theory are mixed up and a pure two-loop calculation is not self-consistent. As 
we have seen, ^mag cannon be determined analytically which means that the correction ~ g^T 
without the logarithmic factor in the Debye mass cannot be determined analytically. 

Also the correction ~ g^T cannot be "experimentally observed". Really, we have seen 
that the invariant physical definition of uid refers to the correlator of Polyakov loops ( |6.11| ). 
The latter displays the exponential behavior (6.13) in the limited range ( |6.20 ). But that is 
tantamount to saying that the correction ~ g^T in the Debye mass cannot be determined 
from the correlator (6.11). To do this, one should probe the distances r ^ (g'^T)~^ where the 



correlator has a completely different behavior being determined by the magnetic scale. The 



correction (6.22) is still observable, however, due to a logarithmic enhancement factor. 
To be more precise, the correlator in the range ( |6.20| ) has the form [27| 
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(6.23) 



where m 



is the lowest order Debye mass ( 6.10 ). 



It is not just 



Cc(r) oc exp{— 2m£ir} 



(6.24) 
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, with m£) being defined in ( |6.22 ), by two reasons. First, at finite r <^ {g^T) ^, the infrared 
cutoff in the loop integrals is provided by rather than iXmag and we have 



In- 



In 



ln(m£)r) 



(6.25) 



Second, there are two kinds of graphs involving the corrections both in the gluon propagator 
and in the vertex (see Figj2^). 
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a) 



b) 



Figure 22: Correlator of Polyakov loops in next-to-leading order. Wavy lines are electric gluons 



with Debye mass (6.10) taken into account. Dashed lines are magnetic gluons. 



The graph in Fig.|22|a would give the correction 



1 In(mz)r) (6.26) 



with logarithmic accuracy. It corresponds to Eq.( |6.24 ) where the substitution ( |6.25| ) is made. 
The second graph also contributes in this order and gives a positive correction with a twice as 
small coefficient | p7| . Exponentiating the one loop correction we arrive at the result ( |6.23| ). 
Note that the Fourier image of the correlator (6.23) does not have a singularity at finite k 
whatsoever. Thus the Debye mass pole in the gluon propagator does not show up as a pole 
in the gauge- invariant correlator of Polyakov loops. Thereby, the quantity ( |6.22 ) is not quite 
physical in spite of its gauge invar iance. 

But, anyway, in the theoretical limit ln[l/(/(T)] — > oo, there is a narrow range of r where 
the condition r <C {g^T)~^ is fulfilled so that the intermediate asymptotic law ( |6.23D holds 
whereas the correction in the exponent oc gm~^ r\ii{mDf) is of order 1 and can in principle be 
singled out in numerical lattice experiment. 

There are two interesting recent proposals to define Debye mass non-perturbatively in a 
gauge-invariant way |153f| . The first one is to consider the correlator of imaginary parts of 
Polyakov lines in a complex representation. If we are interested in a pure Yang-Mills theory, 
an additional requirement for the representation to be invariant under the action of the center 
of the group (such is, say, the decouplet representation in SU{'i)) should be imposed. ^ 
The point is that Pfg ~ -fio is odd under Euclidean time reversion and cannot be coupled to 



^'^ That the leading corrections ^ [g^Tr In(TODr)]" exponentiate is an educated guess. To the best of our knowledge, 
nobody has checked it so far explicitly. 

^^The authors of Ref. |l53|] argued the necessity to consider only the Zjy - invariant representations saying that 
otherwise the correlators give zero after averaging over different "Z^r - phases" . Actually, as we have seen in Chapter 
3, such phases do not exist in Nature and a proper argumentation should be that the correlators which are not 
invariant under Zn transformations just do not have a physical meaning p^ . 



92 



only magnetic gluons. That means that the correlator 



<Pi*o(x)-Pio(x), Pi*o(0)-i^io(0) >T (6.27) 

exhibits the Debye screening falloff even at arbitrary large r. In the intermediate region 
(6.20), the correlator ( |6.27| ) is described by the graph with exchange of 3 electric gluons and 



falls down as oc exp{— Sm/jr}. At very large distances, the dominant graph corresponds to 
the exchange of one electric gluon which is odd under the time reversion + a cloud of light 
magnetic gluons which neutralize the color. Debye mass extracted from the exponential falloff 
of such a correlator at large distances has a gauge invariant meaning, but depends, generally 
speaking, on the choice of the representation (probably, not in the order g'^T ln{l / g)) . Also 
this definition does not work for SU{2) where the Polyakov line is real in any representation. 

Another suggestion was to study the behavior of large spatial Wilson loops in adjoint color 
representation. At large distances, the theory is effectively reduced to the 3- dimensional YM 
theory ( 6.17| ). Adjoint color charges in this theory are not confined but rather screened, and 



the Wilson loop exhibits the perimeter law behavior 

W{C) oc exp{— m* x perimeter(C)} (6.28) 
One can show that m* is of order g^ ln{l/ g)T and coincides with a perturbative correction 



to the Debye mass in (6.22). The problem here is that m* has no trace of the lowest order 



contribution ( |6.10 ). Still, m* certainly has an invariant meaning and is as such an interesting 



quantity to study. 

Free energy. 

This is the most basic and physical quantity of all. May be this is the reason why pertur- 
bative corrections are known here with record precision. 
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Figure 23: Free energy in two loops. 

The correction ~ g^T'^ has been found by Shuryak |152| . To this end, one should calculate 
the two-loop graphs depicted in Fig.^. The behavior of Feynman integrals in this order is 
quite benign and no particular problem arises. 
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Figure 24: An infrared- divergent 3-loop graph in free energy. 



However, starting from the three-loop level, the problems crop up. The contribution of the 



graph depicted in Fig. 24 is infrared divergent: 



24 



/T^/f (6.29) 



That means that a pure 3-loop calculation is not self- consistent and, to get a finite answer, 
one has to resum a set of infrared-divergent graphs in all orders of perturbation theory. This 



is, however, not a hopeless problem, and it has been solved by Kapusta back in 1979 |15J] 
The leading infrared singularity is due to the so called "ring diagrams" depicted in Fig. 25. 
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Figure 25: Ring diagrams. Grey circles stand for a set of one-loop graphs in Fig. 



,20 



The sum of the whole set of ring diagrams has the form 



1 + 



fc2 



fc2 



(6.30) 



The expansion of the integrand in 11^(0, A:) //c^ restores the original infrared-divergent in- 
tegrals. However, the whole integral is convergent being saturated by momenta of order 
k ~ [IT; (0,0)]^/^ ~ gT. Neglecting the /^-dependence in the polarization operator, we ob- 
tain 



T / d-^k 



m 



D 

A;2 



(6.31) 



We see that the correction is non-analytic in the coupling constant which exactly reflects 
the fact that the individual graphs diverge and the orders of perturbation theory are mixed 
up. Note, however, that infrared divergences here are of comparatively benign variety — the 
integral depends on the scale k'^^"'^ ~ gT and we are far from the land of no return k ~ g^T. 
That is why an analytic determination of the coefficient in ( |6.31 ) is possible. 

The next correction has the order ~ g'^T'^ ln{l / g) . It comes from the same ring graphs of 
Fig. 25 where now the term ~ g'^Tk in 11/(0, k) should be taken into account. This correction 
was determined by Toimela [ |159[ |. 

Presently, the correction ~ g^T'^ (without the logarithmic factor) |160|] and the correction 



g T 1 161, 162] are known. This is the absolute limit beyond which no perturbative cal- 



culation is possible — similar ring graphs as in Fig.^ but with magnetic gluons would give 
the contribution oc T^f^^g • 
determined perturbatively. 



the contribution oc T^f^^g ~ g^T^ in the free energy. And, as we have seen, fimag cannot be 
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Collecting all the terms, the following result is obtained [161, 162 



where 



F 



45 



IT 



5/2 



21 



32 



(6.32) 



(6.33) 



^ 4 V 12 ^ 



(6.34) 



F3 = 30 1 + -^ 



(6.35) 
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F4 = 237.2 + 15.97A^/ - 0.413iV| + ^ ( 1 + ^ ) In 



1 + —Nf 
12 ^ 



1 



vr V 6 



33 ^ 



In 



27rr 



(6.36) 



1 + 



6 



-799.2 - 21.96iV/ - 1.9267V| 



+- 



485 



Nf 



1 Nf In 

33 



27rT 



(6.37) 



The expressions ( |6.33| ) - ( |6.37| ) are written for Nc = 3. The coefficients like 237.2 are not a 
result of numerical integration but are expressed via certain special functions. 

A nice feature of the result ( |6.32| ) is its renorm- invariance. The coefficients F4 and F5 
involve a logarithmic //-dependence in such a way that the whole sum does not depend on the 
renormalization scale /i. 



Let us choose fi = 2ttT (this is a natural choice |99|, 2ttT being the lowest nonzero gluon 
Matsubara frequency ) and Nf = 3. In that case, we have 



1 - 0.9a^ + 3.3a^/2 ^ (71 ^ 3.51na^)a2 _ 20.8a^/^ 



(6.38) 



5/2 

Note a large numerical coefficient at as . It is rather troublesome because the correction 

5/2 

~ as overshoots all previous terms up to very high temperatures and, at temperatures 
which can be realistically ever reached at accelerators, makes the whole perturbative approach 
problematic. 
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Take T ~ 0.5 GeV (as was mentioned earlier, this is the temperature one can hope to 
achieve at RHIC [^). Then 27rT ~ 3 GeV and ~ 0.2. (We use a conservative estimate for 



as following from T physics | 163 |. Recent measurements at LEP favor even larger values.). 



The series ( 6.38| ) takes the form 

F = Fo[l- .18 + .3 + .06 - .37 + . . .] (6.39) 

5/2 

which is rather unsatisfactory. We emphasize that the coefficient of as is rather trustworthy 
being obtained independently by two different groups. 



The estimate (|6.39f ) was obtained in the assumption Nf = 3. In pure Yang-Mills theory, 
the behavior of perturbative series is even worse 

1 - 1.19a^ + 5.390^/2 + (16.20 + 6.841nas)a^ - 45.69a^/2j (6.40) 

which means that perturbation theory cannot be trusted at all at ^ .15. One could naively 
expect that deviations of real free energy density from Stephan-Boltzmann value would be 
considerable up to very large temperatures. Somewhat surprisingly, numerical lattice calcu- 
lations show it is not the case. According to recent measurements in pure SU{3) Yang-Mills 



theory |164], the free energy density at T = 5Tc constitutes 85% of the Stephan-Boltzmann 
value. This suggests that the coupling constant in this region is small. Indeed, the estimate 
for the coupling constant obtained in |164| | from measurement of another lattice quantity, the 
spatial string tension, is as ~ .12 at T = 5Tc. If substituting it in the series ( |6.40| ), we find 
F/Fq ~ .87 in a perfect agreement with the data. Actually, the smallness of the coupling 
constant here is an "optical illusion". When comparing the pure Yang- Mills theory with 
QCD with quarks, one should compare the scale of perturbative corrections which is of order 
~ asCv in the pure Yang-Mills theory versus ~ agCp in QCD. This scale is roughly the same 
in pure YM theory and in QCD for the same temperature (measured in units of string tension 
a or in units of Tc). 

At lower temperatures, perturbation theory does not work and non-perturbative effects 
(associated with instantons) become important. The lattice values of the ratio F/Fq rapidly 
decrease with temperature and are rather small at T ~ Tc. 



The last remark is technical. The result ( 6.32 ) was obtained in Euclidean technique. There 
is a real time calculation which correctly reproduces the two-loop term ~ (7^T^ in free energy 



1 17], but nobody so far succeeded in calculating in this way the terms ~ g^T'^ and higher. 



Certainly, real-time technique is not very suitable for calculation of static quantities, and one 
way to get the result is good enough, but, to my mind, it is an interesting methodical problem. 

6.3 Collective excitations. 

One of striking and distinct physical phenomena characteristic of usual plasma is a non-trivial 
dispersive behavior of electromagnetic waves. In contrast to the vacuum case where only trans- 



36 



I am indebted to F. Karsch who pointed my attention to this fact. 
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verse photons with the dispersive law a; = |k| propagate, two different branches with different 
non-trivial dispersive laws uj±{k) and t^\\{k) appear in plasma. The value u;_l(0) = ti^j|(0) = iVpi 
characterizes the eigenfrequency of spatially homogeneous charge density oscillations and is 
called the plasma frequency. 

A similar phenomenon exists also in QGP. The spectrum of QGP involves collective excita- 
tions with quantum number of quarks and gluons. Like in usual plasma, there are transverse 
and longitudinal branches of gluon collective excitations (alias, transverse and longitudinal 
plasmons) and their properties on the one-loop level are very similar to the properties of pho- 
ton collective excitations in usual plasma. A novel feature is the appearance of non-trivial 
fermion collective excitations (plasminos) . But, again, they are not specific for a nonabelian 
theory and appear also in ultrarelativistic e^e~ plasma (in the limit when the electron mass 
can be neglected compared to the temperature- induced gap in the electron spectrum oc eT). 

One-loop calculations. 

The dispersive laws of quark and gluon collective excitations can be obtained via solution of 
a nonabelian analog of the Vlasov system involving the classical field equations in the medium 
and Boltzmann kinetic equation |165]. This way of derivation makes analogies with usual 



plasma (where the Vlasov system is a standard technique) the most transparent. 

I will outline here another way of derivation which is more conventional and more easy to 
understand for a field theorist. This is actually the way the results were originally derived 
I I6t , |167]] (Note, however, that analogous results for abelian ultrarelativistic plasma were first 



derived long time ago in Vlasov technique |16^ ]). 

Consider the gluon Green's function in a thermal medium. As was mentioned in sect. 2, 
at T 7^ 0, different kinds of Green's function exist. To be precise, we are considering now the 
retarded Green's function 

R 



i<e{t)[Al{x),Alm>T (6.41) 



which describes a response of the system on a small perturbation applied at t = at some 
later time i > 0. Q As was discussed in details in Sect. 2.2, the Fourier image of ( |6.41D is free 
of singularities in the upper oj half-plane. The poles of G^^{uj,'k) correspond to eigenmodes 
of the system and exactly give us the desired spectrum of gluon collective excitations. The 
dispersive equation 

det||[G^,(a;,A;)]-i||=0 (6.42) 

splits up in two: 

u? -k^ - Iit{oJ,k) = 



^'^We will meet Vlasov equations again and discuss them in a little more details in the end of this chapter. 

^*As in commonly used gauges oc (5°^, we shall suppress color indices in the following. We have retained them 
here just to make clear that we are dealing with a commutator of Heisenberg field operators, not with a commutator 
of classical color fields. 
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k'^ + Ui{uj,k) = 



(6.43) 



with Ili^ti^^,k) being defined in Eq. (|6.6| ). The solutions to the equations ( |6.43D give two 
branches of the spectrum. 

The explicit one-loop expressions for Ili{uj,k) and Ilt{uj,k) obtained by the calculation of 
the graphs in Fig|2^ in the limit u),k <^T are [ |166| , |167| 



nt(a;, k) = -w, 



nKa;,A:) =3c^2^[l-F(a;/A:)] 



3 2 



A;2 



+ 



A;2 



-F{uj/k) 



(6.44) 



where 



gT 



Nf 



(6.45) 



is the plasma frequency and 



^ , , X , X + 1 



(6.46) 
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Figure 26: Plasmon spectrum in one loop. 



The behavior of dispersive curves is schematicahy shown in Fig. 26. At A; = 0, a;_L(0) 
t<;|j(0) = ojpi. Then the two branches diverge: 



6 

Lol{k<^gT)=Lol + -k'' 



(6.47) 



At k ^ gT both branches tend to the vacuum dispersive law 



u;l{k:$>gT)=k^ + ^Loli 

2e 



00?, {k > gT) = k"^ 



1 + 4 exp 



Hi 



(6.48) 



We see that uj^^{k) approaches the hne uj = k exponentiahy fast. 

The dispersive laws for quark collective excitations are obtained from a similar analysis of 
the quark Green's function. One loop expression for the finite temperature contribution in the 
fermion polarization operator in the limit cj, /c <C T is 



S(a;,k) = ^o'^F{Lu/k) 



k k 



[F{u/k) - 1] 



where ujq is the plasmino frequency at zero momentum: 



-CF 



(6.49) 



(6.50) 



We see that the fermion polarization operator involves two tensor structures ^7'^ — k7 and 
+ k7 which gives rise to two dispersive branches. We will call the branch corresponding 

to the Lorentz-invariant structure "transverse" and the branch corresponding to the structure 
+ k7 — "longitudinal". These terms may be misleading in the fermion case because, 
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in contrast to the plasmons with photon or gluon quantum numbers, these branches are not 
associated with transverse and longitudinal field polarizations. Hence the quotation marks. 
But better names were not invented, and using the words "transverse" and "longitudinal" 
still makes a certain sense because the physical properties of "transverse" and "longitudinal" 
fermion branches are rather analogous to the physical properties of transverse and longitudinal 
gluon branches. 
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Figure 27: Plasmino spectrum in one loop. 

The pattern of the quark spectrum is shown in Fig.^. It is similar to the gluon spectrum 
with one important distinction — at A; ~ 0, ^^(A;) and W|j(A;) involve linear terms of opposite 
sign: 

k 

iolik <^gT)=ujQ-^ (6.51) 

Thus u;J|(/c) first goes down and reaches minimum at some k* . The group velocity of the 
longitudinal plasmino at this point is zero. At large A;, 

ojf{k > gT) = k^ + IujI (6.52) 

and tends from above to the line uj = k exponentially fast. 

Landau Damping. 

The quoted one-loop results for the dispersive laws of transverse plasmons and plasminos 
are gauge-invariant and stable with respect to higher-order corrections. The latter is not true, 



however, for longitudinal excitation branches 169 , 17C ]. We have seen that ujI ^°°^{k) tends to 



the line uj = k exponentially fast at A: ^ gT. This can be easily seen from the analysis of the 
dispersive equations for longitudinal branches which, in the limit k ^ gT, have the form 

wii + A: ~ In — ^ (6.53) 

" k uj\\ — k 

At A; ^ gT, the solution exists when the logarithm is large and — A: is exponentially small. 
The logarithmic factor in Eq.( |6.53 ) comes from the angular integral 

d9 



00 — k cos 9 



(6.54) 
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which diverges at uj = k. This colUnear divergence appears due to masslessness of quarks and 



gluons in the loop depicted in Fig.|20|. But quarks and gluons in QGP are not massless — 
their dispersive law acquires the gap ~ gT due to temperature effects. An accurate calcula- 
tion requires substituting in the loops the dressed propagators. As a result, the logarithmic 
divergence in the integral (|6.54| ) is cut off and the logarithmic factor in Eq.( |6.53| ) is modified: 

Dressing of propagators amounts to going beyond one-loop approximation. Strictly speaking, 
to be self-consistent, one should also take into account one-loop thermal corrections to the 
vertices (this procedure is known as resummation of hard thermal loops |171[] ), but in this 
particular case these corrections do not play an important role. What is important is the 
cutoff of the logarithmic collinear singularity due to effective temperature-induced masses. 
Substituting (lelSSj ) in (|]5|), we see that the new dispersive equation does not at all have 



solutions with real uj for large enough k. This fact can be given a natural physical explanation. 
When k is small compared to gT, there is no logarithmic factor in the dispersive equation, the 
modification ( |6.55 ) is irrelevant, and the dispersive law of longitudinal modes does not deviate 



from the one-loop result. Then the logarithm appears, the modification ( |6.55| ) starts playing 
a role, and, at some k** ~ gT, the longitudinal dispersive curve crosses the line lo = k. At 
this point the longitudinal polarization operator acquires the imaginary part due to Landau 
damping. 

In usual plasma. Landau damping is the process when propagating electromagnetic waves 
are "absorbed" by the electrons moving in plasma. In the language of quantum field theory, 
it is a 2 —> 1 process 

7* + e ^ e (6.56) 

In real time technique, that corresponds to a contribution to the imaginary part of the polar- 
ization operator so that both internal electron lines in the loop are placed on mass shell. At 
T = the standard Cutkosky rules imply positive energies of all particles in the direct channel, 
and the imaginary part appears only due to the decay 7* — > e+ + e~ . At T 7^ 0, Cutkosky 
rules are modified and both signs for energy are admissible. |^ Physically, that corresponds 
to the presence of real particles in the heat bath so that the process ( |6.56|) may go. 

Also in QGP imaginary parts of polarization operators may acquire contributions due to 
Landau damping. The corresponding processes are 

9* + 9^9, 9* + qiq) ^ qiq), q*{q*) + 9 ^ q{q), q*iq*) + qiq) ^ 9, (6.57) 

where g*,q*,q* are plasmon and plasmino collective excitations and g,q,q are the excitations 
with characteristic momenta of order temperature (in this kinematic region, the dispersive 
laws are roughly the same as for free quarks and gluons and the star superscript is redundant). 



■^^We already exploited this fact in Sect. 4. 3 while calculating the damping of pion collective excitations in low 
temperature phase. 
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The kinematic condition for the scattering processes ( |6.57 ) to go is that the frequency of 
collective excitations g* and q* would be less than their momentum. We have seen that the 
condition tj < A; is realized indeed for longitudinal plasmon and plasmino excitations starting 
from some k** ~ gT. At k > k** , the Landau damping switches on and the dispersion law 
acquires an imaginary part. The imaginary part rapidly grows and very soon becomes of order 
of the real part. Prom there on, it makes no sense to talk about propagating longitudinal 
modes anymore. The situation is the same as in usual plasma where longitudinal modes also 
become overdamped at large enough momenta and disappear from the physical spectrum p^] . 
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Figure 28: Plasmon spectrum. 




The true pattern of plasmon collective modes is shown in Fig.^. A similar picture holds 

for plasminos. 



One-loop dispersive curves are gauge-invariant. However, the question whether these curves 
are physically observable is, again, highly non-trivial. It is easy to measure explicitly photon 
dispersive curves in usual plasma — to this end, one should study propagating classical elec- 
tromagnetic waves, measure the electric charge density (say, by laser beams) as a function of 
time and spatial coordinates and determine thereby the frequency and the wave vector of the 
wave. 

But there is no such thing in Nature as classical gluon field due to confinement, and no 
classical device which would measure the color charge density exists. Even more obviously, 
quark fields (which have Grassmanian nature) cannot be treated classically. Hence, one cannot 
really measure the energy and momentum of propagating colored waves in a direct physical 
experiment. 

What one can measure are correlators of colorless (in the first place, electromagnetic) 
currents. Modification of dispersive laws affects these correlators and that can be observed. 
However, a colorless current always couples to a pair of colored particles and, as a result, 
physical correlators involve some integrals of quark and gluon Green's functions which are 
related to quark and gluon dispersive characteristics only in an indirect way. Also, thermal 
modification of vertices is as important here as modification of the Green's functions. 

However, there is one special point on the dispersive curves which can in principle be 
directly measured in experiment. This is the point on the longitudinal plasmino curve where 
its frequency acquires a minimal value and the group velocity turns to zero. Consider the 
problem of emission of relatively soft e^e~ or fj.'^ pairs by QGP. In a "thermos bottle" 



Observability. 
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gedanken experiment, one should make sure that the size of your thermos bottle is much less 
than the lepton mean free path. Otherwise, the leptons are thermalized and their spectrum is 
just Planckian. Btt in heavy ion collisions experiments, QGP is produced in small volume, the 
condition Lchar < path satisfied, and the spectrum of emitted leptons (and photons) 
can provide a non-trivial information on dynamic characteristics of QGP. 

The spectrum of soft dileptons was calculated in []172(|. This is one of very few physical 



problems we know of where the hard thermal loop resummation technique [171| should be 
used (and was used) at full length. The spectrum feels the effects due to quark and gluon 
interactions in the region ^ Pi+i- ~ gT — the spectrum at larger energies and momenta 

is the same as for the gas of free quarks. One particular source of soft dileptons is the process 
<t\_ ^ + I'^l^ . The probability of this process has a "spike" for the momentum of q*^ 
and q*^ coinciding with the momentum k* on the longitudinal plasmino dispersive curve with 
zero group velocity. There are just many plasminos at the vicinity of this point and the phase 
space factor provides a singularity at Ei+i- = uj±{k*) —uj^^{k*) in the spectrum. Another spike 
comes from the process when a longitudinal plasmino with momentum k* annihilates with a 
longitudinal antiplasmino with the opposite momentum to produce a lepton pair with zero net 
momentum and the energy 2LV^^(k*). 

Unfortunately, in the soft region, the main contribution in the spectrum is due to cuts. 
In other words, the most relevant elementary kinetic processes are not q* ^ q* + l^l^ or 
q* + q* ^ I'^l' , but rather q* + g* q* + g* + l^l^ etc. The spikes actually have finite 
width due to collisional damping of collective excitations (the issue to be discussed in the next 
section), and one can hope to see only a tiny resonance on a huge background. Still, such a 
resonance in the spectrum is an observable effect. 

6.4 Damping mayhem and transport paradise. 

Direct decay. 

In the previous section, we discussed the Landau damping contribution to the imaginary parts 
of polarization operators and, correspondingly, to imaginary parts of dispersion laws. It comes 
from the kinematic region uj < k and is physically related to absorption of ingoing excitations 
by thermal quanta like in (6.57). However, we did not say a word about the contribution of 



direct decay processes g* q + q etc. in the timelike kinematic region to > k. 

That was with a good reason. On the one-loop level, the contribution of decay processes 
in imaginary parts is nonzero and is of order ~ g^T. Unfortunately, it depends on the gauge 
and, in some gauges, has even the wrong sign corresponding not to damping of excitations but 
to instabilities |173]. The point is that such one- loop calculation is unstable with respect to 



higher-order corrections. It is very clear physically — quarks and gluons in QGP cannot be 
treated as massless but acquire dynamical masses due to thermal effects. And the decay of a 
plasmon or plasmino into two other collective excitations is not kinematically allowed. The 
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only exception is the process g* ^ q* + q* which in principle may go if |17C] 

Nf > 9cF - 2cy = 6 (6.58) 

But there are at most three light flavors in real QGP and decay processes can be safely 
forgotten. 

Collisional damping. 

Still, damping is there even in the timelike region uj > k due to collisions g* + q* g* + q*, 
q* + q* ^ g* + g* etc. This is also the main source of damping of transverse electromagnetic 



waves in usual plasma [22|. A rough estimate for collisional damping in QGP can be done 
very simply. 

The meaning of damping is the inverse lifetime of excitations. We have 

C ~ (riife)-' ~ ncT*"* (6.59) 

where n ~ is the density of the medium and o"*"* for excitations which carry (color) charge 
has a Coulomb form 



-a'i'^-T- (6-60) 



We took into account the fact that the power infrared divergence for the integral of Coulomb 
cross section is effectively cut off at p_\_ ~ gT due to Debye screening ^ As a result, we obtain 
the estimate 

C'^ ~ g^T (6.61) 

Note that this value for the damping is unusually large. It is much larger than, say, the 
damping of photons in ultrarelativistic e'^e~ - plasma. The latter can also be estimated from 
the formula ( 6.59| ), but 0"*°* is now not the Coulomb, but the Compton cross section. The 
integral has now the form ^f.xdp\/p\ ~ ln(l/e) and the estimate is 

C ~ e^rin(l/e) < e^T (6.62) 

The question arises whether the new anomalously large scale ~ g'^T has a physical relevance 
P^ . We will return to discussion of this point a bit later. 



°and due to Landau damping effects — see more detailed discussion below. 

^^Do not confuse this scale with the magnetic scale which is also of order g^T. The former is related to kinetic 
properties of the system while the latter refers exclusively to static phenomena. 
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Figure 29: 



An accurate calculation of the damping of fast moving [k ^ gT) quark and gluon excita- 



tions in QGP has been done in |170| , |174|| (see also |175| ]). Consider the graph in Figj2^ for 
the quark polarization operator where the lines with blobs stand for quark and gluon propa- 
gator dressed by thermal loops. The imaginary part of the loop in Fig.^ depends on the 
imaginary parts of internal propagators. The imaginary part of the gluon propagator due to 
Landau damping turns out to be of paramount importance. Physically, this contribution just 
corresponds to the scattering processes g* + q* g* + q* and g* + g* g* + g* as can be 
easily inferred if spelling out the exact gluon propagator as in FigJS^ (there is also a similar 
graph with internal gluon loop). 



42 T 



It can be shown |170| that when calculating the leading contribution in ( in the kinematic region k 3> gT which 



is under discussion now, vertex corrections can be disregarded. 
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Figure 30: A graph contributing to Im H{k). Crossed lines are put on mass shell. One of the 
quarks in the internal loop has a negative energy. 



Dressing the quark propagator (indicated by the blob in Fig.|30D is important. If not 
taking it into account, the imaginary part of the propagator is (5 - function and the result for 
Im 'Eji{lo, k) near the mass sheh uj = |k| has the form 

lm^R{k) = --^ujWcfT-/'^ f dpo [ dcos9 [ ^ 6{po - pcosO) (6.63) 

O^vr^ J-oo J-l Jo rfi^ _|_ J.,,,4 Pjl 

where Q is the angle between k and p. The integral has a power infrared behavior at p ^ 
5fT, but this divergence is cut off due to Landau damping effects (the second term in the 
denominator in the integrand). Still, the integral in ( |6.63| ) diverges logarithmically at p ^ gT 
and the result for the damping inferred from Eq.( |6.63|) is infinite. 

The crucial observation is that the dressed quark propagator does not have singularities 
on the real po axis. A self-consistent account of the collisional damping for the quark Green's 
function moves its singularities in the complex plane. As a result, 5 - function in the integrand 
is replaced by a smooth distribution with the width of order ( ~ g^T: 

6{po-pcos9) ^ ^^""0,2 ^,2 (6-64) 

\po — pcosoy + c 

This smoothing cuts off the logarithmic singularity in (|6.63| ) at p ~ g^T. The other source 
for the cutoff could be provided by magnetic screening effects , but the latter is an essentially 
nonabelian phenomenon whereas the cutoff due to smearing out the 6 - function is a universal 
effect which occurs also in an abelian theory. 

Note that the cutoff p ^ g^T corresponds to the cutoff po ~ Q^T. Really, the integral 
in ( |6.63| ) is saturated by the values p and po such that two terms and ~ uj^^^Pq/p^ in the 
denominator in the integrand are of the same order. That means that p^^^"^ ~ p^/Wpj which 
gives the estimate pq ~ g'^T at the lower range of integration. Actually, our calculation 
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corresponds to picking up the pole corresponding to the so called attenuating mode 



of the exact gluon propagator in the graph in Figj29|. This pole presents an extra solution 



of the transverse dispersive equation in Eq.( |6.43 ) when the imaginary part of the logarithm 



in 'n.t{po,p) which is responsible for Landau damping is taken into account, and we assumed 
\Po\ P- The attenuating mode ( 6.65 ) was found in [|170|| . It exists in the region g^T ^ p 'Z gT 



and the dispersion law ( |6.65| ) holds in the region g^T <^ p <^ gT. When p ^ gT, the 
assumption \pq\ <^ p is not fulfilled, and when p ^ g^T, perturbation theory breaks down. A 



similar mode exists also in nonrelativistic plasma (see Ref.|22], p. 135). We will meet the mode 



( |6.65 ) and the frequency scale ~ g^T once again in the end of this chapter when discussing 



chirality non-conservation rate. 

To find the dispersion law, one should add zim S/j(k) from Eq.( |6.63D to the one- loop result 
for (which is real for to > k) and solve the dispersive equation Det||7'^w — 7k— S/j(k)|| = 

0. The solution is complex oj^°^^{k) = uj^ '■°°P{k) — iC,{k) and the final result for C,{k » gT) is 
very simple. 

C'' = asCFT\n{Cyg) 

C^9 = asCvT\n{C3/g) (6.66) 

(damping of gluon fast moving excitations can be calculated along the same lines). Only the 
coefficient of logarithm can be calculated. The constants C, under the logarithm cannot 
be determined. Actually, we will see shortly that these constants are gauge-dependent and 
cannot be defined in a reasonable way. 



Ill 





Figure 31: Polarization operator of soft plasmino. 



The results ( ^.66 ) were obtained solving the dispersive equation with the polarization 
operator ( 6.631) calculated at real uo. A more refined analysis which takes into account the 
modification of Im Ti{k) for complex u: when one starts to move from the real u) axis towards a 
singularity shows that the dispersive equation has actually no solutions and the singularity is 



no longer a pole, but a branching point |176| , |177| , |17^] . But this branching point is located at 
the same distance from the real axis as the would-be pole and brings about the same damping 
behavior of the gluon retarded Green's function G^{t) ~ exp{— ^f} at large real times. 



The modification of the fermion propagator due to higher loops in the graphs in Fig,2£, 30 
was taken into account smearing out the 5 - function singularity in the Green's function as in 
Eq.( |6.64] ). In recent |179| it was found that if doing things quite carefully, namely resumming 
leading infrared divergent pieces in the soft photon loops contributing in the fermion Green's 
function in any order in abelian theory, the infrared divergence in ( |6.63| ) eventually survives. 
That essentially means that G^{uj) does not involve singularities at all at finite distance from 
the real axis. This results in a non-exponential decay of G^{t) at large real time: 



G^{t) ~ exp{-art In(ert)} 



(6.67) 



In spite of the fact that this result is rather unexpected, it seems to be correct. Non-standard 
logarithmic factor in ( |6.67 ) resembles the logarithmic factor in the exponent in ( |6.23| ) . An anal- 
ogy can also be drawn with the well-known logarithmic divergence in the phases of Coulomb 
scattering [ p.80| |. A careful independent study of this question is, however, highly desirable. 
Damping of excitations in another kinematic region k = ( "standing" plasmons and plas- 



minos) was studied in | 171 |. Here the vertex corrections are as important as corrections to the 



propagators, however an accurate analysis of |171] shows that it suffices to take into account 
only one (hard thermal) loop corrections both in polarization operators and vertices. Consider 
for definiteness the damping of standing plasminos. The graphs contributing to the soft plas- 



mino polarization operator are shown in Fig. 31. Using Keldysh technique one can derive in 
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the limit of soft external momentum 



(2vr) 



T^{k,k + q-q)G''{k + q)T^{k + q,k- 



(6.68) 



Here are the retarded vertices presenting certain combinations of the Keldysh components 



Tf'' and Tf^"^: 



6,c=l 



Ibcd 



(6.69) 



where the index 1 is put on the outgoing fermion line. Like retarded propagators, are 
analytic in the upper uj half-plane. Substituting here transverse and longitudinal parts of the 
gluon Green's function and one-loop vertices T^u and F^, calculating the integral (which in 
this case can be done only numerically), and solving the dispersive equation, one arrives at 
the result 



(i{k = 0) = 1.43CF- 



However, the gluon Green's function involves also a gauge-dependent part 



D^!^\q) = (a - I) 



[{qo + ie)2 - q2]2 



(6.70) 



(6.71) 



where a is a gauge parameter and an infinitesimal ie is introduced to provide for the right 
analytical properties. At first sight, this gauge-dependent piece should not affect the position 
of the pole. Really, one can use the Ward identities 



q,rj^{k,k + q;q) = ^ {[G^[k + q)]'' - [G^ik]]-'} 
q^,T%{k-q) = e{T^{k,k-q--q)-T^{k + q,k;-q)] 



(6.72) 



which hold also at finite temperature for retarded propagators and vertices [181, 174| to derive 



S^(")(A:) = -ig^Tcpia - l)[G^{k)]-^ x 



_A_1 

(2^)4(Zol[(9o + ^6)2-q2]2 [(g, 



- «ej 



q 
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G''{k + q)-G^{k) [G^-ik)]-' (6.73) 



Ri 



We see the presence of factors [G^{k)\^^ both on the right and on the left. G^{k) is singular at 
the pole and [G^{k)]~^ is zero. One might infer from this that a gauge-dependent contribution 
to S^(pole) and hence to the corresponding solution of the dispersive equation determining 
the pole position is also zero. 

However, this is wrong [182, 177[. The point is that the integral in (|6.73 ) involves a severe 
power infrared divergence and is infinite at the pole. We have thereby a x cxd uncertainty. 
This uncertainty can be resolved by choosing k not exactly at the pole but slightly off mass 
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shell. Then [G^{k)]~^ is not exactly zero and also the divergence in the integral is cut off 
by an off-mass-shellness. When the distance from the mass shell is small, the final result for 
^R(a)^^ _^ pole) does not depend on this distance and is just finite. In the soft momenta 
region 

j,HW„,__!(£LZi^,o (6.T4) 

This brings about a gauge-dependent part ~ g^T in the damping of soft plasminos. A similar 
analysis with the same conclusion can be carried out for plasmons. 

In Ref. [|I83| , another regularization procedure was suggested where the external momentum 
k was always kept strictly on mass shell, but the divergent integral in Eq.( |6.73| ) was regularized 
introducing an explicit infrared cutoff /i. Then, of course, a gauge-dependent piece is absent 
even after the cutoff is eventually sent to zero. It does not solve a problem, however. The 
physical definition of damping (if any) should be related to the exponential asymptotics of 
G^{t) at large real times. G^{t) is given by the Fourier integral over the real uj axis 



/oo 
Gr[uj) e"'"^* duj 
'OO 



Thus it is determined by the behavior of Gp;{uj) at real (jJ, i.e. off mass shell. The damping 
defined as 

C = - l,m ^ 

t—rCO t 



is gauge-dependent ||177]] . 

The same gauge-dependence shows up in the damping of energetic plasmons and plasminos, 
but in the latter case, this gauge-dependence is parametrically overwhelmed by the leading 



gauge-independent contribution (6.66) involving the factor ln(l/5f). 

Observability. 

The observed gauge-dependence of the damping obviously indicates that it is not a physical 
quantity. This is definitely true at least for soft plasmons and plasminos where the gauge- 
dependent part and the gauge-independent part ( 6.7C| ) are of the same order ~ g^T. 



Indeed, it is not possible to contemplate a physical experiment where this quantity could 
be measured. That should be confronted with the case of abelian plasma where damping 
of electromagnetic waves is a perfectly physical quantity and can be directly observed by 
measuring the attenuation of the amplitude of a classical wave with time. But as we already 
noted, no classical gluon or quark waves exist. This observation refers also for the damping of 
electron and positron collective excitations in the ultrarelativistic abelian plasma. It also has 
an anomalously large scale ~ e^T [with the extra logarithmic factor ln(l/e) for k 3> eT] and 
it also cannot be directly measured. 

One could try to observe the effects due to damping in gauge-invariant quantities like the 
polarization operator of electromagnetic currents. An accurate analysis which goes beyond the 
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conventional hard thermal loop resummation technique and effectively resums a set of ladder 
graphs shows, however, that a self-consistent account of the corrections due to damping in the 
quark Green's functions and in the vertices results in the exact cancelation of the anomalously 
large scale ~ g^T in the final answer yi74| . (for a similar analysis with a similar conclusion in 
scalar QED see pj ). 



However, there is a physical problem where the scale ~ g^T can in principle show up. 
This is the already discussed problem of lepton pair production in QGP. We have seen that 
the spectrum of leptons pairs involves spikes associated with a special point with zero group 
velocity on the longitudinal plasmino dispersive curve. Going beyond the hard thermal loop 
approximation and taking into account the effects due to collisional damping in the Green's 
functions and in the vertices would bring about a finite width for these spikes of order g^T 
and there is a principle possibility to measure this width. This problem has not been studied, 
and it is not clear by now whether the width of the spike can be calculated analytically and 
whether one can single out this spike out of the background. 

What one can say quite definitely is that this width should crucially depend both on modifi- 
cation of propagators due to collisional effects and on a similar modification of vertices and has 
little to do with the (gauge-dependent) position of the pole (or whatever the real singularity 



is |176| , \17% |17^] ) of the quark and gluon Green's function. 

Thus we are convinced that the latter is not a physical quantity probably even for energetic 
plasmons in spite of the fact that the leading contribution ( |6.66| ) is gauge-independent there. 
We just do not know how on earth this quantity could be measured. 

Transport phenomena. 

There is a lot of kinetic phenomena in QGP which are physical and measurable. Indeed, 
nothing in principle prevents measuring the electric resistance of a vessel with QGP or studying 
the flow of QGP through narrow tubes. They do not depend, however, on the anomalous 
damping scale g^T, but rather on a much smaller scale 

(w)-'~/Tln(l/g) (6.75) 



This scale already appeared in ( |6.62| ) determining the damping of electromagnetic waves in 
e^e~ - plasma. And it is also the scale which determines the mentioned physical effects of 
viscosity and electric conductivity, and many others — heat conductivity, energy losses of a 
heavy particle moving through plasma, etc. 

The appearance of the scale g'^Tln{l/g) has a clear physical origin. All the mentioned 
effects are inherently related to the rate of relaxation of the system to thermal equilibrium. 
The latter can be estimated as 

{Trei)-' ~ na'^^^' (6.76) 
It looks the same as the estimate for lifetime ( |6.59 ) but with an essential difference — in 



contrast to (6.59), the estimate (6.76) involves the transport rather than the total cross section. 
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The transport cross section is defined as 

Jda{l-cos9) (6.77) 



^trans 



where 9 is the scattering angle. The factor (1 — cos 9) takes care of the fact that smah-angle 
scattering though contributes to the total cross section, does not essentially affect the distri- 
bution functions ng(p), nq{-p) and is not effective in relaxation processes. For the Coulomb 
scattering in ultrarelativistic plasma, the transport cross section is 

trans 4 [ dPj^P± 9^ , 1 fn 7o\ 

J{gT)^ Pi g 

Multiplying it by n ~ and substituting it in ( |6.76| ), the estimate ( |6.75| ) is reproduced. 

Viscosity and all other similar quantities can be calculated analytically in the leading order 
(probably, magnetic infrared divergences prevent an analytic evaluation of these quantities in 
next orders in g, but this question is not yet well studied). It is interesting that Feynman 
diagram technique proves to be technically inconvenient here , and the good old Boltzmann 



kinetic equation is the tool people usually use (see e.g. [185|). 

Let us make two illustrative estimates which make clear how the relaxation scale ( |6.75| ) 
depending on the transport cross section ( |6.78 ) arises. 



First, let us estimate the electric conductivity of QGP (it is the quite conventional conduc- 
tivity, not the "color conductivity" which is sometimes discussed in the literature, depends on 
the anomalous damping scale g'^T, and is not a physically observable quantity — we do not 
have batteries with color charge at our disposal). Suppose at t = the system was at thermal 
equilibrium so that the quark distribution functions are no(p,x, t = 0) = ni?(|p|). When we 
switch on the electric field, the distribution function starts to evolve according to the kinetic 
equation 

9n(p,x,t) an(p,x,i) an(p,x,i) 

with V = p/|p|. Dots in RHS of Eq. (|6.79|) stand for the collision term which becomes relevant 
at t ~ X ~ Trei ~ Tfree path ~ [g'^T ln{l / g)]^^ . Thus the elcctric field brings about distortions 
of the distribution function which grow up to the characteristic value 

V eTEv 

dn ~ ehi—noTfree path 



.^.^u.ireepa.n ^4 1^(1/^) 

At this point, collisional effects stop the growth (a particle drifting in external electric field 
collides with a particle in the medium, forgets what happened before, and starts drifting anew). 
The density of electric current in the medium is 

i = e v6n d?p ~ E (6.80) 
J rln(l/5') 

The coefficient between j and E gives the conductivity. 
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Let us estimate now the energy losses of a heavy energetic quark in QGP. Of course, free 
quarks do not exist, but a physical experimental setup would be sending into the bottle with 
QGP a heavy meson Qq with open beauty or top. In QGP, the meson dissociates, and a 
naked heavy quark propagates losing its energy due to interaction with the medium. It goes 
out then on the other side of the bottle dressed again with light quarks, but not necessarily 
in the same way as before. When Mq ^> Aqcd, this dressing does not essentially affect its 
energy. A heavy particle containing Q can be detected and its energy can be measured. 

Suppose a heavy quark is ultrarelativistic, but its energy is not high enough for the Cerenkov 
radiation processes to be important. Then the energy would be lost mainly due to individual 
incoherent scatterings. The mean energy loss in each scattering is AE ^ T [T — is a char- 
acteristic energy of the particles in heat bath on which our heavy quark scatters). The mean 
time interval between scatterings is Tfree path ~ [d'^T hi{l / g)]^"^ . We obtain 

-^~^^ = Ca^r2ln(l/5) (6.81) 

(IX Tiree path 

This estimate turns out to be correct up to the argument of the logarithm which in reality is 
energy-dependent 186 |. The numerical coefficient C = Attcf/S was determined by Bjorken 

m 



6.5 Chirality drift. 

Up to now, we were concentrated mainly on the perturbative calculations in QGP. Also 
we discussed the magnetic screening phenomenon which, as was explained in Sect. 6.1, is 
essentially non- perturbative. Magnetic screening scale ~ g^T shows up in static correlators of 
different quantities at large distances. Magnetic screening also affects kinetic characteristics of 
QGP such as plasmon and plasmino dispersion laws, transport coefficients etc. The latter are 
well- defined and analytically calculable in the lowest perturbative order. The problems appear 
(and non-perturbative physics come into play) only when trying to determine higher-order 
corrections. There is, however, a beautiful kinetic effect which is purely non-perturbative. 
This is the axial charge non- conservation in hot quark-gluon plasma, or the "chirality drift" 
presenting the subject of this conclusive section. 

Theoretical picture. 

The primary practical interest of this quantity is related not to QGD, but to the dynamics 
of electroweak theory in early hot Universe. Long time ago, 't Hooft observed |145| ] that the 
baryon charge is conserved in standard model only on the classical lagrangian level. It is not 
conserved in the full quantum theory due to anomaly 

5^j;f = ^G%G% (6.82) 

where is the SU{2) electroweak field. Baryon number non-conservation is precipitated 
by the gauge field configurations with nonzero Pontryagin number, the instantons. At zero 
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temperature, the action of electroweak instantons is very large ~ I'k ja^ and the rate of baryon 
number non-conservation involves the exponential factor ~ exp{— 27r/aiy} which is zero in all 
practical sense. 

Physically, instantons present tunneling trajectories through a barrier in the functional 
space. This barrier is associated with a collective degree of freedom called Chern-Simons 
number 

^2 



Q 
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9_^abc 
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£_ ^abc Ji^ 



(6.83) 



The integrand in (6.83) is not gauge invariant, but Q{t) itself is invariant under topologi- 
cally trivial gauge transformations. The classical energy functional in Yang-Mills theory has 



degenerate minima at integer Q (see Fig. 32). These minima are related to each other by 
topologically non-trivial gauge transformations. In electroweak theory where gauge symmetry 
is broken spontaneously by the Higgs mechanism ^ , high barrier between adjacent minima 
is created. The top of the barrier between, say, the minima with (5 = and Q = 1 (actu- 
ally, it is a saddle point in the functional space with only one unstable mode) is called the 
sphaleron ]188| ] and has the energy .Bsph ~ mw / aw ■ Chern-Simons number of the sphaleron 
configuration is Qsph = 1/2. 



May be, this generally adopted terminology is not the best. In some sense, gauge symmetry is never broken 
as the generator of gauge transformations , the Gauss law constraint, still gives zero acting on the vacuum state, 
and there are no massless Goldstone bosons. However, better words were not invented and, anyway, everybody 
understands what the Higgs mechanism is. 
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Ncs=-1 Ncs=0 Ncs=l Ncs =2 

A=0 A=uVu"^ 



Figure 32: A schematic representation of the (bosonic) potential energy along a particular direction 
(labeled z) in field space, corresponding to topologically non-trivial transitions between vacua. 



It was observed in |189| that, when the temperature grows, the height of the barrier de- 
creases. Baryon number non-conservation processes are thus greatly facilitated. Actually, as 
soon as T ^ mw , the leading contribution in the rate comes not from quantum tunneling pro- 
cesses "through" the barrier, but from the classical thermal jumps "over" the barrier. These 
classical transitions are stochastic in nature and can occur both with increasing and decreasing 
of baryon number. We have the physical picture of Brownian motion. Baryon charge "drifts" 
so that the average square of the fluctuation of the baryon charge [B{t) — 5(0)]^ grows linearly 
with time. Due to the algebraic identity 



Q = ^1 d'^Gl^Gl, = B, (6.84) 

one can study alternatively the fluctuations of Chern- Simons number Q{t). The rate of baryon 
number non-conservation T'^^^^{T) is defined as the coefficient of the large time asymptotics 
of the thermal correlator 

c{t) = < [Q{t) - Qio)f >T r^^^O(r)yt, (e.ss) 

where V is the spatial volume of the system (naturally, the larger is the volume, the larger 
is the net rate), and we subtracted the time-independent part. In the temperature region 
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mw{T) < T « £^sph(r), the estimate for the rate is T^^f^°{T) oc exp{-EsphiT)/T}. At 
some temperature (which is close to the critical temperature Tc where the spontaneously 
broken gauge symmetry SU{2) U{1) is restored), the characteristic energy of the excitations 
is of the same order as the barrier height at which point the exponential suppression disappears. 
An accurate calculation of the rate in the temperature region where the suppression is still 
there was performed in |190|] . 

The question arises what happens at large temperatures T ^ Tc where no barrier exists and 



the quasi-classical approximation breaks down. It was argued in [19C] that the rate behaves 
as 

r^^^°(T) = K{awT)^ (6.86) 



where k is a constant which cannot be calculated analytically. The rate ( 6.86| ) is a purely 



non-perturbative quantity. One can show that the rate r^-^^°(T) defined in ( |6^ IS zero m 
any order of perturbation theory. |^ The rate ( |6.86| ) is large, it leads to the dissipation of 
any primordial baryon charge density (if i? — L is exactly conserved as it does in most models 
of grand unification under discussion) and makes our very existence problematic. I will not 
pursue the discussion of this issue further here and refer the reader to a comprehensive review 

The estimate (|6^ ) can be obtained as follows. As was mentioned, the barrier does not 
exist at high T. There are field configurations A(x) with arbitrary small energy which have 
(5 = 1/2 and interpolate between adjacent minima. Such configurations should have, however, 
large spatial size: the larger is the size, the smaller are characteristic field gradients and the 
smaller is the energy E J d:x.{dA)'^ . ^ One can visualize these configurations as the slices 
r = of the standard BPST instantons in the hamiltonian gauge = 0. It is easy to see 
that the energy of such a slice for the instanton of size p is E{p) ~ l/{awp)- Writing the 
integral over the collective coordinate p with the standard zero temperature instanton measure 
~ / dp/p^ (see Ref. |l95|] ), we obtain 



r-^»(T) ~ /^»p{-^} ~ («T)^ (0.87) 

The characteristic p where the integral is saturated is nothing else as our old friend, the mag- 
netic screening scale. Actually, one can forget at this point, if one wishes, about electroweak 



^"^That does not mean that the correlator C{t) itself is zero in perturbation theory. It is not zero, diverges in 
ultraviolet as a power, and not only in nonabelian theory, but also in standard QED. Already the simplest graph 
with a photon loop gives a nonzero contribution ~ a^T^Aj/yV^ in C{t) which is physically related to the fact that 
a circularly polarized photon has nonzero Chern-Simons number density. These perturbative contributions do not 



grow, however, at large t and do not affect the rate ^9l\ |192 |. 

^^In other words, conformal symmetry of the classical Yang-Mills lagrangian does not allow soliton and sphaleron 
solutions. The latter exist when conformal symmetry is broken and a characteristic mass scale appears. In electroweak 
theory, this scale is provided by the Higgs expectation value. An interesting model is Yang-Mills theory on a spatial 



sphere of finite radius where the latter brings about the scale and sphaleron solutions appear |194| 
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theory altogether. In the high temperature region, it presents "electroweak plasma", and its 
physics is basically the same as the physics of QGP which is of main interest for us here. 
Instead of baryon number non-conservation, we have here the axial charge non-conservation. 
Like baryon charge, it is conserved on the classical lagrangian level in the theory with massless 
quarks and is not conserved in the full quantum theory due to Ua{1) anomaly. The rate of 
axial charge non-conservation in hot QCD can be determined from the asymptotics of the 
Chern-Simons correlator ( |6.85 ) at large Minkowski times. 



The result ( |6.86| ) means that roughly one transition occurs in the volume ~ {g^T)~^ by 
the time ~ (g'^T)^^. The estimate for the characteristic volume seems to be rather solid: 
indeed, the barrier becomes irrelevant and also perturbation theory breaks down at the scale 
p ~ {g^T)~^ . However, the estimate for the characteristic time of the transition t^har ~ {g^T)^^ 
was recently criticized in [196]. The time t ~ [g^T]"^ is the characteristic time of smearing 



out the wave packet with the size p ~ {g'^T) ^ in empty space. But our space is not empty, 
it presents a heat bath densely packed with gluon and quark excitations. It was argued in 



Ref.| 196 | that the time it takes for a quasi-sphaleron to decay in quark-gluon plasma plasma 



is actually suppressed by two powers of g and is estimated as 

Way ~ i/T)-^ (6.88) 

Correspondingly, the new estimate for the rate is 

^ ^ ~ a^T^ (6.89) 



Pchar *decay 



The estimate (|6.88|) and its corrolary ( 6.89| ) can be obtained as follows. Suppose, we have at 



t = a wave packet with the size ~ {g^T)~^. Expand it in Fourier series and solve a classical 
equation of motion (which takes into account the heat bath effects) for an individual mode 
with initial condition 

Afc = (6.90) 

The solution presents a sum of three exponentials 

Ak{t) = aiexj){—Cig'^Tt — iuj±{k)t)} + a2exp{—C2g'^Tt — iuj^\{k)t)} 

+f3exp{-Csg^Tt} (6.91) 

The first two terms correspond to standard transverse and longitudinal plasmons and involve 
a plasmon decay rate factor with ( ^ g'^T while the third one corresponds to attenuating mode 
( |6.65 ) with the decay rate ~ g'^T at A; ~ g^T. The first and the second exponential come with 



small coefficient ~ g^ due to the boundary conditions ( 6.90| ) and will die away at the time 



scale ~ {g'^T)~^ while the second exponential is still alive. It would eventually determine the 
quasi-sphaleron decay rate. 

In our own opinion, these arguments are not conclusive enough. Let us first modify the 
initial conditions and put ~ g'^TA^- Such initial conditions seem to be as reasonable 
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as ( |6.90 ). But then the coefficients ai^2 and f3 would be of the same order and the Chern- 



Simons number of the field configuration ( |6.91 ) would be essentially modified on the time scale 



{g^T)^^ rather than (gi^T)^^. Second, the amplitude of Fourier harmonics in quasi-sphaleron 
wave packet is large and their nonlinear interaction is essential. Even if we started up with 
only attenuating mode at the initial moment, normal plasmon modes would be excited due to 
nonlinear effects and would eventually determine the time decay scale. ^ 

Another question which can be asked here is why we assumed that the plasmon harmonics 
die away with the anomalous damping scale ~ g^T rather than with the scale ~ g'^T which is 
characteristic for transport phenomena in plasma and which, as we repeatedly argued in the 
previous sections, is more physical. The point is that the latter scale depends on the transport 



cross section ( 6.77 ) where the region of small momentum transfer ~ gT is suppressed. In 
this particular problem, we do not see a reason to include such a suppression, however. The 
scale gT is large compared to the characteristic momentum scale g'^T determining chirality 
drift. When a mode acquires the momentum ~ 5T, it is lost as far as the coherent quasi- 
sphaleron wave packet is concerned. 

But all these arguments are extremely heuristic. To solve the problem accurately, it does 
not suffice to play around with small fiuctuations on the vacuum background. Actually, one 
has to expand the field near the quasi-sphaleron background A = A;, + (5A and to solve the 
linearized equations for (5A in the heat bath. 

Let us first see what happens in the absence of heat bath. Vast majority of the modes 
are stable. There are two directions in the functional space decreasing the energy of the field 
configuration. One of them corresponds to rescaling the quasi-sphaleron configuration and can 



be handled introducing the collective coordinate p and integrating over it as in (6.87). The 
mode of main interest for us is the intrinsic sphaleron unstable mode SKunst which corresponds 
to rolling down the slope away from the quasi-sphaleron ridge. In vacuum, the corresponding 
linearized equation is 5 A ~ [A;^[^j^j.]5A with fcchar ~ Q^T. It describes the decay of quasi- 
sphaleron by the time ~ {g'^T)~^. 

Let us now switch on the medium effects. The problem has not been accurately solved 
yet, but, taking insight from the dispersive equation (|6.43| ) on the vacuum background, we can 
tentatively add to /c^har transverse polarization operator nt(u;, fechar)- The new dispersive 
equation looks something like 

-U>^ = fc,\,,+ni(u;,/Cehar) (6-92) 

It differs from Eq.( |6.43| ) only by the sign of lj^ term. Substitute now here Ilt{uj , kcha.r) ~ 
— iWpiw/Zcchar which is true at small uj <C A:char- The solution to the dispersive equation ( 6.92| ) 



^^As B. MuUer pointed out to me, there is also another effect of losing of coherence (or spread) of the quasi- 
sphaleron wave pocket due to dispersion. But a characteristic time scale of this effect is estimated as (rgprcad)^^ ~ 
{dbj-pil dk)kc\iax ~ g'^T . This time is larger than (g^T)^^ so that the all harmonics in the wave packet would die 
away before the coherence is lost. 
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IS 

^ ^ ig^T (6.93) 

This mode corresponds to damping of the fluctuations 6 A with time so that a perturbed field 
tends asymptotically back to the quasi-sphaleron configuration A^. This mode has nothing 
to do with instability. 

The dispersive equation ( |6.92| ) has also an unstable solution, but it corresponds to large 
UJ » fcchar where Ilt{uj , kcha.r) ~ i^pi- The corresponding unstable mode is 

Wunst ~ iujpi ~ igT (6.94) 

Thus medium increases the rate of the quasi-sphaleron decay rather than decreases it. The 
estimate is 

r ~ (6.95) 



This reasoning is quite parallel to the reasoning of Ref.|196|, only the sign of nj(a;, /cchar) in 



the dispersive equation ( |6.92 ) is chosen positive rather than negative. Negative sign would 



transform the damping mode ( |6.93| ) into an unstable one, and the estimate ( 6.891 ) would be 
reproduced. 

By no means we want to insist here that it is the estimate (6.95) rather than ( |6.89| ) which is 



correct. Choosing the positive sign for Ilt(uj,k) in ( |6.92| ) looks more natural, but the thermal 
polarization operator corresponding to fluctuations on a quasi-sphaleron background (with 
the Green's functions on a quasi-sphaleron background in the loops) could essentially differ 
from the vacuum one, and only its accurate calculation would solve the problem. It is possible 
that the the rate is suppressed compared to the estimate ( |6.86D , and it is also possible that it 
is enhanced. ^ The third possibility (which is actually the most attractive) is that plasma 
effects neither increase nor decrease the rate of decay, but leave it basically intact. 

Numerical experiment. Kinetic equations. 

Let us discuss now experimental (lattice) data. A very serious difficulty here is that the rate 
of baryon number (chiral charge) non- conservation is a kinetic quantity defined in real time. 
We are not able currently to calculate numerically field theory path integrals in Minkowski 
space- -time and have to look for some other way to solve the problem. 

As was mentioned earlier, the physical mechanism of chirality non-conservation is the clas- 
sical drift over the barrier. The idea arises to obtain the result just solving the classical equa- 
tions of motion and looking at how the classical Chern-Simons charge ( |6.83| ) is changed with 



^^We should note, however, that the enhancement of the rate as is suggested by the estimate (6.95) does not look 



physically appealing. It impUes that the object with the size ^ [g^T) ^ dissipates faster than it takes light to travel 



it across. Basically, we are presenting this estimate just as an illustration that the same reasoning as in |19(;| with 
only a little natural modification leads to completely different results. 
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time. Performing an average over initial conditions with the classical weight ~ exp{— /Jff^'}, 
one estimates the rate of the drift in hot field theory. This algorithm was suggested in 1197]] 



and successfully tested in some 2~dimensional models |197, 198 |. In [191| and in recent | 192 ], 
the same algorithm was applied to 4-dimensional gauge theories. To discretize the classical 
equations of motion, the system was put on a 3-dimensional hamiltonian lattice, the lattice 
spacing a playing the role of the ultraviolet cutoff. 

The question whether such a classical algorithm is justified or not is under discussion now. 
The main problem here are severe power ultraviolet divergences inherent for a classical field 
theory at finite T. A classical Rayleigh-Jeans ultraviolet divergence in the energy density of 
the photon gas 

E'' ~ / d^p\p\ ni{\p\) tJcI^p TAfjy (6.96) 

[where the classical limit of the Bose distribution is ra^(e) = l/(/5e)] was the reason to develop 
the quantum theory in the first place. Classical contribution to the plasma frequency 

d^p 1 
]P\P\P\ 
is also ultraviolet divergent. 

In the full theory, the divergence is cut off at |p| ~ T (where quantum corrections in the 
Bose distribution function come into play). For the classical approximation to be justified, one 
should always keep ~ Auy <C T. 

In principle, chirality non-conservation rate calculated by the classical algorithm could also 
be cutoff-dependent. Moreover, if the estimate ( |6.89|) of Ref. | 196 | or the estimate ( 6.95| ) is 



(a;^l)2 I ^-±- g^TAuv (6.97) 



correct, one should expect cutoff-dependence of T. In both cases, F depends explicitly on 
the plasma frequency scale which is cutoff-dependent when calculated classically. One easily 
estimates 

^ -pel ^char 9 



Tasy ^ g T ^ Fasy ~ / ci\2 a 



r(6.95) ~ g'T^ =^ rta ~ k^^^co'^, ~ /t^/^aJ/; (6.98) 



However, if the original estimate ( |6.86| ) is not affected by the medium effects, F does not 
depend on the plasma frequency scale, but only on the magnetic mass scale g^T. In this case, 
one could expect that classical numerical results for F do not depend on cutoff, and the whole 
algorithm is trustworthy if the condition 

g'^T Auv 

is satisfied. 

Following the suggestion of |199|], the question of whether classical rate depends on cutoff 



or not, was studied in [192|. Ambjorn and Krasnitz find that the rate depends neither on 
cutoff nor on LOp\. Their answer for the rate is 

T^B^^{T) « l.l{awTf (6.99) 
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This result directly contradicts the estimate ( |6.89 ) and also ( 6.95| ). It was suggested in ]196| ] 
that the result ( |6.99| ) is actually a lattice artifact related to an imperfect lattice definition 
of the Chern-Simons number. It may be, though it is not quite clear why an alleged lattice 
artifact ( |6.99| ) does not depend on the lattice spacing. 

On pure esthetic grounds, it would be much nicer if the plasma effects would not eventu- 
ally influence the quasi-sphaleron decay rate and the original result ( 6.86 ) and the numerical 
calculations of [192] were correct. Our heuristic arguments which take into account not only 
the attenuating but also the plasmon modes and nonlinear interaction between them [see the 
discussion after Eq.( |6.9l|) ] also seem to indicate that. We would rather lay our bets on this 
possibility, but a scientific answer to this question is not yet found. Thus we are in a position 
to discuss what happens if an unfortunate possibility that the rate depends on the scale Wpi 
and thereby on the ultraviolet cutoff, when calculated classically, is realized. Can one still 
hope to suggest a numerical procedure where the rate of the chirality drift in the full quantum 
theory could be evaluated quantitatively ? 

This question was studied in |199[ (see also recent [ ^OOp . An algorithm based on the 
separation of low momentum and high momentum scale was suggested. The idea was to take 
into account explicitly the contribution of high momentum modes in the effective hamiltonian 
and then to solve numerically the classical equations of motion of this effective hamiltonian 
for soft modes. 

The effective theory for soft modes (with momenta of order gT) was first formulated in 
1 171, 181 1 in lagrangian form. This lagrangian involves multiple quark and gluon effective 
vertices of the kind shown in Fig. 33. Characteristic momenta flowing in the hard thermal loop 



in Fig. 33 are of order T. We already discussed at length two-point quark and gluon Green's 
functions which are contained in this lagrangian [see Eqs.daiiD, (|6l9|)]. Both 2-point and 
also 3-point, 4-point etc. effective vertices are nonlocal. 
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p 

Figure 33: An example of HTL graph contributing to the effective 3-gluon vertex. The loop 
momentum is hard p ~ T. External momenta are soft. 
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A nonlocal lagrangian is very inconvenient for a numerical study. Fortunately, it can be 



presented in a local form [201, |165f| . The price is introducing new variables depending on 
X and on an extra 3-dimensional unit vector v. The effective HTL hamiltonian of pure 
gluodynamics can be written in the form 



(6.100) 



where w{x,v) = w'^{x,v)t"- and ujpi = gT^/Wc/?). The variables w"'{x,^) talk to gauge fields 
by virtue of nontrivial equal time commutators 



S,''(x),?i;^(x',v) = i5'^%i5(x - x') 



w;'^(x,v),u;*(x',v') 



Stt 



(ir^'w' - ViVf) 5(x - x')<^(v - v') (6.101) 



{"Df^ is the covariant derivative) and of the Gauss law constraint. Integration over Yl dw{x, v) 
restores original nonlocalities. 

Actually, as was shown in |165], the variables w have a transparent physical meaning. 
They can be related to the deviation of the phase space distribution function of hard gluons 
(5n(x,p) = (5n°(x,p)t° from the thermal equilibrium Bose distribution no(|p|): 



w X 



(5n(x,p)|p|2 c^IpI 



(6.102) 



The equations of motion corresponding to the hamiltonian ( |6.100 ) are equivalent to a non- 
abelian analog of the Vlasov equation system 



V,F^^^ = 2gN,J ^v''6n{x,p 



(6.103) 



v''V^5n{x,p) 



-g^rE 



^^o(|p|) 
d|p| 



(6.104) 



where = (l,v), F^^'^ is the strength of mean gluon field (which is soft), and 5n{x,p) are 
degrees of freedom associated with hard particles. We already met the abelian analog of 
Eq.( |6.10^ ) in Sect. 6.4, but note the appearance of the covariant derivative : in contrast 
to the abelian case, 5n{x, p) carries here color charge. The RHS of Eq. ( |6.103| ) is the induced 
color current of hard particles. 

When also fermions are present, the effective theory can also be formulated in kinetic 



equation language | 203 , 165 ]. To this end, one has to introduce mean fermion soft fields 
and also mixed boson-fermion densities ~< A^il^ >. Such densities and the corresponding 



^^Then the effective action acquires an elegant mathematical structure, being closely related to the 2D Chern- 



Simons theory in light cone variables x± — t ± xv |202| 
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kinetic equations were first introduced and analyzed in |114| where they were used to study 
the problem of goldstino (= phonino) dispersion law in a supersymmetric thermal medium. 

However, as far as the problem of constructing a numerical algorithm to calculate diffusion 
rate in a possible case when classical equations of motion for the tree hamiltonian produce an 
ultraviolet cutoff dependent result is concerned, a simple use of the kinetic equations system 
( 6.103|) , (|6.104|) instead of the classical field equations is not yet a solution. 

The hamiltonian ( 6.100 ) and the equation system ( 6.103 ), ( 6.104 ) were derived in the 
assumption that integration in the hard thermal loop in Fig. 33 etc. goes over all momenta, 
both hard and soft. If we just put the system ( |6.103| ), (6.104) on the lattice with an explicit 
ultraviolet cutoff Auv = soft modes would be counted twice: once in the loop and once 
again in the classical dynamics, //the results of the simplistic classical calculation depend on 
the lattice spacing, so the results of the calculation with the kinetic equations ( |6.103| ), ( |6.104 ) 
would. 

One can hope to obtain a cutoff-independent result only when separating carefully soft and 
hard modes and including only hard ones (with momenta exceeding the separation scale) in the 
quantum loops [ [199| ] . Unfortunately, we were not able to present an explicit and consistent pro- 
cedure of such a scale separation. A simplistic momentum cutoff breaks gauge-invar iance. The 
lattice cutoff preserves gauge invariance but brings about a lot of new rotationally-noninvariant 
structures in the effective hamiltonian which we could not handle. Q 

It is quite obvious that more theoretical and numerical studies of this question are highly 
desirable. 
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On the other hand, i/the classical calculations in the spirit of [ 191 , 192 eventually prove to be cutoff- independent, 



they give the correct result, and all the stuff of effective lagrangian for soft modes and Vlasov equations would play 



here a role of "purple hands on enamelled wall" |204] — very beautiful, but practically irrelevant. 

^°The last technical comment is that if one eventually succeeds in constructing a self-consistent algorithm of scale 
separation, he would simultaneously solve another problem which wc did not mention yet. HTL approximation for 
the effective lagrangian is justified when a characteristic momenta of soft modes is of order gT . In this case, however, 
we need to study the dynamics of modes with momenta of order g^T . In this region, HTL approximation breaks 
down, and one has to resum higher-loop ladder graphs which are all of the same order. In a simple case (for the 



polarization operator of electromagnetic current in abelian theory), such a resummation was done in |174|, but to 
do it for the full effective lagrangian is a formidable problem which is far from being solved now. However, as was 
noted in |19S|| , the dominance of one loop HTL graphs is restored if taking into account only the momenta p ^ gT 



in the loops. Thus, if the separation scale lies in the region gT ^ /i ^ T, HTL approximation is justified. 
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